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Uber das Verfahren von Baranow und seine Erweiterung auf die Ermittlung 
der Kigenschwingungsformen von Schwingerketten 


Von H. Schaefer 


1. Einfiihrung. Das Verfahren von Baranow fiir die Ermittlung der Torsionseigenfrequenzen 
einer Maschinenwelle findet man in Klotters groBem zusammenfassenden Bericht! ausfiihrlich be- 
sprochen. Da ein Beweis fiir die Exaktheit dieses Verfahrens nicht bekannt geworden war (1949) 
wird es dort mit Vorsicht als Naherungsverfahren von erstaunlich guter Genauigkeit gewertet. Dab 
aber der Grundgedanke des Baranowschen Verfahrens, namlich die im nachsten Abschnitt zu be- 
sprechende Reduktion, durchaus als exakt angesprochen werden muf, hat 1952 Haug” betont. Er 
teilt die ersten Schritte eines sehr umfangreichen Beweises mit, dessen Ausgestaltung fiir eine spatere 


Veréffentlichung geplant ist. In einer durch die Lektiire des Haugschen Buches angeregten Arbeit? 
konnte ich 1953 einen anderen, 


vollstandigen Beweis fiir die Exakt- 
heit der Baranowschen Reduktion 
geben. Ferner wurde darauf auf- 
merksam gemacht, daB man den 
Baranowschen Konstruktionen 
mit nur geringer Mihe samtliche 
Eigenschwingungsformen entneh- 
men kann. Zweck der _ vor- 
liegenden Arbeit soll sein, diese 
Erweiterung des Baranowschen 
Verfahrens, die ich seit zwei Jah- 
ren in Ubungen habe erproben 
- lassen, bekannt zu machen. Dabei 
lieB sich nicht vermeiden, eine 
ausfiihrliche Erérterung der Bara- 
nowschen Reduktion voranzustel- 
len. Der hier gefiihrte Beweis 
ist neu und, wie ich glaube, kiirzer 
und anschaulicher als derjenige 
meiner friiheren Arbeit, in der das 
Hauptgewicht auf eine mathe- 
matische Klarung dieser eigen- 
artigen Reduktion gelegt wurde. 


2. Die Baranowsche Reduktion. 
Abb. 1 zeigt eine unverzweigte 
Schwingerkette, das iibliche Er- 
satzsystem einer Maschinenwelle, 
mit den Bezeichnungen, die wir 
in den nachstehenden Rechnun- 
gen verwenden wollen. Die Welle 
ist mit n+ 1 Drehmassen be- 
setzt, deren Tragheitsmomente 


Op, O,,-» O, sind. Die einzelnen 


1 K. Klotter, Ing.-Arch. 17 (1949), S. 49; dort Schrifttumshinweis auf Seite 61 unter [17]. 
2 K, Haug, Die Drehschwingungen in Kolbenmaschinen. Berlin, Gottingen, Heidelberg 1952. 
8 H. Schaefer, Abh. Braunschweig. Wiss. Gesellschaft 5 (1953), S. 141. 


Abb. 4. Die neue, nach der Baranowschen Reduktion entstandene Welle. 
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Wellenabschnitte mit den Lingen 1,, I,,..., 1, seien alle auf dieselbe Torsionssteifigkeit Gl7 = L/y 
reduziert. Die GréBen @, Yj, +++» Pn Sind die Winkelamplituden der Drehmassen, M,, M,, vss M,, die 
Amplituden der Drillmomente in den einzelnen Wellenabschnitten. Denken wir uns die Welle 
links von @, und rechts von @, jeweils ein Stiick fortgesetzt, so sind die Drillmomente M, 
bzw. M,,.1 in diesen freien Enden staindig Null. Endlich sei w die Kreisfrequenz einer Eigen- 
schwingung. 

Wir setzen nun voraus, da® die Eigenschwingungsform héchster Ordnung bekannt ist. Ihre 
Winkel- und Momentenamplituden bezeichnen wir mit gf bzw. Mp, ihre Eigenfrequenz mit po. 
Zwischen zwei benachbarten Drehmassen @,_; und GQ, besitzt die Eigenschwingung héchster Ord- 
nung genau einen Knoten N,, der das Wellenstiick h, in ol, und (1 — ax) h, teilt (Abb. 2). Nun 
sei die Welle in bekannter Weise in 2 n Elementarsysteme aufgeteilt, von denen jedes mit der 
Frequenz @, schwingt (Abb. 3). Die hierzu bendtigte Aufteilung der Tragheitsmomente O,, in By Ox 
und (1 — fx) Ox, mit By = 0 und f, = 1, wie sie Abb. 3 zeigt, gewinnt man aus den Gleichungen 


sat = Ooh +++» = BO, (1 —oy) b, = (1 — Br) On Ogi keys | (2.1) 
= Dent On41(1 — On+1) T4-1 SO O;, (l — On) pe 


Ohne weiteres ist einzusehen, daS man die Winkel- und Momentenamplituden dieser Elementar- 
systeme derart wahlen kann, dafs benachbarte Drehmassen fi, O, und (1 — f;) O, mit derselben 
Amplitude schwingen und dafi gleichzeitig in benachbarten Wellenstiicken a, J, und (1 — «,) i 
dasselbe Drillmoment vorhanden ist. Fast trivial ist demnach die Bemerkung, da man mit diesen 
2n Elementarsystemen die Eigenschwingungsform héchster Ordnung darstellen kann, so wie sie 
in Abb. 3 eingezeichnet ist. 

Wie aber hat man dieselben Elementarsysteme zu benutzen, wenn irgend eine andere Eigen- 
schwingungsform der Welle dargestellt werden soll? So absonderlich die Fragestellung zunachst 
scheinen mag, ihre Beantwortung fiihrt auf die dynamische Deutung der Baranowschen Reduktion. 

Offensichtlich diirfen jetzt zwei benachbarte Wellenstiicke 0; J, und (1 — o,) J, nicht mehr in 
N, festgehalten sein, vielmehr muf die Méglichkeit zugestanden werden, daB beide dort dieselbe 
Winkelamplitude yg_j, die im allgemeinen von Null verschieden sein wird, besitzen. Zwei benach- 
barte Drehmassen f;, 0, und (1 — f;) 0; miissen durch ein zwischen ihnen wirkendes Krafte- 
paar Mf (Reaktionskraft) gezwungen werden, mit gleichen Amplituden zu schwingen, oder anders 
ausgedriickt, die starre Verbindung der benachbarten Drehmassen f{;, 0, und (1 — f;) O, wird 
durch ein Kraftepaar Mf beansprucht. Saimtliche Mf und qj sind Null bei der Eigenschwingungs- 
form héchster Ordnung, bei einer anderen Eigenschwingungsform dagegen von Null verschieden. 
Die geometrischen Bedingungen, daB f, 0), und (1 — f;) 0, dieselbe Amplitude ~,, und die dyna- 
mischen Bedingungen, daB die Wellenstiicke o; |, und (1 —a,) I, dasselbe Drillmoment M;, be- 

g sitzen sollen, liefern die erforderlichen Gleichungen zur Bestimmung 
der Mj und gj. Diese Gedankengange werden in den folgenden Rech- 
nungen ausgefihrt. 
ee a ae Eine elementare Grundbetrachtung sei vorweggenommen. Der 

Schwinger Abb. 5 wird links durch das Kraftepaar \t* = M* coswt 
erregt, das rechte Ende der Welle wird zwangslaufig gefiihrt mit dem 

@ vorgeschriebenen Drehwinkel ©* = ~* cos wt. Gefragt ist nach dem 

Abb. sas het 2 pb al Drehwinkel ® der Drehmasse @ und dem Reaktionsmoment Jt rechts, 
eal das gleichzeitig Drillmoment der Welle ist. 

Bei geeigneter Vorzeichenfestsetzung gelten die Bewegungsgleichungen 


M =yl(o*— G), ag) 
die durch den Ansatz 
. P=geoosat, M —Mecoswt 
in 
M* =M+o0?69q, 
iibergehen. Thre Auflésung nach M und y lautet 
DM = M* —@? Oq’*, 
ee (2.4) 


Dp =—yI M* +¢%, J 
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wobei D = 1 — w/a und 1/ylO = w? gesetzt wurde. Dabei ist @, die Kigenfrequenz des Schwin- 
gers, dessen Eigenschwingung der Lésung (2.4) iiberlagert werden kénnte. 
Ist der Schwinger nach Abb. 6 angeordnet, so lauten die (2.3) und (2.4) entsprechenden Glei- 
chungen (bei Beachtung der oben getroffenen Vorzeichenwahl) 
M* =M—o’@0q, 
p* eka 
DM = M* +@? Oo", 
Dp=yiMetor. | 


Diese Vorbetrachtungen verwenden wir bei dem Zusammenbau 


i J a Abb. 6. Erzwungene Schwingungen 
unserer 2 n Elementarsysteme, wenn eine Eigenschwingungsform be- der Elementarsysteme. 
(4) 44 x 4; (+8, % Pt eat 


5 * * # 
Myr | by : MM ( 0), My Mae he Seer at Niet : Mya 7 
Pr-1 PR-1 Pr Pr 

Fr-7 Py, Y, Fra7 


Abb. 7, Zusammenbau der Elementarsysteme bei einer Eigenschwingungsform Py, die verschieden ist von der Kigenschwingungsform 
héchster Ordnung a : 


liebiger Ordnung dargestellt werden soll. Abb. 7 zeigt einen Ausschnitt aus der Reihe der Elementar- 
systeme. Von links nach rechts gelten die Gleichungen 


Mi_1 = M, + w? (1 — Bu—-1) Ox—1 G1» 


gta=yoh M+ o-, sa 
Mi = M,— w? 6, On Pe » (2 8) 
ger = —y (L—on) Mit pu | 
Mi = My 41 + w? (1 — Br) On Gx » 2.9) 
Oe = 7 Oerilesi Masi + Ge > oe 
Mini = My41 — w? Basi On41Pr-+1> | (2 10) 
Vi = —y (1 — 0441) 1 Magi + Pr41 - 


Die Elimination von Mj aus (2.8) und (2.9) und die von gj_; aus (2.7) und (2.8) fiihrt natiirlich 
auf die bekannten Schwingungsgleichungen der Welle 


Mi41— M, = — 07 O; ve ire, Do cecal) (2.11) 
Pr — Pr—-1 = yb, M;, (k = ik A EO! (2.12) 

mit My = Ma+1 =|) 
Die Gleichungssysteme (2.7) bis (2.10) werden nun umgekehrt. Nach (2.4) und (2.6) lauten ihre Auf- 


lésungen 


DM, = Mi_-1— ? (1 — Bri) x1 Pi» | 


en lis} 
D gra = — 7 Oe. Mea + Pe—1 > ( ) 
D M;,, = Mk + w? By Ox via » 
s é (2.14) 
Dox =y (lL— aon), Mi + pr—-1> 
D M,..1 = ME — 0? (1 — Bx) On Vi » | (2.15) 
Do, = —Y On+1be41 ME + fe » j 
D M,41 = Mi+it+ w* Bri Ox +1 Vie » | (2.16) 
D pr+i =y¥ (i — O41) Mi +1 a Pi . 
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Wegen (2.1) haben die Determinanten aller Gleichungssysteme (2.7) bis (2.10) denselben Wert 


2 


D=1— (2.17) 


a) 
wo? ? 
wobei @y die fiir alle Elementarsysteme gleiche Eigenfrequenz bezeichnet, die héchste der Welle. In 
den Gleichungssystemen (2.13) bis (2.16) mu demnach w = a, ausgeschlossen bleiben. Sie gelten 
fiir alle Eigenschwingungsformen gy, mit Ausnahme derjenigen héchster Ordnung. 


Wir eliminieren M;,,; aus (2.15) und (2.16), ferner q;, aus (2.14) und (2.15) und erhalten 
Mi... — Mk = — w? [(1 — Bx) Ox + Beo1 Oe+1] ve (leva, eee een (2.18) 
Pe — Vir = y [(1 — ax) L, + On+i lesa] Mi (k = i 2p © svostg sh) (2.19) 
mit My = M; = 0: 

Der Vergleich von (2.18) und (2.19) mit (2.11) und (2.12) liegt jetzt nahe. Beide Gleichungs- 
systeme stellen Eigenwertprobleme dar, deren Lésungen Eigenschwingungsformen unserer Welle 
sind. Sehen wir einmal davon ab, daB im zuletzt erhaltenen Gleichungsystem die Kigenschwingungs- 
formen durch Mj und gj charakterisiert sind, statt wie sonst tiblich durch M;, und gy, so geht 
aus der Herleitung des neuen Eigenwertproblems (2.18), (2.19) klar hervor, das die Eigenschwin- 
gungsform héchster Ordnung hier nicht mehr als Eigenlésung auftritt. Das alle Eigenschwingungs- 
formen der Welle umfassende Eigenwertproblem (2.11), (2.12) ist demnach mit Hilfe der bekannten 
Eigenschwingungsform héchster Ordnung iiber die Transformation der Mj, @;, in die Mk, yi um 
eine Stufe reduziert worden. 


In (2.18) und (2.19) fiihren wir nun die Abkiirzungen 
Of =(1— Br) Ox + Bry1O%41, Bo=0, Br=l, (2.20) 
le = (lL—on) he +Ons1he ys (2.21) 


ein und gewinnen dadurch die auf eine naheliegende dynamische Analogie hinweisende Gleichungs- 
form 
Mi,1— ME =—o@* Of pe (k =0,1,2,...,n—]), (2.22) 
Ve — Pia = yuMe. (kes lil en 35 hl) (2.23) 
muti — 2 Oe 
Die Mj deuten wir als Drillmomente einer Welle, die mit den Drehmassen@; besetzt ist. Die gi 
seien die Winkelamplituden dieser Drehmassen, die If die Langen der Wellenabschnitte. Die 
Definitionsgleichungen (2.20) und (2.21) der Of und If zusammen mit Abb. 3 und 4 lehren, wie 
verbliiffend einfach unsere neue Welle gebildet werden kann, wenn die Eigenschwingungsform 
héchster Ordnung der urspriinglichen Welle bekannt ist. Diese Konstruktion verdankt man 
Baranow. 

Die neue Welle besitzt eine Drehmasse und damit einen Freiheitsgrad weniger als die alte Welle. 
Beide Wellen stimmen in ihren Eigenfrequenzen iiberein, abgesehen von der héchsten Eigenfre- 
quenz der alten Welle, die der neuen fehlt. 

Es ist wohl unmittelbar einzusehen, da alle obigen Gleichungen ihre Giiltigkeit behalten, wenn 
man unter (py irgendeine beliebige Eigenfrequenz der Welle versteht, da® also die Reduktion auch 
mit einer bekannten Eigenschwingungsform beliebiger Ordnung vorgenommen werden kann. Man 
iiberzeugt sich aber sofort, daB8 dann alle anschaulichen dynamischen Betrachtungen, die unsere 
Rechnungen bis zum SchluB begleitet haben, hinfallig werden. 


3. Die Transformationen der Eigenschwingungsformen. Unter W wollen wir in Zukunft die alte 
Welle verstehen, unter W* die aus W durch eine Baranowsche Reduktion entstandene. Fiir die Be- 
stimmung der Eigenschwingungsformen von W mit Hilfe des Baranowschen Reduktionsverfahrens 
ist es wichtig zu wissen, wie sich eine Higenschwingungsform von W* in diejenige gleicher Frequenz 
von W transformiert. Diese Transformationen sind sehr einfach den Gleichungssystemen (2.13) 
bis (2.16) zu entnehmen. 


Aus (2.13) und (2.14) ergibt sich 


DG. — Yr) = y by, [o, Ma + (1 — o,) Mi] (3.1) 
und nach (2.12) 


DM; = 0, Mi1+(1—o,) ME (k =1, Paice Np (3.2) 
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Dabei ist M§ = M% = 0. Aus (2.14) und (2.15) folgt 


D(Mi.1— Mi) = — w? O; [Bu pia + (1 — Bx) gf] » (3.3) 
also nach (2.11) 
D = Bx Pk-1 a (1 — Bx) Vi (k = 0, ll. 2 Sceted) n) . (3.4) 
Wegen f, = 0 undf,, = 1 wird 
Dp = > (3.5) 
Do, = O.51 (3.6). 


Die Gleichungen (3.2) und (3.4) zeigen, wie eine Eigenlésung des reduzierten EKigenwertproblems 
(2.18), (2.19) umgerechnet werden kann in die zum gleichen Eigenwert gehérige Lisung des urspriing- 
lichen Eigenwertproblems (2.11), (2.12). 

Mehr als die Drillmomente M, sind gewéhnlich die Winkelamplituden y, der Eigenschwingun- 
gen von Interesse. Ihre Riicktransformation (3.4) von W* nach W laBt sich sehr einfach zeichne- 
risch bewerkstelligen. Um dies einzusehen, schreiben wir (3.4) zunachst als Proportion 

Vi1—Do, _ DH 
Spee © Br ee 


und schlieSen die Sonderfalle k = 0 und k = n, die bereits durch (3.5) und (3.6) erledigt sind, aus. 
Nach (2.1) ist 

Bi(1 — on) b, = (1 — Br) ona ears (3.6) 

so da (3.5) als Proportion von Strecken gedeutet werden 
kann: 

G31 — Dy, __Dey—?; 
— : Sod 
(1 — o,) h, Orr diy om) 


Wir erinnern uns, da gf urspriinglich definiert ist als 
Verdrehwinkel der Welle W an derjenigen Stelle, an der 
ihre Eigenschwingung héchster Ordnung den Knoten 
N,,41 besitzt. Tragt man also die als Strecken dar- 
gestellten Winkelamplituden gyj_; und gj der Welle W* 
in den Punkten Ni: bzw. N41 der Welle W als Ordi- Abb. 8. Zeichnerische Darstellung der Proportion (3.7). 
naten auf, so gibt (3.7) die einfache Konstruktion von 

Dg, nach Abb. 8. Wegen D < 1 ist gy, immer gréfer als Dp. Da aber alley, einer Figen- 
schwingungsform nur bis auf einen konstanten Faktor bestimmt sind, darf man bereits die D 
als die gesuchten Amplituden der Drehmassen ©, von W ansehen. 


Oo, One 7 


Oe Lies (1-06.44) Vag Ore Leyg 


Abb. 9. Die Riicktransformation einer Eigenschwingung des reduzierten Systems W* in diejenige 
gleicher Eigenfrequenz des Systems W. 


Meistens wird es bequemer sein, die Strecken D q,, bereits im System W* abzugreifen und sie 
dann nach W zu iibertragen, wie es Abb. 9 zeigt. 
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Die Dp, sind gleich den Verdrehwinkeln der Welle W* in den Punkten T,, (nicht zu verwechseln 
mit den Knoten der Eigenschwingung héchster Ordnung von W*), also dort, wo bei W die Dreb- 
massen @, sitzen. Wortlich gilt dies auch fiir die 
Wellenenden, die wir nach (3.5) und (3.6) gesondert zu 
betrachten haben. Man stelle sich vor, dab W* tber- 
stehende Wellenenden besitzt, die sich bei einer elasti- 
schen Schwingung der Welle wie starre Kérper drehen 
(Abb. 10). Ihre Drehwinkel sind deshalb gleich den Dreh- 
winkeln der Endmassen 05 bzw. O*_, . 

Im Baranowschen Verfahren wird, nachdem die 
Eigenschwingungsform héchster Ordnung von W* er- 
mittelt worden ist, weiter reduziert auf ein System W**, 
das ‘zwei Freiheitsgrade weniger als W besitzt, und so 
fort, bis nach insgesamt n— 1] Reduktionen ein System 
————w*von zwei Drehmassen entstanden ist, dessen Eigen- 

frequenz gleich der niedrigsten von W ist. Bei jedem 
Reduktionsschritt benutzt Baranow die Kutzbachsche 
' Seileckskonstruktion, die, wie an vielen Beispielen fest- 
Abb. 10. Die Riicktransformation der Winkelamplitude ° : : re 

ene eI gestellt wurde, die Eigenschwingungsform :héchster 
Ordnung und ihre Frequenz mit zufriedenstellender 
Genauigkeit liefert. Auf eine Erérterung dieser Konstruktion darf hier verzichtet werden, da man 

ausfiihrliche Beschreibungen, auch des Baranowschen Verfahrens, bei Klotter und Haug findet. 

Da wir jetzt wissen, wie die Eigenschwingungsformen héchster Ordnung der einzelnen reduzier- 
ten Systeme auf W zuriicktransformiert werden, kénnen wir uns durch simple zeichnerische Kon- 
struktionen samtliche Eigenschwingungsformen von W verschaffen, unmittelbar im AnschluB an 
das Baranowsche Verfahren. Baranow wird diese Transformationen nicht gekannt haben; denn 
sonst hatte er sich nicht mit der Bestimmung der Eigenschwingungszahlen allein begniigt. Ob er 
tiberhaupt einen allgemeinen Beweis fiir die Exaktheit seiner Reduktion gekannt hat, ist nicht fest- 
zustellen. 


a G 


asd 


Abb. 11. Reduktion bei festgehaltenem linken Wellenende. 


Erwahnt sei noch, daB man durch Umkehrung des Baranowschen Verfahrens, von einem Zwei- 


se ie as ausgehend, eine Schwingerkette mit vorgeschriebenen Eigenfrequenzen aufbauen 
cann, 
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4. Festgehaltenes Ende der Welle als Randbedingung. Auch bei dieser Randbedingung laBt sich 
eine Baranowsche Reduktion durchfiihren. Wir erlautern dies an Abb.11. Bild a) zeigt W mit ein- 
gespanntem linken Wellenende gy = 0, Bild b) die Elementarsysteme gleicher Eigenfrequenz wy 
und Bild c) das reduzierte System W*. Jetzt ist x, Null, und die Aufteilungsfaktoren f, beginnen 
mit fj, f,, ..... Beachtet man dies, so hehalten simtliche Transformationsgleichungen fiir Winkel, 
Torsionsmomente, Langen und Tragheitsmomente Giiltigkeit. 

Auch die Baranowsche Seileckskonstruktion la%t sich sofort auf diese Randbedingung iiber- 
tragen. 


5. Zusammenfassung.. Aus anschaulichen dynamischen Betrachtungen heraus wird ein von 
Baranow ohne Beweis angegebenes Verfahren zu Ermittlung der Torsionseigenfrequenzen einer 
Welle entwickelt und seine Exaktheit nachgewiesen. Ferner wird gezeigt, wie man dem Bara- 
nowschen Verfahren mit geringer zusatzlicher Zeichenarbeit sémtliche Eigenschwingungsformen 
entnehmen kann. 


(Eingegangen am 23. Oktober 1954.) 
Anschrift des Verfassers: Prof. Dr.-Ing. H. Schaefer, Braunschweig, Technische Hochschule. 
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Die Abbildung eines allgemeinen Schwingungssystems 
auf eine einfache Schwingerkette* 


Von S. Falk 


1. Einfiihrung, Ziel der Arbeit. Es handelt sich im folgenden um ungedimpfte mechanische 
Schwingungssysteme mit n Freiheitsgraden, fiir welche potentielle Energie U und kinetische Energie ic 
quadratische Formen in den gewahlten Koordinaten x; und deren zeitlichen Ableitungen sind: 


2 Ges S x6, % =41 0, ee = tibet = EOE (1) 
Resi jel p= 
mit 
Ul = (Xy %peee Xn)o El = (Hy, My--- Xn), W = (Gu), B= (bux)- (2) 


Dabei werden die elastischen Eigenschaften des Schwingungssystems durch die ,,Federungs- 
matrix“ 9{ und die Tragheitseigenschaften durch die ,,Trigheitsmatrix“ § ausgedriickt. Die Koor- 
dinaten fassen wir zur einzeiligen Matrix y’, ihre Ableitungen zu 7’ zusammen. Quadratische Ma- 
trizen werden durch groBe, einspaltige Matrizen (Spaltenvektoren) durch kleine gotische Buch- 
staben, einzeilige Matrizen (Zeilenvektoren) durch kleine gotische Buchstaben mit einem Strich 
gekennzeichnet. Ganz allgemein bedeutet der Strich den Ubergang zur transponierten Matrix. 


Die Energieformen (1) lassen sich stets so schreiben, daB Qf und $ symmetrisch, also 2 = QI’ 
(aj, = a,;) und B =’ (bj, = b,:) werden. Beide Matrizen seien zunachst reell. Q ist seiner 
physikalischen Natur nach positiv definit, d. h. die quadratische Form 2 T ist fiir reelle Werte x; 
stets positiv und nur gleich Null, wenn samtliche x; verschwinden, das System also in Ruhe ist. 
Y setzen wir fiirs erste positiv definit oder positiv halbdefinit voraus, d. h. die potentielle Energie U 
sei fiir beliebige Werte x; niemals negativ; zusammengefaBt: 


Y% =" positiv (halb)definit, % = %’ positiv definit . (3) 


Unter dieser Voraussetzung sind die Quadrate w; der Eigenkreisfrequenzen positiv oder gleich 
Null, die Schwingungen des Systems also stabil. 


Wir betrachten nun vier verschiedene Klassen von Schwingungsgebilden, die alle den Be- 
dingungen (3) geniigen und in Tabelle 1 zusammengestellt sind. 


Tabelle 1. Einteilung mechanischer ungedémpfter Schwingungssysteme. 


Klass Federungsmatrix Trag heitsmatrix Mehrfache 
asse positiv (halb) definit positiv definit Eigenfrequenzen sind 
I: (Abb. 1) Vollmatrix Vollmatrix $ moglich 
II: Versetzte Welle Dreiermatrix ©; Zeilen- und Diagonalmatrix D nicht moglich 
(Abb. 2) Spaltensummen sind im 
allgemeinen nicht Null 
III: Durchlaufende Welle Dreiermatrix ©; Zeilen- und Diagonalmatrix $ nicht moéglich 
(Abb. 3) Spaltensummen sind Null 
(gegebenenfalls bis auf erste und 
letzte) 
IV: Reihe von Elementar- Diagonalmatrix } Diagonalmatrix & dglich 
schwingern (Abb. 4) : ae 


* Gekiirzte Fassung der Dr.-Ing.-Dissertation des Verfassers: ,,Die Berechnung der Torsionseigenschwin- 


gungszahlen verzweigter Maschinenanlagen“, Berichter: Prof. Dr.-Ing. H. Sch Mitberi : -Ing 
= Deen Braue eos! > g chaefer, Mitberichter: Prof. Dr.-Ing. 
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Das Schema der Abb. 1 ist nur symbolisch zu verstehen; es vertritt alle denkbaren unge- 


dampften Schwingungssysteme, deren Matrizen der Voraussetzung (3) geniigen. 2{ und B sind im 
allgemeinen Vollmatrizen, d.h. alle ihre Elemente 


sind von Null verschieden. Die Systeme der Klas- 
sen IT, III und IV sind Sonderfalle von I; als Repra- 
sentanten wahlen wir masselose elastische Wellen ° 
mit konstanten Drehsteifigkeiten Cy, (cm kg), die 
mit diinnen Scheiben der Tragheitsmomente J; ; 
(cm kg sec~?) besetzt sind. Die dimensionslosen “P?  A¥semeines unged Teuhoimadoe ne ak 
Winkelauslenkungen aus der entspannten Lage 


seien die Koordinaten x;. Jetzt unterscheiden 4 | ———-— 
wir: Ea en 
4 
Klasse II: Versetzte Welle mit Zahnradern. Die | On 
Federungsmatrix heiBe , die Tragheitsmatrix D. Z 
Man zeigt leicht, daf& bei der getroffenen Koordi- ee SEPM OE EES 


natenwahl D Diagonal- und ( Dreiermatrix ist, 
d. h. D enthalt von Null verschiedene Elemente nur 
in der Hauptdiagonale (d;, = 0 fiir i ~ k) und @ nur 
in der Hauptdiagonale und den beiden benachbarten 
Nebendiagonalen. 


Klasse III: Durchlaufende Welle mit Federungs- 
matrix © und Tragheitsmatrix .§: Auch © ist 
Dreiermatrix, aber spezieller als ©, und 3 Diagonal- 
matrix: 


J . ae is ; In, 


eee —C,, . 0) soodiseoc 0 i 0 0 
=! Cis ‘ce + Cy5 =< (One Ao Fotrt a0 0 0 Js 0 ‘a 
c= | 0 — Cy, Cog + Coq 2 EP ROR A 0 F 5 =| 0 0 Oe (4) 
\ 0 0 ca Cyd, n GaSe nm te al Mb 0 0 a 
Die zugehérigen Energieformen nach (1) sind 
n—1l1 n 
2U= C, x} Se Pas Ci, i411 (4%: — % 4 1)? + Cra, oee = Spx. (5) 


aw Sil 


In © sind samtliche Zeilen- und Spaltensummen gleich Null bis auf die erste und letzte. Ist die 
Welle am linken (rechten) Ende frei, so verschwindet mit C, (C,,) auch noch die erste (letzte) Zeilen- 
und Spaltensumme. Im Fall C, = C,, = 0 tritt genau einmal die Eigenfrequenz Null auf: die freie 
Welle kann unverdrillt als starrer Koérper rotieren. 

Diese Klasse von Schwingungsgebilden hat Grammel? eingehend untersucht und eine Anzahl 
von theoretisch ebenso interessanten wie praktisch wertvollen Satzen bewiesen, inshesondere 
Schranken fiir die Eigenfrequenzen angegeben. 

Klasse IV: Kette von Elementarschwingern. Hier sind Federungsmatrix % und Tragheits- 
matrix & beide von Diagonalgestalt; das Schwingungssystem ist vollstandig entkoppelt, und die 
Eigenfrequenzen lassen sich ohne weiteres ablesen. 

Die Ermittlung der Eigenschwingungszahlen und -formen von Schwingungssystemen der 
Klasse I ist im allgemeinen mit erheblichem Rechenaufwand verkniipft. Die Systeme der Klasse IIT 
dagegen bereiten weniger Schwierigkeiten, und es sind wegen ihres haufigen Auftretens in der 
technischen Praxis (z. B. als Kurbelwellen) in der Ingenieurmechanik zahlreiche praktische Ver- 
fahren entwickelt worden, die fast alle, ob ausgesprochen oder nicht, auf der speziellen Gestalt der 
Matrizen YJ und & — oft auch hier noch C, = C, = 0 voraussetzend — beruhen und deshalb auch 
allein auf diese Klasse beschrankt bleiben. Eine besonders einfache zeichnerische Methode ist die 
von Baranow!, die auch fiir beiderseits eingespannte Wellen gilt. Haug* wies im Jahre 1952 erst- 

1G. Baranow, Z. VDI 76 (1932), S. 184 und Metal Ind. Herald Moscow 11 (1931) S. 60. (In deutschen 
Bibliotheken nicht vorhanden.) 


2 R. Grammel, Ing.-Arch. 14 (1943/44), S, 213. 
’ K. Haug, Die Drehschwingungen in Kolbenmaschinen, Berlin-Gottingen-Heidelberg 1952. 


Ingenieur-Archiy 


316 Falk: Die Abbildung eines allgemeinen Schwingungssystems usw. 


malig darauf hin, daf die noch in dem groBen zusammenfassenden Bericht von Klotier! als Néahe- 
rungsverfahren bezeichnete Seileckskonstruktion von Baranow e x a k t richtig ist, und H. Schaefer 
erbrachte nicht nur einen neuen Beweis fiir diese Tatsache, sondern zeigte iiberdies, wie man der 
Baranowschen Konstruktion auBer den Eigenschwingungszahlen auch noch samtliche Eigen- 
schwingungsformen mit wenigen Strichen entnimmt. Es lag daher der Versuch nahe, dieses so 
einfache graphische Verfahren auch fiir die umfangreichere Klasse I von Schwingungssystemen 
nutzbar zu machen, was in der Tat gelingt. 

Dazu bemerken wir als erstes, daB der in der Tabelle 1 angegebenen Einteilung keinerlei mecha- 
nische, sondern lediglich mathematische Bedeutung zukommt: jedes vorgelegte mechanische 
Schwingungsgebilde, das den Voraussetzungen (3) geniigt, 1aBt sich nach Belieben je nach Wahl der 
(nicht immer anschaulichen) Koordinaten als System einer der vier Klassen darstellen. Oder anders 
ausgedriickt: hat man zunachst irgendwelche Koordinaten x,,.. + %n gewahlt, in denen das System 
als zu einer der vier Klassen gehérig erscheint, so laBt es sich nachtraglich durch eine lineare Koor- 
dinatentransformation immer in ein solches der anderen drei Klassen uberfiihren, wobei samtliche 
Eigenfrequenzen erhalten bleiben, Die Transformation der Klassen I, IJ und III auf IV, also auf 
eine Kette von Elementarschwingern, ist als sogenannte Hauptachsentransformation hinlanglich 
bekannt; sie erfordert jedoch die Kenntnis der — in der Praxis doch gerade gesuchten — Kigen- 
frequenzen, weshalb ihr im allgemeinen nur theoretischer Wert zukommt. Wir stecken deshalb 
unser Ziel weniger weit und begniigen uns mit der Transformation der Klassen I und II auf III, 
zeigen also, wie man mit Hilfe einer linearen Koordinatentransformation ein vorgegebenes mecha- 
nisches Schwingungsgebilde, das den Bedingungen (3) geniigt, auf eine durchlaufende Welle zu- 
riickfiihren oder ,,abbilden“‘ und somit einerseits dem Verfahren von Baranow-Schaefer und zum 
anderen den EinschlieBungssatzen von Grammel zuganglich machen kann. 


J em 
giv em 
eS 2G G 6c 
Y yp ¥3 
_—_—_—SP 
m[kg cmu'see’], c|kg em"). a [cm]; 
2 
J [kg cm see*] ; C[kg em]; yi[em); 4 a 
* 7 fe 
yi lt] 4 Ys 
Abb. 5. Schwingungssystem der Klasse I Abb. 6. Bildwelle (bzw. Bildkette) des Originalsystems 
mit drei Freiheitsgraden. von Abb. 5. 


Vorweg ein einfaches Beispiel: fiir das System der Abb. 5 ist offenbar 


2 U =c (x, — 4%)? + ¢ (x, — xg) + ¢ (x, — x)? +622, 2T= : mx? + ma? + mx. (6) 
Wie man leicht nachrechnet, werden mit den durch 
X= Vio % = Yet Ys. %3 =Ye— Xs (7) 
eingefiihrten neuen Koordinaten y, potentielle und kinetische Energie 
2 U = c(¥1— ¥2— Ys)? + €(¥2 +3 — Yo + ¥s)® + €(¥1 — Yo + ¥s)® + €(¥,— 9s)? 
= 2C( Vin) a a Vs) ives (8) 
Sb 3 ee 
2T= wees + m(¥2+ 3)? + m(¥,— 5)? Sk, +2 yo -- 2m ye. 


Das in (6) stérende Glied ¢ (x, — x3)? ist in (8) verschwunden; die Transformation (7) hat somit 
ihren Zweck erfillt: ein Vergleich mit (4) zeigt, daB die Bildwelle des Originalsystems das Aus- 
sehen von Abb. 6 hat, auch der entsprechende Translationsschwinger (Schwingungskette) ist an- 
gegeben. 

Es entsteht nun die Frage, ob es stets eine solche Transformation gibt, und wenn ja, wie man sie 
praktisch findet. Wir zeigen im Abschnitt 2, daB die Transformation unter den Voranccet ieee (3) 
immer, und zwar nicht nur auf eine Weise, moglich ist. Der Rechenaufwand liegt dabei in der 
gleichen GréBenordnung wie der anderer Verfahren, welche die Frequenzgleichung det (Q(—k %) =0 


1 K. Klotter, Ing.-Arch. 17 (1949), S. 1. 
 H. Schaefer, Abh. Braunschweig. Wiss. Gesellschaft 5 (1953), S. 141. 
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explizit aufstellen, um die Eigenfrequenzen daraus zu berechnen, eine Aufgabe, die das Verfahren 
von Baranow-Schaefer ja gerade zeichnerisch lost. 

Fir die praktische Durchfiihrung der Transformation ist es bequem, mit Hilfe zweier zweck- 
maBig gewahlter Vergleichsgréfen C (cm kg, kg em-! usw.) und J (cm kg sec~®, kg cm™ sec”? usw.) 
alle Elemente der Federungs- und Tragheitsmatrix dimensionslos zu machen und ebenfalls dimen- 
sionslose Koordinaten zu wahlen. Dann sind k; = w? C/J die ebenfalls dimensionslosen »,Kigen- 
werte“ des Systems; die Ausschlige Xity +++ Xin» die zu einem Eigenwert k; gehéren, also die zur 
Eigenfrequenz w; gehérige Eigenschwingungsform, hei®t dann der »Higenvektor“ 7;. 

2. Die Abbildung eines allgemeinen Systems auf eine durchlaufende Welle. Wir fiihren jetzt die 
alten Koordinaten x; als Linearkombinationen neuer Koordinaten yi ein, die wir zum Vektor {) 
zusammenfassen: 


Xi = ey Vik Ve (= 1 2 sot) Coder Li UY ek == (Gees tne (9) 


Dann werden potentielle und kinetische Energie von der Form 
20-—~-AXc=yRURy, 2T=—rBrHpy NK BRg. (10) 


Gesucht sind solehe Koordinaten yi, in denen sich U und T deuten lassen als potentielle und kine- 
tische Energie einer durchlaufenden Welle. Es ist indessen bequemer, das Schwingungssystem der 
Klasse I mit den Matrizen Q{ und 8 zunachst auf die versetzte Welle der Klasse II mit den Ma- 
trizen (§ und D abzubilden; der Schritt von II nach III ist dann leicht nachzuholen. Wir verlangen 
also vorerst nur, da 


RUR=G und KRVSR=D (11) 


wird. Die neuen Matrizen & = (g;,) und D = (dj) erfiillen ebenfalls die Bedingungen (3), ihre 
Elemente sind die Bilinearformen 


Lik = r; 2 ie, wil Gh, = r,B Vk: (12) 


Unsere Forderungen an die gesuchte Transformationsmatrix Ji = (r,,...Y,) sind also nach (11): 
a) G soll eine Dreiermatrix sein, d. h. es miissen alle Elemente g;,, fiir welche die Differenz der 
Zeiger t und k dem Betrage nach gréfer als eins ist (z. B. g., g,, usw.) verschwinden: 


Sn SF N= 0 fir, lt— kl > 1. (13) 


b) D soll eine Diagonalmatrix sein, d.h. die gesuchten Spalten (Zeilen)-vektoren 1; (r;) von 
Ji (MH) miissen verallgemeinert orthogonal zur gegebenen Matrix § werden: 


dz, =t/Br, =—0 fir iXxk. (14) 


Eine hier ibergangene Uberlegung zeigt nun, da sich die Vektoren ry}, r, ... nach willkiirlicher 
Wahl von rj aus einer einfachen Rekursionsformel bestimmen und mit ihnen die Matrizen & und ® 
schrittweise aufbauen lassen: sind die ersten m Vektoren 1, ermittelt, so liegen damit auch die 
m-reihigen Hauptabschnittsmatrizen oben links in ( und ®D bereits fest. Es sind Verfahren be- 
kannt, die zwar auch die Transformation von 2{ und % auf % und ® leisten!, die aber keine Rekur- 
sionsformel nach Art von (15) entwickeln und auch die speziellen Formen © und § nicht zum Ziele 
haben, somit auch nicht auf den mechanischen Fragenkreis des Abschnittes II] anwendbar sind. 
Die grundlegende Rekursionsformel zur Berechnung der Vektoren y; lautet nun: 


Cea , esi 8mm of §m—1, m i 
w= —r, UB 4 ‘Um +g ie (15) 


Sm+1 mm m—1l1,m—1 


fin leon iit. 2, — 0 0nd «do, 1h sp ee Oe 


Zum Beweise nehmen wir an, die Vektoren 1; bis ry, seien bereits bekannt, die Bedingungen (13) 
und (14) also fiir i,k = 1, 2,...m erfiillt, und multiplizieren (15) von rechts mit dem Spalten- 

1 Eine solche Transformation findet sich bereits bei C. G. J. Jacobi (Jacobi, Werke Bd. 6, S. 318—321), wie 
mir freundlicherweise Herr Prof. Bilharz, Wiirzburg, mitteilte. In neuerer Zeit ist das Problem auf andere 
Weise angefaBt worden von C. Lanczos, J. Research, Nat. Bu. Stand. 45 (1950) S. 255 und W. E. Arnoldi, 
Quaterly of Appl. Math. 9 (1951) S. 17. 
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vektor§ rp. Dann stellen wir fest, da die linken Seiten von (15) verschwinden: 


/ 9 
Py ee ee ee eae. 
Sm+1 
ene Come SEL aeeey ae aa) aes. onal mal = 
P = . Sein §m, m—l om — dees ae 
seem B tm dam 
Pp=m: Si =— 8mm + &mm* | +0=0; 
mt+1 mm 
zusammengefaBt : 
te te ee furas pil anager atie (16) 


Das heiBt: der nach (15) ermittelte neue Vektor Ym +1 Steht in der Tat verallgemeinert senkrecht 
auf Y;....Ym, Wie es (14) verlangt. 

Jetzt denken wir uns die Formel (15) fiir irgendeine Nummer q < m — | hingeschrieben und 
von rechts mit Yr,, ,; multipliziert, dann wird wegen (16) 


0=—r1, Wtm41 +9040, also g,n41 = 0 fir q=1, 2,...m—l1. (17) 


Die neu zu berechnenden Elemente der m + 1ten Zeile und Spalte von & sind somit nach (17) 
gleich Null auBer gn, m1 = Zm+1,m UNA Bm 41, m41; damit ist auch (13) erfillt. 


Zum Schlu8 multiplizieren wir noch (15) mit § r,,, 1 und bekommen 


d 
m+i1,m+1 Fi 
erga a oder Cite mee = Sen eet fiir Tay is Ai (18) 


Sm +1 
Der von Null verschiedene Faktor s,, , ; ist dabei noch ganz beliebig; man kénnte ihn nach jedem 
Schritt so bestimmen, da @® in die Einheitsmatrix & tibergeht, also alle d;; = 1 werden, oder 
auch so, daB Zeilen- und Spaltensumme der Matrix ( verschwinden (gegebenenfalls bis auf die 
letzte), so da dann (} von vornherein die spezielle Form © (4) hat. Am einfachsten ist es jedoch, 
zunachst s,, ;; = 1 zu wahlen, dann wird nach (18) 


(Pats eae te 1s Ps fiir ae dee (19) 


Das heifit: die Hauptdiagonalelemente d;; der Matrix Q brauchen tiberhaupt nicht berechnet zu 
werden, sie sind auBer dem ersten d,, gleich den negativen Nebendiagonalelementen von (; diese 
selbst mitissen demnach negativ sein, denn die d,; sind als quadratische Formen der positiv definiten 
Matrix § stets positiv. 

Wir schreiben nun, (19) benutzend, die Rekursionsformeln (15) noch einmal ausfiihrlich hin; es 
ist 


tr, =v, frei wahlbar, damit ist d,, =1,S1, zu berechnen, 
r= —r1 ABS ae dag = — £15 | 
atk 
r, =—1, UB + a ue a i> dss =— gy, | (15a) 
t= — 1, US4 + By + ey, day = — 843 | 
trckaeer ac Cen eee j 


und so weiter. Damit sind ® und D ermittelt. Es ware natiirlich iiberfliissig, die zugehorige ver- 
setzte Welle wirklich aufzusuchen, da wir sie ja mit Hilfe der Faktoren s; sogleich wieder auf 
eine durchlaufende Welle abbilden wollen. Zu diesem Zwecke fiihren wir noch einmal neue 
eran 2; = Y;/s; baw. y; = s; z; ein, oder, wenn wir die s; zur Diagonalmatrix G zusammen- 
assen, 


D625 SRY i C3. s "Gee (20) 
Damit werden potentielle und kinetische Energie in den neuen Koordinaten z; 
2U=yAr=h Gy=7 SG GS4=2' C3, (21) 


2T=ESt=) DH= 74 SDS5 =F 33. (22) 
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und die neuen Matrizen © und }, die wiederum den Bedingungen (3) geniigen, sind 
C=C’ GG, Cik = Chi = Si Sik Sk » (23) 
v=ODS, tig = di; 8}. (24) 
[ Die Elemente Cik in (23) sind nicht mit den Drehsteifigkeiten C;, in (4) zu verwechseln.] Jetzt 
verlangen wir, daB die Zeilensummen — und damit wegen der Symmetrie auch die Spaltensummen — 
von © gleich Null werden bis auf die letzte. Die Matrizen © und § entstehen nach (23) und (24) 
aus ( und , indem deren Zeilen und Spalten der Reihe nach mit s,,... s, multipliziert werden. 


Ziehen wir in & aus jeder Zeile den gemeinsamen Faktor s; wieder heraus, so verbleibt zur Berech- 
nung der GréBen s; das homogene Gleichungssystem 


$11 $1 1 S12 Se = Us 
S21 $1 1 S22 So + 8o3 $3 =0, 
§32 Se + 833 $3 T S34 $4 = 0, (25) 


0) 16 01 0/0) ehja\"e')e) “a! 8! (6 
SON ee Le 18 ei fare, 68210) "e: Silene sei el 0) a 6) (G auley 6. 0. 9. 0), o/h @:@. @. iq 8) 0) Sle) 0: 76 


—-—-—— > 


Sa en 2 Sn —2 ew ikl Sips aR yelp Si = OF 


durch das nach willkiirlicher Annahme von s, auch die iibrigen s; festgelegt sind. 


Ist nun eine versetzte Welle oder ganz allgemein ein System der Klasse II mit den Matrizen & 
und Q von vornherein gegeben, so kommt man ohne Benutzung der Rekursionsformeln allein mit 
den Faktoren s; nach (25) zum Ziel. Fiir den allereinfachsten und praktisch wohl auch wichtigsten 
Fall der nur einmal versetzten Welle findet man in der technischen Literatur sogenannte ,,Reduk- 
pronstormeln’ ; dort ist s,—---—=s, = 1, 5,41) =--- = ssl. 

Nun ist aber noch keineswegs gesagt, daB die nach (25) und (23) berechnete Matrix € sich wirk- 
lich als Federungsmatrix einer durchlaufenden Welle deuten 148t; wir haben vielmehr, wie ein 
Blick auf (4) zeigt, noch zweierlei zu verlangen: 

b,) die Nebendiagonalelemente von © miissen samtlich negativ sein, 

b,) die letzte Zeilensumme o, darf nicht negativ sein; denn in (4) isto, = — C,1,n + Cr—1n 
+ C, = C,, und das ist entweder positiv oder bei der rechts freien Welle gleich Null. 

Die Bedingungen b,) und b,) sind jedoch unter den Voraussetzungen (3) von selbst erfiillt, was 
wir jetzt zeigen wollen. Die Unbekannten s; verhalten sich wie die mit dem Faktor (— 1)' +! ver- 
sehenen n — | reihigen Unterdeterminanten U;, die nach Streichen der i-ten Spalte des Gleichungs- 
systems (25) ubrigbleiben: 


Sti Sikes bey sO Br Gras Ul (26) 
Diese GréBen U; lassen sich jetzt durch die Hauptabschnittsdeterminanten H; und die Neben- 


diagonalelemente von & ausdriicken als 


n—l n—l 
U; => ET meey (— ye 7. Sn gdot, == H;_; V; mit Ay = lie fl Par si ee (i = MA. eee n) (27) 
n—l1 
Nach (19) sind sdmtliche g, ,,; negativ, also miissen die Produkte V; = (—1)'+1 JT g,,,41 von 
=i 


a 
einerlei Vorzeichen sein. Falls nun 2{ und damit ( positiv definit ist, sind samtliche Hauptab- 
schnittsdeterminanten H; von ( positiv; ist aber Q{ und damit ( positiv halbdefinit, so ist jeden- 
falls H,, = det & gleich Null, aber nicht etwa schon eine der Hauptabschnittsdeterminanten H,,... 
H,,—1; denn angenommen, die erste verschwindende von ihnen ware H,, dann miBte wegen 


ss oral, Ocal H, — Zio41 Tis = 0 Boost Hy —1 (e < n) (28) 
H, , ; negativ sein, weil H, 1 > 0 angenommen und g,, 941 nach (19) ¥ 0 ist. H, +1 <0 aber steht 


im Widerspruch dazu, daB 2{ und & positiv halbdefinit sind, also gilt 
(Hn > 0 falls Q{ positiv definit , 
; Wee = 0 falls { positiv halbdefinit “| 


und da in (27) nur die nach (29) stets positiven Hauptabschnittsdeterminanten H,,...H,—1 auf- 
treten, folgt nun, daB die U; ebenso wie die V; und damit nach (26) auch die s; von Null verschieden 
und von einerlei Vorzeichen sind, was wiederum bedeutet, daB die Nebendiagonalelemente c;, von 
€ ebenso wie die g;, von & nach (23) negativ bleiben, womit die Bedingung b,) erfillt ist. 


Tha ee His ee 0 (29) 
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Nun hangen die Hauptabschnittsdeterminanten H* von © und H; von (, wie leicht einzu- 
sehen, folgendermafen zusammen: 


T= i, hes (30) 


Desa 
und da gezeigt ist, daB alle s; A 0 sind, gilt nach (29) auch 


Te oe fin > 0 falls Y{ positiv definit , | (31) 
[a =0 falls 2{ positiv halbdefinit . J 

Jetzt denken wir uns in der Determinante von © die ersten n — 1 Spalten zur letzten addiert: die 

Elemente der neuen letzten Spalte sind dann gerade die Zeilensummen der Matrix ©, die ersten 

n — 1 von ihnen verschwinden daher — gerade so haben wir ja die s; aus (25) bestimmt — und 

das letzte Element rechts unten ist gleich g,. Die Determinante nach der letzten Spalte entwickelt 

ergibt nun 


det © = Nahe =On et ’ On SSg = ’ (32) 


und es folgt ohne weiteres, da nach (31) Hy; _, auf jeden Fall positiv ist: 
On > 0 falls Y{ positiv definit; On = 9 falls Yf positiv halbdefinit . (33) 


Also ist auch die Bedingung b,) erfiillt. 

Damit ist gezeigt: Unter den Voraussetzungen (3) lassen sich mit Hilfe der Rekursionsformein 
(15 a) und der Gleichungen (25) die gegebenen Matrizen Q{ und ¥ stets tiberfiihren in zwei andere 
Matrizen © und & mit folgenden Eigenschaften: 

1) ist Diagonalmatrix mit lauter positiven Hauptdiagonalelementen, 

2) © ist symmetrische Dreiermatrix mit positiven Hauptdiagonal- und negativen Nebendiagonal- 
elementen; die ersten n — 1 Zeilen- und Spaltensummen von © sind gleich Null, die letzte Zeilen- 
und Spaltensumme ist nicht negativ. 

Die Matrizen © und ¥ lassen sich somit immer deuten als Federungs- und Tragheitsmatrix einer 
durchlaufenden Welle, die links frei und rechts eingespannt oder frei ist, je nach dem ob Y positiv 
definit oder positiv halbdefinit war. Mit anderen Worten: Besitzt das Originalsystem den Eigen- 
wert Null, so ist die Bildwelle beiderseits frei; wenn nicht, so ist sie rechts eingespannt. 

Sind nun Eigenschwingungszahlen w? und -formen 3; nach Baranow-Schaefer ermittelt, so er- 
halt man die Eigenschwingungsformen 7; des Originalsystems durch die einfache Riicktransforma- 
tion (20) 

jee mC, mit © (sets tae esate (34) 
Der Ausgangsvektor r, in (15a) und der erste Koeffizient s, in (25) sind frei wahlbar; der Vektor s, 1; 
enthalt genau n willkiirliche Konstanten, also gibt es auch oo” voneinander verschiedene Bild- 
wellen zum gleichen Original. Man kénnte versuchen, eine oder einige der Komponenten des Aus- 
gangsvektors als unbestimmte Parameter mit in die Rechnung aufzunehmen, um sie an gegebener 
Stelle nach praktischen Gesichtspunkten zu verwerten, etwa so, da einige Tragheitsmomente oder 
Drehsteifigkeiten der Bildwelle vorgeschriebene Werte annehmen, es zeigt sich aber, da selbst 
ein einziger mitgefiihrter Parameter schon nach dem ersten Schritt die Rechnung derart anschwellen 
14Bt, daB fiir die Praxis sich ein solches Vorgehen von selbst verbietet. Man wird daher am ein- 
fachstenr, = (1, 0, 0,..0) oder ahnlich und s, = 1 wahlen, um méglichst viel Rechenarbeit einzu- 
sparen. 

Erstes Beispiel. Wir wahlen das System der Abb. 5. Dividiert man alle Elemente der 
Federungsmatrix durch c, die Elemente der Tragheitsmatrix durch m, so werden die dimensions- 
losen Matrizen 2{ und § nebst der Reziproken §—! 


2 —l —l1 3/4 0 0 4/3 0 0 
Dd |) Ail et 0 
—l —Il 3 0 


oOo = 
— © 
(or) 
i 
I 
= SS 
fom 
i=) 


0 1 


Wir beginnen die Transformation etwa mit r; = (1, 0,0), schreiben diesen Vektor links unten 
neben die Matrix 9 (siehe Hauptschema) und bilden den Vektor x/ 9{ = (2, — 1, — 1), der rechts 
neben r,, also unter 2 hingeschrieben wird. Rechts neben Q{ tragen wir den Spaltenvektor r, ein 
und ermitteln aus ihm und dem Zeilenvektor ry} { das innere Produkt , Ut, =e = 2-.Ferner 
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wird vB Y, = d,, = 3/4, und der erste Rechenschritt ist beendet. Die stark ausgezogene Linie mit 
der Ziffer 1 im folgenden Hauptschema grenzt alle Zahlen ein, die wir bis jetzt hingeschrieben 


haben. 


iI 2 
Yy Yo V3 
aes ies = 1 Ta ie ah ala 3/40 0° 
aN ae —1 2-1 UF earlae edo D={0 2 0 
5 a en | 3 0 i= 12 0 01/2 
px [0 0 ae i ee | 22 0\ @ 
Ac 1 s239 1 2 =—2 Sei 
NS pe es ES, Bye ee- =e) OF 2 07/4, 
Aus (15 a) ist nun Ce —r ABS 3+¢, rj/d,, zu ermitteln; den Vektor r, 2 finden wir bereits 


im Hauptschema vor, wir brauchen ihn nur noch mit 8—! zu multiplizieren, das hei®t hier insbe- 
sondere, seine erste Komponente mit dem Faktor 4/3 zu versehen und die beiden anderen wegen 
der Einsen in §~! einfach stehen zu lassen. Danach wird 


—r, =(=]2, 1, 1) 
—1, UB =(—838, 1, 1] 
Cee ee 7 Re = 
ues “Ua. OO 
Ua oa (aad UF "pay tee § 1) 
Jetzt tragen wir c unter a in das Schema ein, berechnen den Zeilenvektor ae = (— 2, 1, 2) 
und das innere Produkt r,t, = g,, =—2. Rechts neben r, ist der Spaltenvektor x, hinzu- 
schreiben; mit ihm werden die beiden Produkte r/ 1. = g,, =—2 und 1, %Ut, = go = 3 er- 
rechnet. Damit sind auch alle Zahlen innerhalb der mit 2 bezifferten Begrenzungslinie bekannt. 
Wir bemerken, da tatsichlich g,, = g,, = — 2 ist, die Matrix () muf ja symmetrisch werden. 


Will man das nicht als Kontrolle verwerten, so braucht man rj (1, = gj. gar nicht auszurechnen, 
sondern nur hinzuschreiben, da ja g,, schon vorher bekannt ist. Diese Bemerkung gilt natiirlich 
fiir jeden weiteren Rechenschritt. Nach (19) ist d,. =—g,, = 2. Der dritte und letzte Vektor 
wird nach (15a) 


fae Oe = (Serie 3/2ag* 3/2) (4 
J 


io oS 
Sie eee a, 
dy; ch Y, ( 8/3, 0 9 0) 
es 0, 1/2 , —1/2) 
Damit 148t sich das Hauptschema voll ausfiillen. SchlieBlich ist noch d,, = — g,. = 1/2, und die 


Matrizen (S} und ® sind berechnet. 
Das Gleichungssystem (25) wird nun 
2s, —2s, == (5 
—2s, +3s, —1/2s,=0. 
Wahlt man s, = I, so folgt leicht s, = 1 und s, = 2. Ks sind also lediglich die dritten Zeilen und 


Spalten von & und D nach (23) und (24) mit s, = 2 zu multiplizieren, das ergibt die gesuchten 
Matrizen © und § 


2 2 0 3/4 0 0 
0. —1 1 +6 0 0 2 
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Die zugehoérige Bildwelle zeigt die Abb. 6, die Seileckskonstruktion von Baranow-Schaefer Abb. 7. 

Die genauen Eigenwerte sind k, = (11 —/85)/6 = 0,2965, hk, = (11 + //85)/6 = 3,375 und 

ar ay k,; = 4,000 mit den LEigenvektoren 

a7 3 = (15 08970; 14); 4, = (1; 0,29; — 1,14) 

und 3, = (1; —1/2; 1/2). Die Trans- 
formationsmatrix (34) ist hier 


It @ 0 
0 RS = (Ler, 1-ry, 2-15) = Oa 
0 1—t!1 


Transformiert man zum Beispiel den 
Eigenvektor 4s zum Eigenwert k,= 4,000 
zuriick, so wird 4, = (ht © 3)’ = (1, 0, 
—1). 


a Ist die Tragheitsmatrix keine Diagonal- 
matrix, so braucht man nicht etwa die 
Reziproke $1 wirklich zu_berechnen, 
sondern nur %} in das Produkt 8 = OH’ 
zu zerlegen, wo % eine obere, also ©’ eine 
untere Dreiecksmatrix ist. Dann aber laBt 
sich r/,9( 8-1 ohne Kenntnis von $7! auf 


H=27 -——+ 


ee a 
HE, 27 €,[cm) a § [em] 


2 


Abb. 8. ,,Vollstandig‘* verzweigtes 
Abb. 7. Seileckskonstruktion nach Baranow-Schaefer fiir das System von Abb. 6. Wellensystem. 


einfache Weise ermitteln. Rechentechnische Einzelheiten hierzu findet man ausfiihrlich dargestellt 
bei Zurmiihl1, dessen bekanntem Lehrbuch auch die hier gewahlte Matrizenschreibweise ent- 
nommen ist, und in einer Arbeit des Verfassers 2. 

Wenn wir auch die allgemeinen Transformationsformeln hergeleitet und uns anhand eines Bei- 
spieles von ihrer Brauchbarkeit tiberzeugt haben, so bleibt doch, wie ein Blick auf die letzte Spalte 


J J J J 
Z Z Sete i” 
B\ AC ANC. C 


7- 3ma/ 
Abb. 9. Zerfallende Bildwelle des Originalsystems von Abb. 8. 


der Tabelle 1 zeigt, noch eine Frage offen: Die durchlaufende Welle, wie tberhaupt jedes Schwin- 
gungssystem der Klasse III hat lauter einfache Eigenfrequenzen, das vorgelegte Originalsystem 
der Klasse I aber kann durchaus mehrfache haben, wie z. B. das von Abb. 8 mit der (n — 2)-fachen 
Kigenfrequenz w? = @2 = ---- =w,—2 = C/J, ferner w, 1 4 w? und w2 = 0. Wie soll dann 
die zugehérige Bildwelle aussehen? Die Antwort gibt Abb. 9: die Bildwelle zerfallt in mindestens 


' R. Zurmihl, Matrizen. Berlin-Gottingen-Heidelberg 1950. 
2S. Falk, Abh. Braunschweig. Wiss. Ges. 6 (1954), S. 166. 
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n — 2 Teile, und zwar so, dai} jeder einzelne von diesen die mehrfache Eigenfrequenz des Original- 
systems als seinerseits einfache Kigenfrequenz besitzt, womit der anfangliche scheinbare Wider- 
spruch gelést ist. Ganz allgemein ist die héchste vorkommende Vielfachheit einer Eigenfrequenz 
des Originalsystems gleich der Mindestanzahl der Teilwellen im Bildsystem. Der Zerfall der Bild- 
welle kiindigt sich bei der Transformation von selbst an, und zwar dadurch, daf fiir irgendeinen 
Rechenschritt die Formeln (15a) den Nullvektor tp +1 = 0 liefern, dann verschwindet nach (19) mit 
8p +1, p auch c, , 1, p in der Federungsmatrix ©, und die entsprechende Drehsteifigkeit wird Null. 


Nach dem Zerfall wird die Transformation fortgesetzt mit einem Ersatzvektor ee 1, der auf 
allen bisher errechneten Vektoren Yj, - - -Yp verallgemeinert senkrecht stehen mu8; ein solcher 
Vektor ist zum Beispiel 


he , P tS Ff 7 
eo > j i tj (’ beliebig) . (35) 
JJ 


j=l 


Sollte auch hier wieder der Nullvektor herauskommen, so ist ein anderer Vektor fh’ einzusetzen. 


Zur Fortsetzung der Transformation brauchen wir weiterhin das Element dp 4.1,p+1, das man 
aus (19) natiirlich nicht gewinnen kann; denn (19) ist nur fiir solche Vektoren giiltig, die aus den 
Formeln (15a) stammen, nicht aber fiir den Ersatzvektor (35). Multiplizieren wir jetzt (35) von 
rechts mit § f, . 1, so wird 


Roe oe Opel pei) tps (36) 


und da fy’, weil frei wahlbar, im allgemeinen einfacher gebaut ist als ie 41, kommt man nach (36) 
bequemer zu d,,1,p +1 als auf direktem Wege. Fir d,+2,,42 usw. gilt dann wieder (19). Zer- 
fallt die Welle ein zweites Mal, so sucht man nach (35) wiederum einen Ersatzvektor, der auf allen 
bereits berechneten Vektoren r; verallgemeinert senkrecht steht, und so fort, bis die Gesamttrans- 
formation beendet ist. Der Zerfall der Bildwelle ist vom praktischen Standpunkt aus natiirlich zu 
begriiBen, denn die Behandlung zweier Wellenstrange mit n — p und p Freiheitsgraden ist erheb- 
lich einfacher als die einer einzigen Welle mit n Freiheitsgraden. 

Das Vorhandensein mehrfacher Eigenfrequenzen des Originalsystems ist fiir den Zerfall der 
Bildwelle zwar eine hinreichende, aber keineswegs notwendige Bedingung; wesentlich ist vielmehr 
die Wahl des Ausgangsvektors rj. Ohne die etwas langwierigen Beweise seien die genauen Be- 
dingungen zusammengestellt : 


Fall 1. Das Originalsystem hat lauter verschiedene Eigenfrequenzen. 


la) Der Ausgangsvektor ry, war so gewahlt, daB er auf keinem der n Eigenvektoren 7; des 
Originalsystems verallgemeinert senkrecht steht: dann zerfallt die Bildwelle nicht. 


1b) 1; steht auf n — t der Kigenvektoren y; verallgemeinert senkrecht; dann liefern die Formeln 
(15a) den Vektor 1;,; = 0; die Transformation bricht ab, und zwar hat die erste Teilwelle des 
Bildsystems genau t Freiheitsgrade. Je nach Wahl des Ersatzvektors f/, 1 kann die Welle noch 
weiter zerfallen. Eine zum Eigenwert k; gehérige Eigenschwingung geht stets folgendermafien von- 
statten: Alle Teilwellen ruhen bis auf die eine, die k; als einen von mehreren oder als einzigen 
Eigenwert enthalt. 


Fall 2. Das Originalsystem hat nur v <n voneinander verschiedene Eigenfrequenzen. Jetzt 
steht der Ausgangsvektor ry, von vornherein auf n—v Higenvektoren 1; verallgemeinert senk- 
recht, also bricht die Transformation spatestens an der Stelle v ab: die erste Teilwelle hat hochstens 
v Freiheitsgrade. Ist ti zu noch mehreren Eigenvektoren orthogonal, so hat die erste Teilwelle 
weniger als v Freiheitsgrade. Je nach Wahl des Ersatzvektors kann die Welle noch mehrmals zer- 
fallen. Die zu einem einfachen Eigenwert gehérige Eigenschwingung geht genauso vor sich wie 
unter Fall 1. Fiir einen mehrfachen Eigenwert k; mit der Vielfachkeit 6 aber gilt: die 6 Teilwellen, 
die k; als ihrerseits einfachen Eigenwert besitzen, schwingen gleichzeitig unabhangig voneinander. 
Insbesondere kénnen einige von ihnen in Ruhe bleiben, eine aber muB mindestens schwingen, da- 


mit auch das Originalsystem schwingt. 


ZusammengefaBt: hat das vorgelegte Originalsystem eine Eigenfrequenz von der maximalen 
Vielfachheit Sax, so zerfallt die Bildwelle in mindestens Onox Teile, wie auch der Ausgangs- 
vektor i gewahlt war. Sind aber die Eigenfrequenzen des Originalsystems alle voneinander ver- 
schieden, so kann die Bildwelle zerfallen je nach Wabl von r,, muf es aber nicht. 

22 
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Zum Schluf dieses Abschnittes machen wir noch eine technisch wichtige Bemerkung zur Schwin- 
gungsberechnung von Kurbelwellen: wie Grammel } gezeigt hat, treten bei Kurbelwellen auBer 
den Haupttorsionen noch sogenannte Nebentorsionen auf, d.h. die Federungsmatrix ist durchaus 
nicht die Dreiermatrix ©, sondern eine gewisse Vollmatrix Y{; die Kurbelwellen gehéren somit, 
wenn man die Nebentorsionen beriicksichtigen will, zu den Schwingungssystemen der Klasse I, 
nicht III. Nach der Transformation jedoch kann man sie wieder mit den iiblichen fiir Kurbelwellen 
ohne Nebentorsionen entwickelten Verfahren behandeln. 


Dp 
a, 
= 8 
le Gas 
Ms (RS)g 

Abb. 12. Die Transfor- Abb. 13. Die Transformations- 

mationsmatrix (Jt G)ges. matrix (Jt G)ges. einer zweimal 

einer einmal verzweigten verzweigten Welle. 

Abb. 10, Einmal verzweigtes_,Wellensystem, schematisch. Welle. 


ay 


eo 


Tp 

B 

NK Se 
Fay ete | 12, 
Up=LZy 


wh =27° ig Abb. 14, Die Federungsmatrix 2{ einer zweimal 
d — SS verzweigten Welle, schematisch. 
A \ 8 
i { ae 3. Die Abbildung durchlaufender Wellen in sich 


selbst. Wie jedes beliebige Schwingungssystem unter 

den Voraussetzungen (3) 1a8t sich natirlich auch 

eine durchlaufende Welle von n Freiheitsgraden auf 

B co"-fache Weise ihrerseits auf eine andere durch- 

St pe * laufende Welle, also gewissermafen ,,in sich selbst“ 
abbilden. Eine solche Abbildung erscheint auf den 

| | 7 | ersten Blick praktisch wertlos, liefert aber eine recht 


m 


A ia 


7 
| { brauchbare Methode zur Behandlung verzweigter 
Wellensysteme, wie sie etwa im Motorenbau haufig 
Abb. 11. Zweimal verzweigtes Wellensystem, schematisch. 2 £3 ° 
Maubiditet cas cine duconlauhenae Walle auftreten, Obwohl von Grammel2 fiir fast alle in 
der Praxis vorkommenden Falle ein einfaches Ver- 


fahren zur Berechnung solcher verzweigten Wellen entwickelt worden ist, moge die besondere Art 
der Abbildung hier kurz erlautert werden. 


Beginnt man die Transformation einer durchlaufenden Welle (allgemeiner irgendeines Schwin- 
gungssystems der Klassen II und III, auch I, sofern$ Diagonalmatrix ist) nach den Formeln (15a) 
mit dem i-ten Kinheitsvektor: 1; = ¢}, so wird die neue Koordinate y, gleich der alten Koordinate 


1 C, B. Biezeno u. R. Grammel, Technische Dynamik, Bd, 2, S. 349 und S. 410, 2. Aufl. Berlin-Goéttingen- 
Heidelberg 1953. 


> C. B. Biezeno u. R. Grammel, a. a. O. Bd. 2, S. 385. 
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*;, wovon man sich leicht tiherzeugt (siehe auch das gerechnete erste Beispiel). Diesen Umstand 
machen wir uns folgendermafSen zunutze: wir denken uns die einmal verzweigte Welle nach Abb. 10 
kurz vor der Verzweigung durchschnitten, setzen also vortibergehend C,_; , = 0. Nun bilden wir 
die durchlaufende Welle B so in sich selbst ab, daB die verknupfende Koordinate x, gleich y, wird, 
beginnen also die Rekursionsformeln mit t; =€, und hangen nach der Transformation die in sich 
umgeordnete Welle B mit Hilfe der Drehsteifigkeit C,_1,, wieder an den Wellenteil E an. Die 
Gesamttransformationsmatrix (}t ©)ges. hat dann die Gestalt von Abb. 12. mit der r-reihigen 


_Einheitsmatrix (,. Abb. 11 zeigt schematisch, wie man bei einer doppelten Verzweigung vorzu- 


gehen hat; die zugehérige Transformationsmatrix (Abb. 13) enthiilt in der Mitte die Einheitsmatrix, 
ferner die beiden Teilmatrizen ({ ©) 4 und (KR ©), die sich beim Umordnen der abgeschnittenen 
Teilwellen A und B ergeben. Daf diese Méglichkeit bei der doppelten Verzweigung iiberhaupt be- 
steht, liegt daran, da die Rekursionsformeln (15a) sich ja auch riickwarts lesen lassen, d.h. die Dreier- 
matrix (§ muf} nicht unbedingt von links oben nach rechts unten, sie kann auch in umgekehrter 
Richtung aufgebaut werden. Die Federungsmatrix Q{ des zweimal verzweigten Originalsystems 
von Abb. 11 hat die in Abb. 14 schematisch dargestellte Form. Es ware nun unzweckmaBig, die 
Transformation oben rechts oder unten links zu beginnen, weil dabei der ganze bereits dreireihige 
mittlere Teil der Matrix Qf iiberfliissigerweise in sich umgeordnet werden miiBte. Man laBt daher 
den mittleren fertigen Teil stehen — in (R ©)ges. erscheint fiir diesen Bereich die Einheitsmatrix, 
die alles ungeandert la8t — und beginnt an den Verzweigungsstellen die Transformation von innen 
nach aufBen, wie die Pfeile in Abb. 14 andeuten. Bei drei- und mehrfach verzweigten Wellenan- 
lagen kénnte man auf ahnliche Weise vorgehen, indem man verschiedene Zweige abschneidet, um- 
ordnet und wieder anhangt, und diesen Vorgang solange wiederholt, bis eine einzige durchlaufende 
Welle entstanden ist. Es ist indessen ratsam, bei einer mehr als zweimal verzweigten Welle die 
Transformation insgesamt durchzufiihren, was kaum mehr Rechenarbeit erfordert, dafiir aber auf 
jeden Fall iibersichtlicher wird. 


YT Tz D3 
Tei/A 47 ssmess Ze, Teil B 
Cb 


@, Ds 


Ca 
-24y a 
L @ | G ~ ra Kurbeltriebe G | C | 
a J 2d ay en Uy) 


Abb. 15. Zweimal verzweigte Wellenanlage. 


Zweites Beispiel. Es sind die Eigenschwingungszahlen der doppelt verzweigten Maschinen- 
anlage nach Abb. 15 zu ermitteln. Man kénnte hier Federungs- und Tragheitsmatrix des Gesamt- 
systems aufstellen und die Transformation in eins durchfiihren. Es ist aber zweckmaBiger, die 
soeben hesprochene Methode anzuwenden; wir denken uns die Wellenstrange der Drehsteifigkeiten 
C, und C; entfernt und betrachten jeden der abgetrennten Teile A und B fiir sich. Fir den Teil B 
werden mit den Gréfen der Abb. 15 in iiblicher Weise Federungs- und Tragheitsmatrix zusammen- 


gestellt, und zwar ist 


8 2 0 2 0 20° 02.0, 0° 0 
2 2 —1 0 0 Umer Tee thee 20 
%e=|0 —1 1 0 On me aa Oo Ol Ona) 
2 0 0 2 —1 One0s Ont 0 
0 0 0 —l1 1 O00), 0 I 


Der Ausgangsvektor 1, ist natiirlich nicht beliebig — er ware es nur dann, wenn die Transformation 

fiir das ganze System insgesamt durchgeftihrt wiirde — sondern wir miussen r, = (1, 0, 0, 0, 0) 

wahlen, damit die neue Koordinate y, gleich der alten x, wird; denn sonst kénnte man die in sich 

abgebildete Welle B nicht einfach wieder an das Mittelstiick E anhangen. Nach den Rekursions- 
22* 
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. . / / a 
formeln (15a) werden nun die weiteren Vektorenr,. ..1;: 


dy, = 20, —t; %s 8; * 299/20 2 0 Ome A 
8-0 //dy = 8/20 r; En so 0 0 er 

16 = 0 —2 2 0 

dyy = — 8 = 8, — 1, UnB;z! = 19/000 aes e yeaa, 
(WN Lg teat = 0 —4 0 —4 0 ae 

B10 °t, [dy = — 8/201, a Opa 0 0 0 0 

Ls = 0 0 —2 0 —2 

dg = — 839 = 8, —t1, UsB," == QO —2 a 2 
Se, ade = 0 0 —2 0 —2 -{- 

Bog Vq/4yg = — 1-1, = 0 2 he 2 0 

cs = 0 0 0 0 0 
Als Ersatzvektor errat man hier leicht, wenn man nicht Formel (35) benutzen will, etwa t= 


1,1, —1, —1). Es wird dann d,, =f; Of, = 4, und damit fahrt man fort: 
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dy =4(#—éa =) ¥ rie By? = 0 1 0 
EO ie OUP. = 0 1/2 1/2 —1/2 —1/2 
da 6,7, 1. ce “= 0 —l1/2 1/2 1/2 —1/2. 
Das Hauptschema ist 
We | 8 2 0 2 0 il 0 0 0 0 
2 2 —I1 0 0 0 —2 0 1— 1/2 
| 0 —1 1 0 0 0 0 —2 1 1/29 
| uy 0 0 2 —l1 0 —2 0 —l1 1/29 
Ry | 0 0 0 —1 1 | 0 0 —2 —1—1/2 
(1 0 0 0 0 8 2 0 2 0 8 —8 0 ; 0 0 
| 0 —2 0 —2 0 —8 —A4 2 —4 2 — 8 16 —8 0 0 
0 0 —2 0 —2 0 2 —2 a 0 —8 ome 0 0 
; 1 1 ee he =D 0 1 Qnty Melee: - at ae acces car E~ a 
(OR ee ee 2 O32) Sele ce 0 0 One 5/2 


Die dreireihige Teilmatrix oben links in pg hat bereits die Endform (4); denn alle Zeilen- und 
Spaltensummen sind gleich Null. Die zweireihige Restmatrix unten rechts jedoch stellt in der 
jetzigen Form die Federungsmatrix einer versetzten eingespannten Welle dar, man sieht aber ohne 
Rechnung, da nur die fiinfte Zeile und 


ae Spalte von Gg mit dem Faktor 2 multi- 

ef bu pliziert zu werden braucht, damit auch die 

4J 4 vierte Zeilen- und Spaltensumme gleich Null 

yas wird; es stimmen also alle Elemente c;, 
(C;) 8 | 8C ZC. 8c Y undi;; mit den g;, und d,; iiberein auBer 
Yy Ys C54 = C45, C55 und t;,, und zwar ist nach (23) 

Ys 2, und (24) ¢54 = Cys = 2 854 = —2, C55 = 27g,5 

yn By = 10 und i,, = 27 d,, = 4. Fir 10 schrei- 


Abb. 16, Bildwelle des Teiles B des Originalsystems von Abb. 15. 


nach (4). 


ben wir 2-+8; die Federkonstante des 
eingespannten Wellenstiickes ist dann 8 C 


Abb. 16 zeigt den transformierten Wellenabschnitt B. Da nun der Teil A des Gesamt- 


systems die gleiche Form hat wie Teil B, so liefert er auch die gleiche Bildwelle, womit sich die 


) 
) 
) 
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zweite Transformation ertibrigt. Das Gesamtbild des Originalsystems ist somit das der Abb. 17. Die 
beiden eingespannten Teilwellen ergeben nach leichter Rechnung oder auch mit Hilfe der an- 
gedeuteten Baranowkonstruktion (Abb. 18) jeweils die Eigenwerte (3 —y5)/2 und (3 + )5)/2, 


so da®B also vier Eigenschwingungszahlen der Gesamtanlage bekannt sind unabhangig von der 
Beschaffenheit des mittleren Teiles E! 


ie 


7 Kurbetriebe 
FY YY SS X Sd, AS re Se See 
A> A; G, &, 


Abb. 17. Bildwelle des Gesamtsystems von Abb. 15. 


4. Die Abbildung einer Reihe von Elementarschwingern auf eine durchlaufende Welle. Fede- 
rungs- und Tragheitsmatrix eines Systems der Klasse IV sind beide von Diagonalgestalt; die 
Eigenschwingungszahlen lassen sich als wm? = C;/J; unmittelbar ablesen. Die Transformation nach 
(15a) und (25), die jetzt kaum noch nennenswerte Rechenarbeit erfordert, gestattet somit, eine 
durchlaufende Bildwelle mit vorgeschriebenen Eigenschwingungszahlen zu konstruieren, eine Még- 
lichkeit, die in der Mechanik nur theoretischen Wert zu besitzen scheint, in der Elektrotechnik je- 
doch von einiger Bedeutung ist. Auch H. Schaefer! hat hierfiir eine Methode angegeben. 


5. Verallgemeinerungen. Von einigen der bisher gemachten Voraussetzungen wollen wir uns 
zum Schlu®B wieder befreien: 

1. Die Federungsmatrix 2{ des Originalsystems sei indefinit; einige Werte k; = w? C/J sind dann 
negativ, in den Lisungsfunktionen treten somit Ausdriicke der Form oj «;t und @in @; t auf: die 
Schwingung ist nicht mehr stabil. Die dann 
ebenfalls indefinite Matrix ® berechnet man ; 
nach den Rekursionsformeln (15a) wie sonst [om] al 
auch. Vor der Berechnung der Gréfen s; 
aus (25) jedoch fihrt man eine sogenannte 
,. Spektralverschiebung“ der Eigenwertskala 
durch, d. h. man bildet die neue Matrix 


&%* =G+tD (¢ positiv reell). (37) 


Die Eigenwerte des, Paares (§*,®D sind 
simtlich um t gréBer als die des Paares 


, D bzw. Y,B, also nicht negativ, wenn vA 
man nurt grof genug wahlt, z. B. sicher 4 2 


dann, wenn 


beet | hgh | sage 0 ay 74); (38) 


4J 


> 


nach einer bekannten Abschatzung fiir H=4J S| 
den maximalen Betrag von Eigenwerten. 
Nach (38) bildet man zunachst die MatrixS), 
indem man die Zeilen von ( der Reihe ie aC aoe 
nach durch d,,,... dyn dividiert, sucht > HE. ZTE, fom) °F Elem) 
das dem Betrage nach gréBte Element Abb. 18. Seileckskonstruktion nach Baranow fiir die Teilwelle B, von Abb. 16 
in § heraus und multipliziert es mit der oe 
Ordnungszahl n der Matrix (. Dieser Wert wird zweckmafig nach oben zu einer ganzen Zahl auf- 
gerundet. 

2. Die Federungsmatrix 9 sei negativ (halb)definit. Es handelt sich also nicht mehr um ein 
Schwingungssystem im Sinne des Wortes, da keiner der Werte w? positiv ist. In diesem Fall rechnet 


1 H. Schaefer, Vorlesungsmanuskript 1954. 
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man einfach mit dem Paare —Q(,%, dessen Eigenwerte gleich —k,, also samtlich nichtnegativ sind, 
und verfahrt wie sonst. 

3. Die Matrizen 2{ und % seien jetzt komplex Hermitesch, 8 auSerdem positiv definit, wie es 
zum Beispiel bei den Anwendungen in der Elektrotechnik haufig vorkommt: 


xy —A, B=%’, BG positiv definit. (39) | 
Der Querstrich bedeutet Ubergang zu konjugiert komplexen GréBen. Alle bisher abgeleiteten | 


Formeln und Folgerungen gelten wortwortlich auch hier, wenn man unter den Bilinearformen g;, 
und d;, die mit den konjugiert komplexen Vektoren t; gebildeten Ausdriicke 

fy, == 0, 2%, und dy Tt, oY, (40) 
versteht. Die Hauptdiagonalelemente d;; und g;; sind als Hermitesche quadratische Formen stets 
reell, somit sind es nach (19) auch die Nebendiagonalelemente von (}. Die Matrix {i dagegen wird 
im allgemeinen komplex sein. Da nun schiefhermitesche (im Reellen schiefsymmetrische) Matrizen 
durch Multiplikation mit der imaginaren Einheit i in Hermitesche iibergehen, rechnet man mit den 
Paaren il, % (ik reell) oder i, iQ (k reell), je nachdem, ob nur eine oder beide Matrizen des 
Paares Y{,%8 schiefhermitesch sind. 

4. Sehen wir schlieBlich ganz von der physikalischen Bedeutung der Matrizen 2{ und § ab, so 
geniigt es, wenn entweder Q{ oder  negativ oder positiv definit ist; man braucht ja nur 1/k bzw. 
— 1/k als neue Eigenwerte einzufiihren. 

Zusammenfassend laBt sich dann sagen: Die Transformation des gegebenen Paares Y(, 8 auf das 


reelle Paar ©, § und somit die zeichnerische Lésung des Kigenwertproblems nach Baranow-Schaefer | 


gelingt stets, wenn die folgenden beiden Bedingungen erfillt sind: 

a) Jede der beiden Matrizen des Paares ist entweder Hermitesch oder schiefhermitesch (insbe- 
sondere im Reellen symmetrisch oder schiefsymmetrisch). 

b) Eine der beiden Matrizen des Paares ist positiv oder negativ definit. 

Die Voraussetzung b) laBt sich im allgemeinen nicht entbehren, da die quadratischen Formen 
in den Nennern der Rekursionsformeln (15a) nicht Null werden diirfen; doch ist es durchaus még- 
lich, da bei gliicklicher Wahl von ry, die Transformation wie im Regelfalle vonstatten geht, auch 
wenn keine der beiden Matrizen 2{,9 definit war. 

5. Bislang haben wir nur von der freien Schwingung mit der Bewegungsgleichung 


St + Wr = 0 (41) 


gesprochen. Oft lautet die Aufgabe allgemeiner, namlich 


BE-XE Ad (42) 
mit beliebiger rechter Seite p. Die Komponenten von p sind z. B. konstant, wenn die gewahlten 
Koordinaten x; fiir die statische Gleichgewichtslage des Systems Null sind, diese aber nicht mit 
der entspannten Lage zusammenfallt (sogenannte »,Vorspannungen’). Greifen auBere periodische 
Krafte oder Momente am System an (,,Erzwungene Schwingungen‘), so sind die Komponenten 


von ) Funktionen der Zeit, gegebenenfalls auch des Ortes usw. Fiihren wir nun in (42) die Trans- 
formationsvorschrift (20) ein, so wird 


BRSZ4ARSGS3Z=—p. (43) 
Dies von links mit @’ 8h’ multipliziert ergibt 
SWBRSILSG NARS =S' Hp, (44) 
oder nach (11), (23) und (24) 
Sj +03=S Rp, (45) 


und das heifit: der Vektor ) am Originalsystem transformiert sich in den Vektor @’ SR’ an der 
durchlaufenden Bildwelle. Da nach der Transformation ©’ 3’ p leicht zu berechnen ist, lassen 
sich damit auch sdmtliche fiir die inhomogene Schwingungsgleichung der durchlaufenden Welle 
entwickelten Verfahren auf die inhomogenen Systeme der allgemeineren Klasse I ausdehnen, 


(Eingegangen am 23. Oktober 1954.) 
Anschrift des Verfassers: Dr.-Ing. Sigurd Falk, Braunschweig, Celler Str. 93a. 
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Bestimmung von Balkenschwingungen mit Hilfe des Matrizenkalkiils * 
Von H. Fuhrke 


1. Einleitung, Bezeichnungen. Wenn man von komplizierten mechanischen Systemen die 
hoheren Eigenfrequenzen bestimmen will, mu8 man das fiir die Rechnung beniitzte Ersatzsystem 
dem wirklichen sehr weit annahern. Ein Mittel, diese Anniherung bei geeignet gegliederten Trag- 
werken sehr weit zu treiben, ist der Matrizen-Kalkiil, der es erlaubt, die Rand- und Ubergangsbe- 
dingungen besonders tibersichtlich zu formulieren. Dabei kénnen, wenn das Tragwerk aus Balken 
besteht, auBer Biegeelastizitat und Massentragheit (Ziff. 2) auch Schubelastizitat und Rotations- 
tragheit beriicksichtigt werden (Ziff. 3). Eine Unbequemlichkeit, die auch der gewéhnliche Matrizen- 
Kalkil nicht vermeidet, daB namlich bei ausgearteten Systemen (sehr grofe Massen, sehr steife 
Stiitzen) in der Enddeterminante Differenzen groBer Zahlen auftauchen, laBt sich durch den Uber- 
gang auf Determinanten-Matrizen (Ziff.6) vermeiden. Auch Eigenschwingungsformen (Ziff. 9) 
sowie die Verformung und Beanspruchung bei erzwungenen Schwingungen (Ziff. 10) und unter 
statischer Last (Ziff. 11) lassen sich mit Hilfe der Matrizen ermitteln. Ein Beispiel, die Bestimmung 
der Biegekritischen eines Maschinenlaufers, zeigt den Rechnungsgang (Ziff.12). Uber die Er- 
weiterung der Methode fiir ebene und raumliche Rahmentragwerke soll ein zweiter Teil berichten. 

Den positiven Achsen x, y, z eines Rechtskoordinatensystems seien positive Verschiebungen u, 
v, w und der positiven x-Achse ein positiver Torsionsdrehwinkel y zugeordnet. Schnittkrafte seien 
dann positiv, wenn ihre Vektoren und die AuBennormale nz der betreffenden Schnittflache beide 
in positive oder beide in negative Koordinatenrichtungen fallen. Ein homogenes Balkenstiick sei 
als Feld bezeichnet, mehrere Felder, z. B. zwischen zwei Lagern, seien zu einem Abschnitt zusam- 
mengefaBbt. Auf die Stelle i folge in x-Richtung das Feld i oder der Abschnitt 7. 


Feld: i mit i; Ty, my =u; 1; 


cag 


Stelle: (i +1) @ 


Le 
az 
Abb. 1. Abb. 2. | 


2. Matrizengleichung fiir die Querschwingungen eines homogenen Balkenstiicks. Aus den dyna- 
mischen Grundgleichungen ys iv = Q’, M’ = Q fiir das Balkenelement und dem Elastizitatsgesetz 


w’’ = — M/ET fiir den gebogenen Balken ergibt sich die Differentialgleichung fiir die Schwingungen 
des homogenen Balkens: EI w'Y = — wv wo, die mit dem Eigenschwingungsansatz w(x, t) = w(x) e*” 
iibergeht in 4 ) 

; wv —(4) ie 0, mit aon, (1) 


wo w nunmehr nur noch ortsabhangig ist. Hierbei bedeutet 1 Masse der Langeneinheit, | die Linge, 
E den Elastizitatsmodul und J das Flachentragheitsmoment des Balkenfeldes. Schubelastizitat 
und Rotationstragheit sind zunachst nicht beriicksichtigt. Die vier Konstanten des Lésungs- 
ansatzes zu (1) lassen sich durch die vier mechanischen GréBen 


Verschiebung w, Neigung y = w', | 


2 
Moment M —— EI w"’, Querkraft Q = — EI w'” | (2) 


ee yee 
am Randi (x =0, Abb.2) ausdriicken, wenn man die vier Grundlésungen Wo} 7 * om 7 


anu 2 oc 
x, sin * geeignet kombiniert : 


1 A A Ie eit FE Pee 
w(x) = w; 5 (sof 7% + cos 5 x) + yl (gin 7 * + sin = x) 


eG 

"EI 222 

* Auszug aus der Dissertation des Verfassers, Darmstadt 1953, Referent: Prof. Dr.-Ing. K. Marguerre, 
Korreferent: Prof. Dr.-Ing. K. Karas. 


A 
cos a 


F p (3) 
i p a ee ; 
(sj 7 x cos 5 x) Q; EI 973 (gin ie Bip — fey Fi x] . 


—M 
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Durch Ableitung erhilt man hieraus die iibrigen GréBen y, M und Q an der Stelle x. Setzt man 
x =1 und fihrt die Funktionen! 


M4 28 j2 
C=5 (soja + cosa) =1 Ericeries re 


1 ; : Aft a8 
s =} (Gina + sind = Ey a gee ee 
1 A Re 
¢ = 5's (Gof — eos) = a GI gh ta 
: : 1 MM a8 
s = 515 (Gina—sina) = ~+— +a te 


ein, so erhalt man die Gleichungen 


PM, LQ, 
wisi = Cw, +S8ly;+ e(— ar)‘ s(— =I. 


5 
EM. ; Pm), 5(_Pa) ®) 
s0L6, 5 Br Mew +A sly; + = EI | EI)” 
OOha 


PM; 18Q, 
— iT =i Su, +Bely+ats(— a | c( oe 


Dieses Gleichungssystem verbindet die vier (dimensionsgleichen) mechanischen GréBen w, 1 y, 
a M, ae Q am linken Rand mit denen am rechten Rand. Wir kénnen es als Matrizengleichung 
schreiben und wollen es dabei zugleich insofern verallgemeinern, als wir zulassen, dal} die Bezugs- 


groBen lund EI nicht die des gerade betrachteten Feldes sind ?. Nennen wir 


rs = RP =e B a 
0 = 4 — —— VM, = M,, — (= C- 6 
Ly; Yi» EI EI Q; Q; ( ) 
und fihren 
EI. Le 
— = i? —— a 7) 
EI S l P ( 
ein, so kommt 
is a ie 2 3 | 
w C pS = oe w 
2 
» os C gE le | oi : 
= <a Ata Aha s M ( ) 
M Ge mar S p 
= Ma Ma Ma = 
Ree ee, rma pe © foe g Ji Or: 


Symbolisch kann man dafiir schreiben 
hina] = [ Ri eae (8’) 
indem man die vier GréBen w, yp, M und Q als Komponenten eines (Spalten-)Vektors \) auffaBt. 
Die Funktionszeichen C; usw. in (8) bedeuten > Goj A; + +++ mit 4 = 02: H/ET; . 
Die Matrix {i sei im folgenden kurz als Rayleigh-Matrix bezeichnet, da sie sich wesentlich aus 


den Rayleighschen Funktionen aufbaut. Sie ist symmetrisch zur Nebendiagonalen; der Wert ihrer 


* Die Funktionen (4) sind, von den /-Potenzen im Nenner abgesehen, als Rayleighfunktionen bekannt; 
vel. K. W. Wagner, Schwingungen und Wellen, S. 373, Wiesbaden 1947. 

> Es ist durchaus méglich, EI; und |; als BezugsgroBen zu verwenden; man mu dann aber bei mehr- 
feldrigen Systemen besondere Ubergangsmatrizen einfiihren [s.-W.Schnell, Z. angew. Math. Mech. 35 (1955), S. 272.] 
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Determinante ist | K| = 1. Sie ist gleich ihrer Kehrmatrix Si{-! bis auf die Vorzeichen, die schach- 
brettartig verteilt sind, wobei das erste Element (links oben) mit + beginnt!. 


Die Funktionen (4) lassen sich durch gut konvergierende Reihen darstellen; fiir 2 —1 und 
8 Stellen Genauigkeit reichen z. B. die 3 ersten Glieder aus. 


3. Matrix fiir Balkenquerschwingungen bei Mitberiicksichtigung yon Schubelastizitat und Rota- 
tionstragheit. Wir wollen die Matrix { fiir Balkenbiegeschwingungen noch erweitern auf den Fall. 
da Schubelastizitat und Rotationstragheit nicht vernachlassigt werden diirfen. Beide Einfliisse 
haben zur Folge, dafs in der Differentialgleichung auch die zweite Ab- 
leitung w’’ auftritt, und da8 infolgedessen die charakteristische Glei- 
chung zwei Lésungen J, und A, aufweist, so daB wir nicht mit den ein- 
fachen Funktionen (4) auskommen. Dagegen bleibt der iibersichtliche 


Aufbau der Matrix } erhalten. ( 
Aufere und innere Dampfung lieBe sich durch eine komplexe 
»»Masse‘ ww (1 —ik,/m), einen komplexen ,,E-Modul* E (1 + iw k,) \ 


und einen komplexen ,,Gleitmodul** G (1+-iwk,) erfassen?. In der 7’ 
folgenden Herleitung sei die Dampfung aber nicht mitgenommen. 
Bei der Aufstellung der Differentialgleichung miissen wir unter- 
scheiden zwischen der Neigung wy des nach der Verformung schief Abb. 3. 
stehenden Balkenquerschnitts gegen die Vertikale und der Neigung 
w’ der Balkenmittellinie gegen die Horizontale (Abb. 3). Ihre Differenz (w’ — y) ist die Schub- 


verformung. Die Elastizitatsaussagen lauten dann? 
—M=EIy’, Q=GFs (w’—y). (9) 
Die Bewegungsgleichungen sind 
ww=Q, pwiap=Q—M’. (10) 
Fiir die harmonischen Schwingungen haben wir also die vier symmetrisch gebauten homogenen, 
linearen Differentialgleichungen 


[ee Ce: / ht, 
ey 0 Re i eGR eI Py (11) 
QW =—auw, M’—Q=a' psy, 


die man unmittelbar durch den e-Ansatz lésen kénnte. Es ist aber tibersichtlicher, drei der Unbe- 
kannten zu eliminieren. Zundchst erhalten wir, wenn wir die erste und die dritte Gleichung (11) 
in die einmal differenzierte zweite Gleichung (11) einsetzen, 


1, Ob 1 ee: 


genau so ergibt sich, wenn wir die vierte Gleichung (11) differenzieren und die erste und die dritte 
einsetzen, 
wo wig 
M’ + —7-M +o? pw =0. (12 b) 
Aus den beiden Gleichungen (12a) und (12b) folgt* schlieBlich durch Elimination von M die Differen- 
tialgleichung fiir w 


2 52 
ry 4, 22H (EL ) a ee ee) = 0. 13 
ey (Gr +i um ii eT okey Me (13) 
Mit J (1), wobei wir statt A hier lieber A, schreiben wollen, sowie mit 
I \e 72 
5 ——=! he Aa und pa? (14) 
PGi OF, al P 


1 Wegen weiterer Eigenschaften dieser Matrix vgl. die Dissertation von W. Schnell, T. H. Darmstadt 1954 
,Krafteinleitung in versteifte Zylinderschalen“. ZFW 1955. ; 

2K. W. Wagner, a. a.O., S. 262 und 370; L. Cremer, Vierpoldarstellung und Resonanzkurven bei 
schwingenden Staben, SB. d. pr. Akad. d. Wiss. Phys.-Math. Kl. 1934. I. ws 

3 Vel. z. B. K. Marguerre, Neuere Festigkeitsprobleme des Ingenieurs, S.57, Berlin, Gottingen, Heidel- 
berg 1950. 

4 Z.B. S. Timoshenko, Vibration Problems in Engineering, S. 338, New York 1937; W. Prager- u. 
K. Hohenemser, Dynamik der Stabwerke, S. 118 f., Berlin 1933. 
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lautet (13) 
wiv 4 Ab B00) uw” a le (15) 


Dieselbe Gleichung erhalt man auch fiir y, M und Q. 

Laf®t man den SchubeinfluB und die Rotationstriigheit beiseite, setzt also G's = co und ip = 0, 
so wird 6 = 0 und } = 0, und aus (15) wird die Gleichung (1). Nur einen der beiden Kinfliisse zu 
vernachlassigen, bringt mathematisch keine wesentliche Vereinfachung. 


Der e-Ansatz w —c,e | fiihrt auf die charakteristische Gleichung 


M4 +24 (6 + 8) A? —At (1 — A156 9) = 0, (16) 
deren Wurzeln + J, und + iA, sind mit 
oe ae Se A (6 — 0) + = 12 (6. 4.8). (17) 
Zwischen den Wurzeln bestehen die Beziehungen 
22 —}? =/)8 (6 + 8) und "2? 2 = (lA 00) (18) 


Die weitere Aufgabe besteht nun darin, die vier Integrationskonstanten des e-Ansatzes durch 


die vier Randwerte von w, y, M und Q auszudriicken. Das wird am einfachsten, wenn wir von dem 
Ansatz 


Q =, Gof A, a + c, Gin A, + + cg cos daz + ¢, sin daz (19a) 
fiir die cas ausgehen. Aus der dritten Gleichung (11) erhalten wir 
i oa Q’ = as A, Sin A, — 7+ CoA CofA > 2h ond seins are, Aecos As i ae (19b) 


wobei 1 —1 und EI = EI die speziellen GréBen des betrachteten Feldes sein sollen. Aus der 
zweiten Gleichung (11) kommt 


y=lw' +6Q0= alee + A) (e, Gof a * +c, Sind, 4) 
— (A? — 6 Ad) (c cos A, — + c, sind, =\|, (19¢c) 
schlieflich liefert die erste Gleichung (11) 
WA ee Lapis 7 G + dA‘) A, (« Sin A, + + cy Co} A, 3 


(Ai — OA Ae (: sin Ay | — ¢4 608 Ay =| (19d) 
Wenn wir in den Gleichungen (19a) bis (19d) « = 0 setzen, erhalten wir 
ae = 1 2 : 
Qi =o +45 » Vi ia [(At + OA) 1 — (Az — 644) ¢s] » 
; Oe , (20) 
Ce 7 As Cy + Ap Cy) » M; = 7a [Ai + 6 AG) Ay 2 — (AZ— 6 Af) Ap ey] - 
Wr fithren noch die Abkiirzungen ein: 
= 1 = 1 
8 AFAR Ba EE — OF 
eae 12(6 — 9) 
is t= 2] meee 
1 EEA ES 5) * Tao) A (21) 


Ad _ 1, ~_ BO) 
PPL 2 ce V4+a(6— 9 |’ 
zwischen denen die Beziehungen 


Apa A, ls undsy AAG aes. (22) 


1 
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bestehen, und erhalten 
¢ =MA Yi + AL Q;, 6, =A! A, y, + A, 0;, | 
ak ABS ae as ay 23 
= 7 Arm +7 AM, , = 9 Ai — 7 Ay Mi. | ) 


Setzen wir endlich diese Konstanten in (19 a) bis (19 d) ein, so entstehen daraus fiir x =I vier 


Gleichungen zwischen den GréBen w, wy, M und Q am rechten Rand i und denen am linken Rand 


sea 


Die so gewonnenen Gleichungen lassen sich wieder zu einer sehr iibersichtlichen Matrizen- 
gleichung zusammenfassen, wenn wir die folgenden neun Funktionen cinfiihren; 


fir 6=0 und 9 =0: } 
Cy = A, Goj A, + A, cos, | 
Cyr = A, CofA, + A, cos, | 
Spo is o(A, Sin A, ne sin J) 


> (Gof Ay + cos Ay) , 


Sites i, A ind, io — = sin A, —> oe (Git Ay + sind) , 


Ales 
Sim = 7, 2 ind, 2 ae sit Ay (24 


: Rene as 3 73 (Sof Ay — 05 4) 


1 : 
> FF (Sint Ay — sin Ay) . 


ST area 205 e As Sin Ay =, 71 1 sin ho) ] 
Sey 0( 3, Gin 4, —7-sin 2) yl 
J 


SIL = 4a ig A Ay Gin A, —A, A, sin A,) | 
Mit (24) erhalten wir dann die Matrizengleichung fiir die Querschwingungen eines homogenen 


Balkenfeldes unter Beriicksichtigung von Eigenmasse yw und Rotationstragheit 1 i, Biegesteifigkeit 
EI und Schubsteifigkeit GF,: 


Tw | Cy Sr c sp | | | 
YY | Ao sit Ci Sir c y 
F As — | 4 4 : =e 9 (25) 
M Ab € Agsmr Cig Sr M 
PO eg (AG Stir Ap € Acre Cr 4 Q aR 
oder kurz 
[Yi4i] = [Us] ye. (26) 


Die Matrix 1 ist ebenso wie { in (8) symmetrisch zur Nebendiagonalen. Ihre Determinante ist 
1, ihr Kehrwert entsteht durch Vorzeichen-Umkehr bei den Elementen mit ungerader Zeigersumme. 
Bei Beziehung auf ein ,,fremdes‘‘ EJ und | treten zu den einzelnen Elementen dieselben a- und f- 
Potenzen wie in (8). Als Sonderfalle erhalten wir 


a) fir yw; = 0, d.h. fiir das masselose Balkenfeld, mit A, = 0, die Elastizitatsmatrix, gleich 
erweitert im Sinne von (7) und (8): 


cea : (27) 
Oueae ll B 
BOs. 0 1 Ji 
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b) fir 1; > 0, wobei aber gleichzeitig 4; 1; > mj gehen soll, d. h. fiir das zu einer Einzelmasse 
zusammengeschrumpfte Balkenstiick die Massenmatrix: 


if 0 0 0 | 

0 {I 0 0 28 

Te Wie 2 ornament 0 oo 
im; Ww? B/EL 0 0 lL; 


(Die Matrizen © und Nt lassen sich auch unmittelbar aus den Elastizitatsaussagen und den dyna- 
mischen Grundgleichungen fiir den Balken mit Einzelmasse finden.) 


4. Abschnittsmatrix. Besteht ein Balkenabschnitt aus mehreren homogenen, jedoch vonein- 
ander verschiedenen Feldern und werden an den Ubergangsstellen i + 1 bis k —1 keine auferen 
Krafte z. B. durch Lager, eingeleitet, so gilt an jeder Ubergangsstelle \); links= vj rechts, und man 
erhalt mit (19) die Abschnittsmatrix 2{ nach den Regeln der Matrizenmultiplikation? (Abb. 4), 


aus 


(C] = [Mea]... [Wear] (WT, (29) 
womit [ye] = [20]: [92] (30) 
Feld:  (k-1) + (irt) i fous ; 


0; Y ae 7H 
Stelle: (k-7) --- (U7) 4 Y) fi 
Abschnitt | Stelle: (+1) 
Abb. 4. Abb. 5. 


wird. Hiernach la8t sich z. B. aus (27) und (28) die Matrix & fiir ein Balkenstiick gewinnen, das 
aus einem elastischen, masselosen Feld und einer anschlieBenden konzentrierten Einzelmasse be- 
steht? (Abb. 5), als Summe zweier Matrizen geschrieben: 


0 0 0 0 
0 0 0 0 
ee f B3 Bi : 
g; = MN, &; ae a TS 0 —p Do St See +. ¢;. (31) 
pe pe /1 
ei clcmese 


Die Abschnittsmatrix fiir einen Balken mit beliebiger Massen- und Steifigkeitsverteilung kann man 
naherungsweise bekommen, wenn man den Abschnitt in n Felder teilt und die n SR- oder S- 
Matrizen nach (29) ausmultipliziert. In beiden Fallen muf§ man mit einer mittleren Steifigkeit 
rechnen*. Fiir die $i-Matrix ist 4; = m;/l;, fiir die T-Matrix ist my = 1/, (m; + mj.1) zu nehmen. 

Die Z-Matrix hat den Vorzug, daB sie einfacher gebaut ist als die $i-Matrix, insbesondere lat 
sich der Kigenwert w* abspalten, (31), was fiir das numerische Rechnen sehr vorteilhaft ist. Lange 
homogene Abschnitte miissen aber in mehrere Felder unterteilt werden, wenn man ausreichende 


Genauigkeit erzielen will, waihrend die {-Matrix beliebig lange homogene Abschnitte exakt er- 
fabt. 


5. Randbedingungen, Eigenwertdeterminante. Gleichung (30) lautet fiir einen Balken mit den 
Randern | und k ausfihrlich 


oy Ay, Aj, Aj Ayy w 
yp = A» Ay, Ag, Ay, y . (32) 
M As, Az» As, As, M 

LQ Ais ik Ay Ago Ay; Ag, | LQ jl 


* R. Zurmihl, Matrizen, Berlin, Gottingen, Heidelberg 1950. 


° W.T. Thomson, J. appl. Mech. 17 (1950), S. 337, wo diese Matrix mit 6 = 0 und @ = 0 als Ausgangs- 
punkt dient. 


° Siehe z. B. W. T. Thomson, a. a. O., Area Moment Principle. 
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Von den je vier Komponenten der Vektoren ), und {), lassen sich je zwei aus den Randbe- 
dingungen bestimmen. Bei den vier Lagerungsarten der Abb. 6 sind je zwei Komponenten Null. 


Setzt man die jeweils geltenden Randbedingungen in (32) ein, so kann man diejenigen zwei 
homogenen linearen Gleichungen 
herausziehen, fiir welche die linke 4 hl 
Seite Null ist. So ergibt sich )}-— ere 3) 4) 
z. B. fiir den links eingespannten, 
rechts gelenkig gelagerten Ral- Finspannung Gelenk Vertikaltihrung Frej 
ken 


0 | 0 
0 Y, 
Me == [2]- ce 
Qk. | L Q1 
vet y. S| 
eee ee ae) | P| 
a ie a | s | =) See 


Das Gleichungssystem (33) hat nur dann eine Lésung fiir p, und Q,, d. h. der Balken fiihrt nur 
dann freie Schwingungen aus, wenn seine Determinante verschwindet. Die Frequenzgleichung 
lautet also 
Aj, Ay 

22 “toa 

Ist Q{ noch geschlossen darstellbar, so kann hieraus die Eigenfrequenz bestimmt werden. Ist 

z. B. fiir den homogenen Balken 9{ = St, so wird im Beispiel (33) 


Sis 
Cec 


Ist eine geschlossene Lésung nicht mehr méglich, so muB die Determinante 4 fiir verschiedene 
Frequenzen q@ (in der Umgebung der zu erwartenden Eigenfrequenzen) berechnet werden. Die 
Nulldurchgange der Kurve A= A(w) oder A = A(A) ergeben dann die gesuchten Kigenfrequenzen.” 


A= =0. (34) 


=Sc—Cs=0, oder mit (4) te MN SGA (35) 


6. Determinanten-Matrix. In der Frequenzgleichung (34) tritt eine Determinante zweiten 
Grades aus Elementen der Matrix 2{ als Eigenwert-Determinante auf. Nun besteht ein Balken 
oder Balkenabschnitt im allgemeinen aus mehreren (in sich homogenen) Feldern. Die Feldmatrizen 
sind nach (29) auszumultiplizieren, und aus vier Elementen der Produktmatrix soll eine Determi- 
nante zweiten Grades gebildet werden. Fiihren wir diese Aufgabe fiir zwei Matrizen in allgemeinen 
Zeichen durch, d. h. bilden wir die Determinante 


4=\ 6 Gs io 

aus der Produkt-Matrix © = %- Y, so wird 
A = S Bai Aic Bor Aka — Bsi Aic Bak Aka - (36b) 
Von den 32 Gliedern hie Oe Npeleee heben sich die 8 Glieder mit i = k weg. Unter der 
Summe kénnen wir 4;, Aza heraussetzen, der Rest 1aBt sich zu der Determinante aa ae = 


sammenfassen. Unterscheiden wir nun bei der Summation zwischen i < k und i> k, so kénnen 
wir fiir (36 b) weiter schreiben 
4 3 


Bui Bae 


3 4 B A B k 
at ai Da ; Ane Amas 36 
A ae By; Bor ee ene By; Box | a 89) 
Pea = 


1 §. Timoshenko, a. a. O., S. 345. ree 
2 | Holzer’sches Verfahren“ in der Theorie der Torsionsschwingungen. Die Ubertragung dieser Gedanken 


auf das Biegeproblem mit Hilfe der Matrizenschreibweise hat neuerdings auch E. Pestel gegeben. (Abh. 
d. Braunschw. Wiss. Ges., 1954, S. 227.) 
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vertauschen wir im zweiten Glied i und k, sowie die beiden Spalten der Determinante, so konnen 
wir beide Glieder wieder zusammenfassen und erhalten! 


A= s by Binoy (Aic Aka — Ale Aia) » (36d) 
Vi h— 2 Boi Bor 
ick 
oder mit (36 a) ; 
Cie Caza) ae, Bai Bak " Ape Asa 37 
Pe Cha > Saceal Bee Biel ide aete ae oo 
ah 


Um mit den hier auftretenden Determinanten zweiten Grades bequem und ibersichtlich rech- 
nen zu kénnen, miissen wir sie durch ein kurzes Symbol kennzeichnen. Wir brauchen hierzu auBer 
dem Zeichen fiir die Matrix der die Determinante entnommen ist, vier Indices, zwei fiir die Zeilen 
und zwei fiir die Spalten. So kénnten wir fiir die auf der linken Seite der Gleichung (37) stehende 
Determinante etwa Ge schreiben. Wir kénnen mit weniger Indices auskommen, wenn wir die 


b 
Reihenfolge der Zeilen- und Spaltenkombinationen durch ein Schema festlegen, etwa das folgende: 


Zéilen oder Spalteaia’ (12, dau daleesml moana 
og) (38) 


Kombinationsindex : 1 2 a 4 5 6 


Hiernach bezeichnet das Symbol C4, diejenige Determinante, deren Elemente nach dem Schema 
(38) durch die p-te Kombination zweier Zeilen und die g-te Kombination zweier Spalten gefunden 
werden. So ist z. B.: 


C C 
A 12 14 
Cr = | Op Sy ° (38a) 


Wir wollen uns im folgenden dieser Schreibweise bedienen, nicht nur weil sie kiirzer ist, sondern 
weil wir dann den Matrizenkalkiil auf diese Determinanten tibertragen kénnen. Es sei aber als Nach- 
teil vermerkt, daB zwischen den Kombinations-Zeichen p, q und den Indices der Determinanten- 
Elemente a, b, c, d kein formal zwingender Zusammenhang besteht. Der Zusammenhang ist viel- 
mehr durch das bis zu einem gewissen Grad willkiirliche lexikalische Schema (38) gegeben. 

Mit der Verabredung (38) kénnen wir nun jede der in (37) auftretenden Determinanten durch 
ein Symbol entsprechend (38a) darstellen. AuSerdem wird dann aus der in (37) stehenden Doppel- 
summe iiber t, k eine einfache Summe; denn in (37) wird fiir die Kombinationen i, k in den zwei 
Zeilen a und b der Matrix B die gleiche Reihenfolge vorgeschrieben wie in den zwei Spalten c und d 
der Matrix 9{. So kénnen wir fiir (37) schreiben 


(AY eek . A A 
C4, = 2. Baa (39) 


Diese Formel, die genauso gebaut ist wie die Produktdefinition der Matrizen, legt es nahe aus 
den nach (38) definierten Determinanten C3, eine neue Matrix mit 6 mal 6 Elementen zu bilden, 
die wir als die erste Determinanten-Matrix der Matrix © oder kurz als 4-Matrix &4 bezeichnen. Aus 
(39) folgt dann die Matrizengleichung 


C4 = 84-4 oder = (BW A)4A = B4- 4, (40) 


in Worten: Die A-Matrix eines Matrizen-Produktes ist gleich dem Produkt der A-Matrizen.2 

Bei der Determinanten-Auswertung (36 b) zeigte sich, daB ein Viertel der auftretenden Glieder 
sich gegenseitig wegheben. Nun kénnen gerade diese Glieder sehr groB werden gegen die iibrig- 
bleibenden. Das wiirde zur Folge haben, daB bei der Zahlenrechnung das Resultat als Differenz 
von zwei nahezu gleichen Zahlen sehr unsicher wird, und da® ein Grenziibergang zu dem unbe- 
stimmten Ergebnis 00 — oo fiihrt, wenn man nicht die betreffende GréBe, mit der man zur Grenze 


ibergehen will (z. B. die Steifigkeit ¢ eines elastischen Lagers), als Parameter bis zum Schluf8 mit- 
nimmt. 


1 Vel. W. Schnell, a. a. O. 
2 In den technischen Anweisungen ist diese Matrizenformel bisher wohl noch nicht bekannt gewesen. 
In der Algebra wird sie im Zusammenhang mit dem Laplace’schen Entwicklungsgesetz einer Determinante 


ae (sae aan bewiesen; z. B. H. Hasse, Aufgabensammlung zur Hoh. Algebra, Goschen 1082, 
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Der durch die Beziehung (37) oder (40) gegebene Zusammenhang gestattet uns nun, den letzten 
Schritt zur Ermittlung der Eigenwert-Determinanten 4: den Ubergang zur Determinanten, be- 
liebig weit vorzuverlegen, gegebenenfalls bis zu den einzelnen Feldmatrizen. Damit kann nicht 
nur ein Teil des Rechenaufwands ge- 
spart werden, sondern es 1a8t sich vor 


\ 
allem auch die numerische Unsicherheit — aa aes — 


des Ergebnisses beseitigen und ein 


_ Grenziibergang (z.B.c—> 00) schonin =z A ; 5 . , ; 

den Feldmatrizen durchfiihren. ] Cn : [Me 3 Cry 5 |& Og 
; Die A-Matrix & enthalt simtliche — (ea ie cA cA od oe 
méglichen Determinanten zweiten Gra- = Sememe) Ee Es 
des der Matrix ©. Ihre Elemente C+, gfe Csi Ci C33 ee he ies 
sind also die Kigenwert-Determinanten : cA ci ca : a4 = if i; y 
eines durch © erfaBten, aus beliebig il ~ bee Sal bo SEA a iad 
vielen Feldern bestehenden Balkens bei L— Cs4 We) bs Ge Os ips 
allen méglichen Kombinationen der a oA cA 7A a 4 .: 7 
rechten und linken Randbedingungen. be vi HY 64 65 | 6 
Fir jede Spalte gelten zwei bestimmte 
Abb. 7. 


rechte, fiir jede Zeile zwei bestimmte 
linke Randbedingungen (Abb. 7). 

Die dritte und vierte Spalte oder Zeile wiirde fiir w = 0, Q = 0 und y = 0, M = 0 gelten, was 
mechanisch nicht denkbar ist. Die Frequenzgleichung fiir den links eingespannten, rechts gelenkig 
gelagerten Balken hei®Bt hiernach z. B. C}, = 0. 


7. A-Matrix der Rayleigh-Matrix. Wir wollen aus der Rayleigh-Matrix { (8) die zugeordnete 
A-Matrix Si4 entsprechend der Verabredung (38) gewinnen. Dabei treten durch Determinanten- 
bildung Kombinationen der Funktionen (4) auf, die sich durch die vier neuen Funktionen 


: 1 ej 4 eel 2 4 26 78 
B = 575 (Gof sin A — @in / cos A) = el as Thi 

SA ee (ieee 24 2° are 
S* = 553 Sin A sin = o> 6% + i014 Tee 

1 ; ; 2? 24 2° aie 7 
A* = 57 (CofA sin 2 + ©111/ cos A) =1— 517 = aT 3t4 

2 28 2° 

E* =1+ (1 + Go{/ 00s ) Se a ia 


darstellen lassen!. Auch diese Funktionen kénnen in gut konvergierende Reihen entwickelt werden. 
Die der Feld-Matrix $j zugeordnete A-Matrix lautet damit 


[ 2 B2 Bes, sah amet Ee | 
ee ee 
3 
— 15 B 2 E*—] p A* pat | 2p St Bh 
P 
— 135, 5" —4 Be ee ae Eas 
4 = | tig ee == = i (42) 
AF S* ei B* rea) E* BoA® = S* 
2 
o at ee Ag Ay * pea | B A* 
moe ee a a 
2 o O o a 
ApoE ape —Aps ARs Ape] || 


1 Diese Funktionen sind im wesentlichen die bei W. Prager und K. Hohenemser, a. a. O., S. 323 tabulierten 
Funktionen. 
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M4 ist wie Ht symmetrisch zur Nebendiagonalen. Entsprechend Abb. 7 sind die Elemente Ry; 
die Eigenwert-Determinanten fiir den homogenen Balken bei allen méglichen Kombinationen 
der rechten und linken Rand- 


bedingungen'. Die Lésung 


g i ZY +. *. . . . 

;— Ah [/__— B*-0, fir den beidseitig einge- 
a ; spannten Balken heiBt hier- 
| TgA=-tgd cea TgA=tgA nach z.B. 

= h* — 0 oder 


1=E "30; ae S*=0; CofAcosA = 1. 
" "\ 


y y : Ferner ist zu erkennen, 
j 7 pee ce nS DX Ue daB die in Abb. 8 zusam- 
mengestellten vier Paare je 
die gleichen Kigenfrequenzen 
haben. 

Wie fiir die Matrix kénnen wir auch fiir die in Ziff. 3 gefundene Matrix (25), bei der Rotations- 
tragheit und Schubelastizitaét mit beriicksichtigt sind, die zugeordnete A-Matrix finden. Die hierzu 
aus (24) durch Determinantenbilden abzuleitenden Funktionen sind jedoch naturgemaB kompli- 
zierter? als die entsprechenden Funktionen (41) fiir die Matrix 4. Es seien hier nur die beiden 
Sonderfalle €4 und INt4 angefiihrt, die sich auch aus (& (27) und 9)t (28) unmittelbar ergeben: 


| 1 Bile ee B20 PP la(s +6) Bila (z's + 9) 


Abb. 8. 


; heal i . aS, s 
woe [20 
Pos 10 One r Ra |? (43) 
10 0 0 0 1 Bla 
eer ye 0 0 1 Ue 
r ine 0 0nn0 0 01 
One 1 0 0 0 0 
{Bul 
0 0 The 2H 0 0 
a 0 0 ea 0 0 
Me : (44) 
— m; wo? —= 0 0 O it 0 
EI 
m,0;o1l' B 1 
io —_— i a 0 0 ey — 
(EI)? eT SS gi 


8. Elastische Stiitze. Eine elastische Stiitze kann durch eine Stiitzen-Matrix © erfaBt werden, 
die zwischen den Vektoren {); und ));,, vor und hinter der Stiitze vermittelt: 


Dit1 = O° ye - er) 

Die Stiitze habe eine Quer-Federkonstante cg und eine 

GE Dreh-Federkonstante cy. Das zwischen den Stellen i und 
i + 1 liegende Balkenfeldi habe auBerdem die Einzelmasse m 

z4 mit der Drehmasse Q, jedoch keine Lange: ]->0 (Abb. 9). Es 


ist dann 


Stelle: +7 


Wef1 — Us Vee h ie ee (45b) 
und es gelten die Bewegungsgleichungen 
mw;=Q:4,;—Qi— cow, O yi = M;— Mi4,—em yi, (45¢) 


oder fiir harmonische Schwingungen 


Abb. 9, 


— asl — (cm —- (0) w?) Wi M; = Qins = (m 2 — co) WwW; — Q; . (45d) 
Schreiben wir zur Abkiirzung 
1 —B 
a OM 9) =I» =, o—mat) —Io, (46) 


1 Ri, = Ry, usw. sind mechanisch die Folge des Maxwellschen Vertauschungssatzes, siehe Diss. W. Schnell. 
2 Siehe Diss. des Verfassers. 
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so lassen sich die Aussagen (45b) und (45d) in Form der Vektorgleichung (45a) schreiben mit 


eae le Omg 104 
| 2 aren 
Se ho yt 6 (47) 
Die A-Matrix der Matrix © wird dann 
an Qo bs 0°= 204 
i ot Ott 020 
0 Rf ene an 
> ae 
oo at ost One 6-0: | (48) 
emi On ie 20 oT <0 
Tie oe 0rd Al 


Mit cy —> 0 und cg -> 0 wird aus (47) die Massenmatrix Jt (28) aus (48) Wt4 (44). 

Ist ein Balkenfeld i von der Linge 1; elastisch gebettet mit der iiber I; konstanten Bet- 
tungsziffer k;, so daf} die drtliche Federkraft dK = — wk dx ist, so gilt die Bewegungsgleichung 
pw = Q’ —kw oder fiir harmonische Schwingungen (k —w? uz) w = Q’. Damit gilt auch hier die 
Differentialgleichung (1) und die Matrix (8), jedoch mit 

wo u—k 
Jn —ay  t- (49) 
Mit Hilfe der Matrizen © oder G4 1af®t sich ein Balken mit beliebig vielen elastischen Lagern in 
sehr einfacher und iibersichtlicher Weise erfassen. Wir wollen das am Beispiel eines Balkens mit 
drei elastischen Stiitzen zeigen. 

Die Abschnittsmatrizen 2{ und 8 fiir die Abschnitte 2 und 4 zwischen den Stiitzen seien aus 
beliebig vielen Feldmatrizen durch Matrizen-Multiplikation (29) gewonnen. Dann gilt wieder nach 
(29) fiir den ganzen Balken mit den Randbedingungen Abb. 6, Fall 4, wobei die Stiitzen mit 1, 3, 5 


bezeichnet seien, 


a rom Fay | 
L9 J6 L941 LO 


wobei © durch Ausmultiplizieren der fiinf Matrizen gewonnen wird. Die Eigenwertgleichung heift 
dann entsprechend (34) 
Cy, Cp 


6 
41 42 
Fiihren wir diese Rechnung in allgemeinen Zeichen durch, wobei wir J'y = 0 annehmen und fiir 
Ig kurz I" schreiben wollen, so wird unter Verwendung der Verabredung (38) 


A= =O, (51) 


3 An Bs; : S Ais By: 
KSA By Ty 
— Ay, Bay rs i— Ale By * D's | 
ee = 544, Bs 
oe pean +2 At, Bi te 
— 2» Ai, By I; | =e (Ay, Bis + Ax Be.) : Ds 
ae a By, Ls = Ajo Bay 1 zs 2 ze eg By a i > (52) 
Se es aye + (Aj, Bes + Ad, BS) Ls 
feeds Bae I, L's | + 2 At, By, ids 
203) sed Ee en + (48 Be + 48 Bh) Tals 
+A Tals +A Bu ToT) 4 (ag Bh + ag, Bh) Ly DT 
— Aj, By i3l's | ae 0 hd be 
aes TH, 
+ 0 cram 
= 
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Hier zeigt sich deutlich das in Ziff. 5 Gesagte: Bei der Zusammenfassung verschwinden 4 von den 
12 Gliedern, darunter gerade das Glied mit der héchsten T,-Potenz. Dies hat die sehr unangenehme 


Folge, daB, wenn wir mit Zahlen rechnen, d. h. nicht mehr nach [-Potenzen ordnen kénnen, die | 


Rechnung insbesondere bei groBen Werten von J; sehr unsicher wird, und da®B der wichtige Fall 
starrer Lager, also J’; > 00 mit diesem Verfahren nicht erfaRt werden kann, wenn wir nicht [}; = J” 
als Parameter bis zum Schlu8 mitnehmen. 


Die Lésung fiir den Sonderfall starrer Lager liefert der Faktor mit der héchsten ver bleiben- | 


den /-Potenz. Hier also 


Ail, BY, + Ay, By = 0. (53) | 


Wenn wir an Stelle der in (50) verwendeten Matrizen gleich die zugeordneten A-Matrizen, 


entsprechend (38) beniitzen, so kénnen wir nach (40) 


Gs -B4. S4- 94. Gt = 64 (54) | 


bilden, und die Eigenwertgleichung ist dann entsprechend Abb. 7 


C4 =0. (55) | 


Fiir die Rechnung brauchen hier die Matrizen nicht voll besetzt zu sein. Wir nehmen nur die not- | 
wendigen Elemente mit (bei der Gewinnung von 9(4 aus 9 und B4 aus 8 kénnen wir uns ebenfalls | 


darauf beschranken) und erhalten fiir (54) ausfiihrlich 


ee io ae een inCreate ar fo Ee SIE noo - 
oe woee | HBA cBS | Ba. BA, Bax BA |e rian eee ae 

cok Sars lle ice ace gee ovo konto AA ee ee ee a 
DS Saw bese 2 epee Ae a Bln eg tLe We A eee 6 
oe ae Pl gesee Lee: Py. ..d.| AdTe tA re Sere | 

Bel Ged |) BAL BS BasBAt BAG BA ae aii. a4 aoe oe | (56) 


Fiir (56) kommt durch Ausmultiplizieren der Matrizen unmittelbar 
ea Si (BA Ad BAA). Pe 
+2 BY AL-T, A (BeaaS Be age le . 
+3 Bi 44D, + (Bf, 44, + BS Ad) TT "4 
+ & BY AST, T, + (BA 44, + BY Ad) PP, T, =0. 


Im Gegensatz zu (52) treten jetzt die hdheren /-Potenzen, deren Koeffizienten dort verschwanden, 
gar nicht mehr auf. Der Balken mit starren Stiitzen kann jetzt sehr bequem als Sonderfall erfaBt 


werden: Eine starre Stiitze an der Stelle i bedeutet I; > 00; in der Frequenzgleichung sind dann 


alle [’;-freien Glieder klein gegen die Glieder mit J’; und kénnen unterdriickt werden, wahrend die 
Summe der mit /’; behafteten Glieder (J; selbst ausgeklammert) fiir sich Null werden mu8, wenn w 


ein Eigenwert sein soll. Das bedeutet, daB die A-Matrix ©4 (49) degeneriert zu der A-Matrix fiir | 


die starre Stiitze 


TO WO 
PO ad OOW 
0700 OmOENO 
Ce | 
ee 0s. 0220s Om OREO Ns (58) 
1008000. 20 
0. Fl Om O uO ec0s 


Nehmen wir etwa im Beispiel oben in der Mitte ein starres Lager an, wahrend die Au®enlager 


elastisch bleiben mégen, so ist in (56) die mittlere Matrix durch ©? zu ersetzen. Von (57) bleibt 


dann nur der rechte Teil mit i= Dubrig 
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Sind alle drei Stiitzen starr, so wird aus (56) 


* ee ee ee Oe m7 qpPep eeistol mygaiel ae r 
A Ra / 
en eee eB IMA a, 2 | be <A, | 


Cie etnies ela ta 


ieee es eer Oe ae velo slliem ce op sy uel) ay iiew lie e (hye © 6.’ Il! esceve” Sete ve 


ee IA AR HAE Mom tit) be ee. 
und A = Cj, = By, Ay + BY, Ad, =0, wie nach (57) mit I, =I, =I, > o. 


9. Eigenschwingungsformen. Das homogene, lineare Gleichungssystem, das zur Eigenwert- 
gleichung fiihrt, z. B. (33) hat 1- 00 viele Lésungen, wenn w eine EKigenfrequenz w, ist. Man kann 
daher von n Unbekannten dieses Systems eine frei wahlen (z. B. = 1 setzen) und hat dann noch 
n — | Gleichungen fiir die restlichen Unbekannten. Auf diese Weise lassen sich die Komponenten 
des Vektors ),; am rechten Rand des Balkens finden. Im Falle (33) z. B. folgt aus w, = 0; M, = 0 


und y, = 1, also M, = 0 und y, =1: 
— Any Aw, (60) 


Aus der Matrizengleichung (30) ergeben sich die Komponenten des Vektors {);, insbesondere w;, 
an beliebigen Stellen i, wenn Y =(w,) die Matrix fiir das zwischen den Stellen 1 und i liegende 
Balkenstiick ist. So erhalten wir z. B. fiir den links eingespannten, rechts gelenkig gelagerten homo- 
genen Balken die Eigenfunktion aus (8) mit (60): 


1 A,,(l) 7 x x\ CA 
w(x) = Aj(x) 1 — Ayy(x) 0) Bee s(r a Sa ( a P ss 


sin 4 + CofA — Sin A= cos A 


l : 
= = ; il 
7 ie ee » mit tgA=Zgd_ nach (35) (61) 
Fiir den Balken mit elastischen Zwischen- —Adsohp/7: 3 Qt 
stiitzen gilt dasselbe, wenn man unter Qf jetzt Feld: 4 (t-1) + -b (b-1) a 


das (die Stiitzenmatrix enthaltende) Produkt 


aller Matrizen zwischen 1 und i versteht.  #lle:c sel ae oe 
Aber auch bei starren Zwischenstiitzen labe (links) (rechts) 
sich die Eigenfunktion ohne Schwierigkeit ASE. 10. 


bestimmen: Aus der Bedingung w, = 0 1 
(Abb. 10) kann die fehlende Komponente von {), (hier Q.) gewonnen werden. Es ist fiir py. = 1 


= = = = ATS 
‘tp, = Vs v1. n= Us a ai ie 
Aus ty, =%U-Yafolgt * 
Wp =—V, 
— 7 Aint. EG 
yo = Ano l— Ang Ge pe ae (62a) 
Ay, 7 As. 
M, Ae ag! = ae 


Os; wird hier nicht gebraucht. Oo1 die Querkraft links von der Stiitze b, folgt aus der Bedingung 
We = Bip Po + Bys Mo + By Qo1 = 0 zu 


On = Are By + Aus By; ie (62b) 
Ay, By, 


Damit ist ) an jeder Zwischenstelle i bekannt, so ist z. B. 
w; = Vip Yo + Viz M, Si Va Qui» (63) 


wenn die Matrix § =G,_, Bi-2 - - - Bs ist. 
23 
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Wenn w keine Eigenfrequenz ist, werden die Randbedingungen und Ubergangshedingungen 
nicht simtlich erfiillt sein. So ist im letzten Beispiel die Bedingung 


— = — = AA BS ee eee 
M. = Byg Po + Bs Ms + Bog Qs = = a ime “~i=0 (64) 


noch offen. Sie ist nur fiir ~ = «, erfiillt und liefert daher die Kigenwertgleichung 
A(w) = Af, By, + Az, By, = 9, (65) 


die, wie es sein mu}, mit (53) ibereinstimmt. 


M, ist dasjenige Erregermoment, das eine Schwingung mit der Anfangsneigung ig = Aleta 
zwingt. Es ist Null fiir =«,. Statt der Nulldurchgange der Frequenzfunktion A(w) kann man 


also auch die Nulldurchgange von M,(w) zur Bestimmung von W = We bentitzen. 


10. Erzwungene Schwingungen. Wie Ziff. 9 zeigt, kénnen wir, im Rahmen der hier behandelten 
Matrizenmethode, erzwungene Schwingungen dann auf sehr einfache Weise erfassen, wenn die 
erzwingende Ursache (Kraft oder Verriickung) an einem Rand auftritt. Die (i — 1)-Felder eines 
Kragtragers z. B. seien zu einem Abschnitt mit der Matrix 2 zusammengefaBt. Am freien Ende 1 
greife eine mit der Frequenz 2 pulsierende Kraft an, deren Amplitude P; ist. Die Randbedingungen 
lauten dann 


M;=0, ene Pie Spe 
I (66) 
oe ies 
und wir erhalten, wenn wir die letzten zwei Gleichungen von tj; = W- YH, verwenden 
Va ast iy = Toh Ag3 
eS me Q; Pr (67) 
Ferner ist z. B. 
Val A PD Ass 
iA Miata Ota ete aa (68) 
66 
Fiir den masselosen Kragtrager mit Endmasse wird 
h 1 
ie Es iis ee te (em 


wobei h = 13/3 EI die EinfluBzahl des Kragtragers, w. = P;h die statische Durchsenkung 
ind 97) On, tnt 1//mh die Abstimmung ist!. 

Bei einem symmetrischen Tragwerk kénnen wir den Schwingungsansatz auf eine Symmetrie- 
a halfte beschranken. Die Symmetriestellen werden 
v dann Randstellen, so daB wir nach der eben be- 
schriebenen Art auch Erregerkrafte erfassen kén- 

Wt, | Q.1 | | Q.., | GY nen, die an einer scan ete angreifen. 
Anders liegen die Verhaltnisse, wenn die er- 


3g regende Kraft (oder Verschiebung) nicht an einem 
y Rand wirkt, sondern innerhalb eines Abschnitts. 
eR Greift an der Stelle i die Erregerkraft P; an, 


so gilt fiir die Querkrafte unmittelbar links 
bzw. rechts von i die Gleichgewichtsbedingung (Abb. 11) 


Qi1 == Qi, = P; > (70a) 


wahrend w;; = w;,, pit = pir und Mi, = M., bleibt. Wir kénnen also schreiben 
my 
— |0 
Yit =Yir — Prey) = (70b) 
| 


" K. Klotter, Techn. Schwingungslehre Bd. 1, S. 65 u. 167, 2. Aufl. Berlin, Gottingen, Heidelberg 1951. 
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Die Stelle i mége zwischen beliebig viclen, beliebig gearteten Feldern 1 bis (k — 1) legen (Abb. 12). 
Die Feldmatrizen aller Felder rechts bzw. links von i seien zusammengefaBt zu den Abschnitts- 
matrizen 2{ bzw. B. Dann gilt 


A 
Ye =B- Yi, Yir =W-H,. (7c) <i ine ines 


df L a 7 
Die Verbindung zwischen diesen Glei- x ote iy o = D 
_ chungen wird durch (70b) hergestellt, B 94- 
und wir erhalten durch Einsetzen Kiwi: 
lO! 
= bp a 
Yr =O 1 A, — P, o| J=BAH, — P: by, (70d) 
IL 1 
' wobei b, der vierte Spaltenvektor der Matrix % ist. Setzen wir 
S- x=, (70e) 
so kénnen wir fiir (70d) schreiben be 
©: ), =t. + P.b,. (71) 


Fiuhren wir nun die Randbedingungen fiir die Rander 1 und k ein, so wird aus (71) z. B. fiir beid- 
seitig freie Enden 


em | fF Rae 
| 21 | Ga a 
©.) Vi] =| Pel p,. | Bost (72a) 
| 0) | o i Be 
e0e. 194 | Bas) 


Aus diesem Gleichungssystem greifen wir die zwei Gleichungen heraus, die das Verschwinden 
zweier Komponenten von {), ausdricken 


Ae 


= = (72b) 
Cy w, + Cy, Y, = P; By. 
Hieraus erhalten wir = 7 
w,=— = Cs2 Boy 4 Y, = C5, By, (73) 
Cox Cy. By Cox Cy By, 


Fiir die hierbei auftretenden Zahlerdeterminanten, bei denen die erste Spalte aus Elementen der 
Matrix © =%- 2 besteht, kénnen wir allgemein die folgende Entwicklungsregel aufstellen: 


Gre Bra = Ba; Bea 
Coe Boa By; Boa 


Entwickeln wir namlich die links stehende Determinante nach Elementen von Y{ und %, so kénnen 
wir die entstehenden 8 Glieder zu der rechtsstehenden Summe wieder zusammenfassen, Dabei 
tilgen sich 2 der 8 Glieder (von den rechts auftretenden Determinanten pie Se Ouee mit t = d zu 
Null). Die restlichen 3 Determinanten sind nach Ziff. 6 Elemente der Matrix B4. Die Beziehung 
zwischen den Indices mn von B4 _, und den Indices a, b, i, d der zugehérigen Determinanten-Ele- 
mente laBt sich nicht unmittelbar formulieren; sie geht aus dem Schema (38) hervor. 


Mit (74) erhalten wir fiir )), aus (73) 


Ase (74) 


es) 
j=1 


Pi 
Wises eee (Be. Flee =f (ayes As Ba Be A;,) 9 
w Pi 4 1 { HS 
, = Yr ra C4, (Be, Ay, af Be As sis dob As) ? ( ) 
M, =i) > 


0,=0. 
An einer beliebigen Stelle d zwischen den Stellen 1 undi erhalten wir dann die Komponenten des 
Vektors jg mit Hilfe der zwischen ] und d vermittelnden Matrix ® (Abb. 12) aus 


ha =D- yp. (76a) 
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Setzen wir die Beziehungen (75) ein und fassen die Elemente aus 2{ und D zu Determinanten 
zusammen, so wird 

P; 
Sra 

C1 


Dy Dy» 


Dy Dy» 5 vi 
31 32 


11 12 


+ Bé [Du Mie) + Be ) (6b) 


Sie == 


(Be 


und entsprechend wa, M, und 03. 


Auch hier kénnen wir fiir die jeweils im Zahler auftretenden Determinanten, bei denen die 
zweite Zeile aus Elementen der Matrix 2( = %- D besteht, eine allgemeine Entwicklungsregel auf- 


stellen: 


21 22 


4 


Dae Deg 
= 12 Dia 


> Fai UT, 

Aye Ava : a 
Wie bei der Regel (74) fallen auch hier von 8 Gliedern 2 weg, von den rechts stehenden Determi- 
nanten wird diejenige mit i = a zu Null. Die restlichen 3 Determinanten sind Elemente der Matrix 


M4 (Ziff. 6). Hiermit kénnen wir (76b) auf die folgende Form bringen: 
_ | Bi, (Di Fis + Da Fis + D5 Fi) 
Fi) + B4 (D4, Fe, + Dé, Foy + D5, Fos) 
Cin | + BS (Df Fs. + Di Fs + Di, Psa) 


Pepsi 


| Dire Dy 


(78) 


Wd = 


aI 


und entsprechend wa, Ma und Qu. 
Fiir D = Y und § = & (Einheitsmatrix) erhalten wir hieraus den Vektor \);, unmittelbar rechts 
von der Erregerstelle i 


wi = —Fa (Bi, Ai, + Bd Ad), 
v;,= Pi A 4A A ga 
ees C4, (Bes Agi — Ds A‘) ? 
i) (79) 
p= 4 (Bes Any + By; Aa) » 
0; = Pi A 4A A A 4A 
es C4, (2, tke =f: Bé, Axi ae Bee Az) * 


Der Vektor {);, unmittelbar links von der Erregerstelle i, unterscheidet sich von Yi; nur durch den 
Querkraftsprung, der in den Gleichungen (70a) und (70b) zum Ausdruck kommt. Fiir eine Stelle j 
zwischen i und k erhalten wir, wenn die Matrix <) den Abschnitt zwischen i und j erfaBt (Abb. 12), 
mit (70b) 


IV Gr 


7 = 3° da = 3 on — PB (80) 


Setzen wir {);, aus (79) in diese Gleichung ein, so erhalten wir 

; i Ans (Bi, Is, eS Be. Ts ah BY, Is) 

(¥s); = maa : a Ax (Ba Ts ae Be Is3 Ts By, Is4) > (81) 
| + Ad, (Bg, Len — Be, Teg + BE, Is,) | 


66 


L a 


wobei sich fir s = 1, 2, 3, 4 die Komponenten w;, Wp M,, Q; ergeben. 

Die Gleichung (81) hat [wie die Gleichungen (73) und (76b)] eine sehr unangenehme Eigen- 
schaft : Im Zahler treten Produkte yon Elementen der Matrizen 94, B4 und <j auf, wobei die 
Matrix % die Matrix J enthalt (Abb. 12). Liegt nun z. B. im Abschnitt der Matrix $ ein quer- 
elastisches Lager, das wir durch eine Matrix © (47) mit cy = 0 erfassen, so kommt in den Zahler- 
produkten die Federzahl c dieses Lagers zweimal vor: innerhalb der Matrizen J und B4. Die so 
entstehenden c?-Glieder im Zahler heben sich wieder weg, wenn wir die eckige Klammer von (81) 
ausrechnen. Abgesehen vom unnétigen Rechenaufwand wird dadurch die Rechnung sehr unsicher 
wenn ¢ grof} ist und fiir starre Lager mit c—> oo miiBten wir c als Parameter bis zum Schlu8 Henen 
lassen [vgl. auch (53)]. Um diesem Ubel zu begegnen, wollen wir (81) auf eine Form bringen, bei 
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der eine Uberlappung von Matrizen nicht mehr eintritt. Dies ist wieder durch eine Entwicklungs- 
regel ahnlich den Regeln (74) und (77) méglich. 


Zunachst wollen wir die drei in runden Klammern stehenden Ausdriicke von (81) als Determi- 
nanten dritter Ordnung schreiben: 


P; thes, Ik Tea | | Zsa Iles They Tks Ihe Tq | 
(¥s)j a C4 aa Bs, B35 Poa 7 An Bsy Bs Bz, ells Ale, By» Bs, Bs, | . (82) 
ee a1 Bay By, | | a1 Bag By, Byy Bag Be 


Fur die hier auftretenden Determinanten dritter Ordnung gilt allgemein die folgende Entwick- 
lungsregel: Wenn B = St: § ist, so gilt 


eat h et  e e [ead Be 

sa s Sc vy | 4 , |~ sa sb Sc 

Bi Ba Bye| = = ihre ae Via Tin Tiel; (83) 
UB ie Bar Bic ick tie ee Wie Tis ig 


Die Indices i, k durchlaufen dabei die sechs Kombinationen des Schemas (38). Die unter der 
Doppelsumme stehenden Determinanten zweiter Ordnung sind Elemente der Matrix §. Von den 
sechs Determinanten dritter Ordnung aus Elementen der Matrix <} werden stets 3 zu Null, und zwar 
diejenigen mit i = s oder k =s. Die drei nicht verschwindenden sind Unterdeterminanten der 
Matrix o: oder mit dem Vorzeichen (—1)P*4 versehene Elemente der Kehrmatrix Sie = [eels da 


die Gesamtdeterminante J — |S | = 1 ist. Somit erhalten wir aus (82) mit Hilfe der Entwick- 
lungsregel (83) 


P; | ‘AG (Sin Ls Fg Kgs in — TS De) | 
mj = — ca cw eels dl pyre 1G Praca Cau Ee) (84) 
ea = BG ii We KG. Ls ~— 1Gy, er) ei 


und entsprechend ¥,, M;, Q;. 
Mit Hilfe der Gleichungen (78) und (84) kann nun der Vektor i) fiir jede beliebige Stelle des 


Stabwerks bestimmt werden, und zwar auch dann, wenn das System beliebig viele elastische oder 
starre Zwischenlager enthalt. Bei starren Lagern kann dabei der Grenziibergang c—> oo bereits in 
den Feldmatrizen vorgenommen werden. 


11. Statische Einflufzahlen. Bei komplizierten Rahmengebilden wird man oft genétigt sein, 
das wirkliche System durch ein angenahertes mit n diskreten Massen” zu ersetzen und steht dann 
vor der Aufgabe, die statischen EinfluBzahlen fiir die Angriffsorte dieser Massen zu gewinnen. Dies 
kann natiirlich nach einem der iiblichen statischen Verfahren geschehen: Die EinfluBzahl hj; ist 
die Verschiebung (Drehung) an der Stelle a durch die Last 1 (das Moment 1) an der Stelle 6. Wir 
kénnen fiir diese Aufgabe aber auch das Matrizenverfahren nutzbar machen. 

Im statischen Fall geht die Matrix fiir ein homogenes Balkenstiick mit 4 — 0 wher in & (27), 
wenn wir die Schubelastizitat mitberiicksichtigen. 

Auch hier gilt wieder die Gleichung (29) fiir einen aus mehreren Feldern bestehenden Abschnitt. 
Wir kénnen nach Ziff. 6 die Determinanten-Matrix einer beliebigen Feld- oder Abschnittsmatrix 
aufstellen und erhalten so z. B. aus © die 4-Matrix (©4 (43). Ferner gilt die Matrix (47) fiir elastische 
Lager mit m = 0 bzw. O = 0, es gelten auch die im II. Teil noch zu behandelnden Matrizen fir 
Rahmentragwerke. 

Auf das so durch Matrizen dargestellte elastische System soll nun an einer Stelle die Last 1 
(oder das Moment 1) wirken. Wir wollen uns hier auf den Fall einer duBeren Querlast P; = 1 
beschranken. 

Wie bei den erzwungenen Schwingungen (Ziff. 10) 1aBt sich eine AuBere Kraft P; dann sehr ein- 
fach erfassen, wenn i eine Randstelle ist. Symmetriestellen eines symmetrischen Tragwerks sind 
ebenfalls Randstellen in diesem Sinne. 

Greift die statische Kraft P; = 1 nicht an einem Rand, sondern innerhalb eines Tragwerks an, 
so kénnen wir wieder genau so verfahren wie im Fall der erzwungenen Schwingungen. Mit P; = 1 
gelten alle Gleichungen von Ziff.10 auch hier, wobei die Matrizen D, F und G, & jetzt reine Elasti- 
zitatsglieder reprasentieren. Insbesondere erhalten wir wieder 1), aus (78) und 4); aus (84), also 


1 KR. Zurmihl, Matrizen, S. 36. Berlin, Gottingen, Heidelberg 1953, 
2 H, Fuhrke, Stahlbau 23 (1954), S. 181. 
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jede beliebige der vier mechanischen GréBen w, py, M und Q an jeder beliebigen Stelle des Stab- 
werks. Fiir die Stelle i kénnen wir die Gleichungen (79) verwenden; so wird z. B. 
1 Bis At + Boo An 

Hore Ce 


hi = (85) 

Hat man auf diese Weise das Schema der Einflu8zahlen h;; bestimmt, so 1aBt sich die Eigen- 
wertgleichung — in anderem Sinne als bisher — als Matrizengleichung formulieren. Bezeichnet 
* die Matrix der EinfluBzahlen, Jt* die Massen-, ©* die Einheitsmatrix und A = 1/@? den Eigenwert, 
so ist! 

H* M* —A C* = 0 (86) 
die Eigenwertgleichung, die sich z. B. mit Hilfe des Verfahrens von Hessenberg? in eine Polynom- 
Gleichung n-ten Grades iiberfiihren laBt, deren n Wurzeln dann die n Eigenwerte des Systems 
sind. Da die Matrix )* symmetrisch ist, sind die n Wurzeln sémtlich reell und positiv. 


12. Zahlenbeispiel: Die biegekritischen Drehzahlen einer Liufergruppe mit fiinf Stiitzen. Fiir 
eine starr zusammengekuppelte Laufergruppe (Abb. 13, Tab. 1)’, die fiinffach gelagert ist, sollen 


Abschn.: 9 
Feld: 


8(D) 7 &(C) 5 4(6)3 2A) 7 
2 bo Te BBY 5) 


3 
2 


a 
Stelle: 10 9 


1) die vier untersten Eigenfrequenzen der Biegeschwingungen berechnet werden, wobei die Stiitzen 
als vertikal starr anzusehen sind. Es soll 2) untersucht werden, welcher KinfluB auf die Lage der 
Eigenfrequenzen zu erwarten ist, wenn die Stiitzen nicht starr, sondern elastisch sind. 


Wir verwenden die folgenden Matrizen: 


fur die Felder A,, C,, D, die Rayleigh-Matrix (8) 
fursdie Welder 7B; Cy, C;, G.. D, die Elastizitats-Matrix [(27) mit 6 = 0] 
fiir die Felder A,, A,, C, die Massen-Matrix [(28) mit 0 = 0] 
fiir die Felder A,, A;, Dg die Produkt-Matrix &,- ©, aus [(27) mit 6 = 0] 
fiir die Stiitzen 1, 3, 5, 7, 9 die Stiitzen-Matrix © [(47) mit Jy =0 und m = 0| 
Tabelle 1 
5 Steifigkeit 
ae Feld Linge (m) | Masse = ee 
A I if 0,720 0 2 
2 0,665 0 7,0 
2, 0 0,255 — 
3 4,000 1,770 36,0 
4. 0 0,255 — 
5 1 0,665 0 7,0 
wD 0,750 0 1 
B 1,045 0 0,8 
G 1 0,942 0 0, 
2 0 0,204 = 
3 0,815 0 a2 
4 1,138 0,255 5,0 
5 0,915 0 Os52: 
D 1 15335 0 0,44 
2 1,610 0,357 Deo 
at 0,715 0 0,55 
2) 0,135 0 0,07 


Vel. z. B. I. P. Den Hartog, Mechanische Schwingungen, S. 143 f. oder R. Zurmiihl, Matrizen, S. 348 f. 
? R. Zurmiihl, Praktische Mathematik, S. 136, Berlin 1953. 


_ * Das System ist hier stark vereinfacht. Dieselbe Laufergruppe wurde im IPM Darmstadt mit einer Auf- 
teilung in 47 Felder berechnet, wobei der interessierende Frequenzbereich in 50 Intervalle unterteilt wurde. 
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Wir ermitteln die Abschnittsmatrizen fiir ~ = 150, 200, 250, 300, 320, 340, 400. Die Matrizen 
fiir die masselosen Felder sind frequenzunabhangig, wir brauchen sie nur cinmal aufzustellen. Fiir 
die Funktionen (4) der Si-Matrizen verwenden wir die dort angegebenen Reihen, die sich auf der 
Rechenmaschine bequem bestimmen lassen. Tabelle 2 zeigt das Rechenschema fiir die Abschnitts- 


matrix 2 bei w = 320/sec. Dabei sind die Bezugsgréfen 1 = 1m und EI = 10! mt gewahlt. Fir 
die Lager ware dann nach (46) [9 = 10-' c mit c in t/m in G4 (48) einzusetzen. Wir setzen besser 
Ig = mit ¢ in ¢/0,1 mm um angenehmer rechnen zu kénnen. 

Aus den Abschnittsmatrizen bilden wir die zugeordneten Determinanten-Matrizen [nach (38)], 


ee a (Tabelle Je Die Frequenzgleichung lautet mit M = 0, Q = 0 am linken und rechten 


xe le 0 (87a) 
wobei 


X4 = G4. 4. G4. 64. S4.H4. S4. 4. 4 (87b) 


ist. Ferner ist das fiir y; = 1 am rechten Rand erforderliche Erregermoment M,) am linken Rand 
(M,, = 0 fiirw =.) [vgl. (64) ] 
A 


Vi : s 
“ac Ay, By, Cy Dy, 


10¢mt. (88) 


Tabelle 4 zeigt diese Rechnung fiir ~ =320, wobei innerhalb der einzelnen Matrizenelemente 
nach Potenzen von J’ = c geordnet ist. 


Tabelle 2: Rechenschema fiir die Abschnittsmatrix YX. 


Ci == HO3 ial == Tis Ue CS AaOie nl 


Bs = 43 S = 1,271724 
74 = 32,262 c = 0,545096 
48 = 1040,85 s = 0,173094 
1 2,37014 5,08690  0,24257 ~—0, 30811 
1 1,39610 2.37014  0,14148 024257 
: 39,5190 50,1967 2.37014 + ~—«5,08690 
2.6112 i , | 230497 39,5190 1,39610 2.37014 | | 
1 l 1,385 0,2253 
1 0 1 0,3161 
i 0 0 1,385 
2,6112 1 f : y : ‘ 
l 1,415 0,3374 0,2452 2,37014 5,08690  0,24257 ~—«0, 30811 
0 1 0,7202 0.6817 1,39610 237014 0,14148 —0~,24257 
0 0 1 1,415 39,5190 50,1967 2,37014  —5,08690 
0 0 0 1 29,2387 52,8019 2,029503,17468 | 
5 
24,8487 38,3240 1,74008 —_3,14610 1 1,38500 0,22530 
49,7897 74,5169  3,23196 ~—«6,07033 0 fl 0,31610 
80,8917 124,911 5.24188 957907 0 1 1,38500 
29,2387 52,8019 2,02950 317468 || 2,6112 3,61651 158830 | 
33,0638 84,1174 25,1196 
65,6405 165,429 48.8902 
X= | 105,905 271,589 80,1839 
37,5284 104,779 31,1314 


Tee | wm = 320 | 


1 Die Matrizen und 4-Matrizen brauchen dabei nicht voll besetzt zu sein, wie das die Tabellen 2 bis 4 zeigen. 
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Tabelle 3. Determinanten-Matrix! 4 fiir o = 320. 


— 51,80 , , —43,00* 
RTS ; —17,36* 
+ 307,6 ; 831 
+ 307,6 ; = 13.41 
4+ 669,4 : 427,38 
+ 904,24 ; 453,38 


Tabelle 4. Eigenwert-Determinante A fiir elastische Lager. 


6. Zeile der Matrix (S4- 94) (= Dé,-c + Dé; +1): 


— 9,251 — 30,98 — 21,71 — 21,71 — 15,95* — 1,804* |-¢ 
+ 2,867 ae Bly — 1,457 — 1,456 — 11,45 — 22,27 mal 


6. Zeile der Matrix (G4- D4- G4): 


— 15,95* — 1,804* 0 0 0 0 s¢ 
— 20,70 — 52,85 — 21,71 — 21,71 — 15,95 — 1,804 °C 
+ 2,867 + 95,392 — 1,457 — 1,456 — 11,45 — 22,27 1 


——— 


6. Zeile der Matrix (64-D4- €4-€4- ©4 BY4- SA): 


—— 716,04 — 43,09* 0 0 0 0 : e 
+ 4248 + 1961,7 — 28,81 — 28,86 — 76,03 — 43,08 $ Le 
-+ 12335 2570 + 3786 a7 3785 + 4055 + 1863 165 
— 837.6 + 4087 + 5733 —- 5732 + 9940 + 7820 c 
= SOO ae NOT = 1160 == 1760 —— AS == SBT al 

6. Zeile der Matrix 

(G4- 4- 64.64. G4- B4- S4- Y4): A/10? = 
+ 864,5 + 6603* -c! + 6,603 ag 
— 187,70 - 10 — 338 - 103 ~ 3 — 337,1 -ct C> Ow: 
+ 6986 - 103 == 1393 = 102 -C == ISyey! -¢ [af a a cual 
4+ 17580 -108 + 256,7-103 |-c 7243 | |e = eee 
— 2414 - 108 + 157,8 - 103 olf + 17740 “Cc 

ee ee ee eee Aa on} Oy == &'pXi) 
1. Spalte Pacerie sian at aa [v= 20] 

Fir wm = 320 wird also bei starren Lagern J = + 6603, oder Mj, (y, = lmm/m) = — 47,4 mt; 


bei elastischen Lagern ergibt sich fiir A ein Polynom fiinften Grades. Wir kénnen hieraus A fiir eine 
bestimmte Federzahl ¢ bestimmen, oder wir kénnen auch diejenige Federzahl c als Wurzel der 
Gleichung A= 0 ermitteln, fiir diew) = 320 Eigenfrequenz wird. So erhalten wir hier ¢ = 550 t/mm 
(die tibrigen Wurzeln sind nicht reell), Fiihrt man die Rechnung fiir weitere Frequenzen w durch, 
so lassen sich aus den Kurven A(w) oder M,)(w) die Nulldurchgange und damit die Eigenfrequenzen 
bei starren oder elastischen Lagern finden; oder man kann c(w) aufstellen, was fiir die Beurteilung 
konstruktiver Anderungen von besonderem Vorteil sein kann. 


13. SchluBbemerkung. In dem vorliegenden Teil I wurden die Grundziige der Anwendung des 
Matrizenkalkiils auf Balkenbiegeschwingungen entwickelt. In einem folgenden Teil II soll das An- 
wendungsgebiet der Matrizenrechnung erweitert werden auf die Bestimmung von Rahmen- 
schwingungen, bei denen im allgemeinen Fall die Biegeschwingungen gekoppelt sind mit Langs- 
oder Torsionsschwingungen. Diese Kopplung fiihrt zu Matrizengleichungen sechsten Grades. Sie 
ermoglichen es, offene und geschlossene Rahmen zu erfassen, unter gewissen Voraussetzungen auch 
Stockwerkrahmen und Rahmenfundamente. 


(Eimgegangen am 18. November 1954.) 
Anschrift des Verfassers: Dr.-Ing. H. Fuhrke, Fellbach (Wiirtt.), Heinrich-Heine-Str. 8. 


Fir den Fall starrer Lager werden nur die mit * bezeichneten Elemente benotigt. 
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Uber ein Verfahren zur Bestimmung der Biegemomente von Platten 
unter Einzellasten 


Von S. Woinowsky- Krieger 


1. Einleitung. Die Differenzenrechnung bietet in ihrer Anwendung auf Probleme der Platten- 
biegung den wohlbekannten Vorteil einer weitgehenden Anpassungsfahigkeit an die vorgeschriebenen 
Randbedingungen der Platte; demgegeniiber steht ihr Nachteil, keine Ergebnisse von allgemeiner 
Giltigkeit zu liefern. 

Es handele sich insbesondere um die Bestimmung der maximalen Feldmomente einer diinnen 
Platte im Falle einer Einzelkraft, die sich etwa gleichférmig auf eine gewisse Druckflache verteilen 
mége. Will man diese Aufgabe unter Anwendung eines numerischen Verfahrens, insbesondere der 
Differenzenrechnung lésen, so miissen vorerst spezielle Annahmen iiber die Begrenzung der Last- 
flache gemacht werden. Die auf dieser Grundlage gewonnenen Ergebnisse lassen sich dann aber 
nicht auf Druckflachen einer anderen Form oder GréBe iibertragen. 

Die im folgenden entwickelte einfache Methode gestattet es nun, die einer numerischen Rechnung 
entnommenen und fiir einen Einzelfall der Belastung giiltigen Momentenwerte zum Aufbau ziem- 
lich allgemeiner Ausdriicke fiir Biegemomente in der Mitte einer rechteckigen oder kreisférmigen 
Lastflache nutzbar zu machen und eine wiederholte Anwendung des numerischen Verfahrens auf 
diese Weise zu umgehen. 

2. Grundlage des Verfahrens. Die Belastung einer isotropen diinnen Platte von konstanter 
Biegesteifigkeit D bestehe aus einer Einzelkraft P, die sich gleichférmig und zwar zunachst auf 
die Flache eines Rechteckes mit den Seitenlangen s und t verteilen mége. Zur Vereinfachung des 
Problems werde angenommen. da} die Biegeflache der Platte 
mindestens eine, den Mittelpunkt der Lastflache enthaltende 
Symmetrieachse besitzt. SchlieBlich miissen noch s und ¢ als 
klein im Verhaltnis zu den GrundriBabmessungen der Platte, 
sowie zu der kleinsten Entfernung der Lastflache vom Platten- 
rand vorausgesetzt werden. 

Im Sonderfall der nach Abb. 1 belasteten, allseitig frei- 
gestiitzten Rechteckplatte laBt sich zeigen, da die beiden 


WA 
W 
Biegemomente m, und m, im Schwerpunkt der Druckflache NN NN 
durch die Relationen! N N 


mw & 
4a sin — 
SoCo) a a | 
Be eee ).| 
she) 


1) 
She nee he An (u + Y) | Y 


Abb. 1. Rechteckige Platte mit einer Einzellast. 


bestimmt sind. Wahrend » hierin die Querdehnungszahl bezeichnet, gelten fiir die Gréfen ~ und py 
in guter Naherung die Ausdriicke 


p(k) =k arctg = +- + arctgk j 

(2) 
; 1 il 

p(k) =k aretg ko are tg k 
mit k =t/s. Die von gy, y und dem Logarithmus-Glied abhingigen Momentenanteile entsprechen 
dem Spannungszustand einer gewissen Kreisplatte, die in ihrem Mittelpunkt eine gleichformig 
iiber die Rechteckflache s X t verteilte Einzelkraft P tragt. Die (weiter nicht bendtigten) GréBen 
A und ys sind dagegen Korrekturglieder, die dem Unterschied zwischen jener Kreisplatte und der 
gegebenen Rechteckplatte, sowie der speziellen Stellung der Last auf dieser letzteren Platte Rech- 


1S, Woinowsky-Krieger, Ing.-Arch. 21 (1953), S. 335, 
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nung tragen. Dementsprechend hingen J und y wohl von den Verhaltniszahlen b/a und é/a, nicht 
aber von der Art der Lastkonzentration ab. 


Allen weiteren Betrachtungen sei die etwas allgemeinere Form 


1 1h B 

mm DP a gt pug) : 
1—v)P 

mx, — My es = (u ) 


der Beziehungen (1) zugrunde gelegt, wobei zwischen der frither eingefiihrten GréBe A und der 


neuen % die einfache Relation Pete 
(SUD) Ses 


besteht. Die Ausdriicke (3) gelten zunachst nur fiir die freigestiitzte Rechteckplatte nach Abb. 1. 
Hine Anderung der Grenzbedingungen dieser Platte oder gar ihrer Gestalt wiirde aber nach dem 
Vorhergesagten nur die GréBen x und yw beeinflussen, wahrend man die GréBe a als eine nunmehr 
willkiirliche Vergleichslange beibehalten kénnte. Von ihrer Wahl bleibt allerdings die GroBe ~ mit 
abhangig. 

Aus diesen Uberlegungen ergibt sich, wenn nach den Biegemomenten in der Mitte einer in ge- 
wissen Grenzen beliebigen Rechteckflache s x t gefragt wird, der folgende Vorgang zu deren Er- 
mittlung. Man berechnet zuerst die Biegemomente nach einem numerischen Verfahren fiir ein 
nach praktischen Gesichtspunkten gewahltes Wertepaar s,/a und t)/a und die gewiinschte Stellung 
der Last. Die Werte s, und t, sollten von gleicher Gré®enordnung sein und am besten etwa 10 bis 
20% der kiirzesten GrundriBabmessung der Platte betragen. Bei Anwendung der Differenzen- 
rechnung wird man auch fir ein ganzzahliges Verhaltnis von s) und tf, zu den Maschenweiten des 
értlichen (dichteren) Netzes in den entsprechenden Richtungen sorgen miissen. Hinsichtlich der 
Querdehnungszahl y sind die folgenden Falle zu unterscheiden: 

1) Die Zahl y geht in die Randbedingungen des Plattenproblems nicht ein. 

2) Die Zahl y ist in den Randbedingungen mit enthalten. 

Nur im Falle 2) wird man aufer vom Wertepaar s,/a, t)/a auch von einem bestimmten nume- 
rischen Wert von y ausgehen miissen. 

Die Rechnung wird nun fiir die Krimmungen der Platte 0?w/dx? und 0’w/dy? im Schwerpunkt 
der Lastflache gewisse numerische Werte liefern und die zugehérigen Biegemomente ergeben sich 
dann aus den Relationen 


Ow O*w 
ee ae ae (es : oh (5) 
Ow ow 


Nachdem das Ergebnis (5) feststeht, ist auch der linke Teil der Gleichungen (3) fiir das Wertepaar 
$,/4, to/a bekannt. Die zugehérigen Werte gy und wy) von pm und yp sind durch die Ausdriicke (2) als 
Funktionen von ky = t)/sy ebenfalls gegeben und nun kénnen die GréBen x und yw auf der rechten 
Seite der Gleichungen (3) eindeutig bestimmt werden. Die Biegemomente am Lastort fiir beliebige 
Werte s/a und t/a findet man endgiiltig aus den Gleichungen (3), die nunmehr rechterhand lauter 
bekannte GréBen enthalten. Die Zahl y geht nur im Falle 2) in die berechneten Werte 0’w/0x? und 
0’w/dy” wirklich ein, Wie ein Vergleich der Gleichungssysteme (3) und (5) sofort lehrt, sind auch die 
GréBen x und y nur in diesem Fall 2) von der Querdehnungszahl y mit abhangig. 

Ist die wirkliche Aufstandsflache der Last nicht ein Rechteck, sondern ein Kreis vom Halb- 
see c, so bediirfen die Gleichungen (3) einer Umformung, die sich aus der folgenden Uberlegung 
ergibt. 

Im Mittelpunkt einer eingespannten Kreisplatte vom Halbmesser a mit zentrischer, quadrati- 
scher (s = t) beziehungsweise kreisformiger Druckflache der Last gelten die Ausdriicke ! 


rep reaes saad bey ! ae a, 


Se 2 ’ 
(6 
Sie fe e 
Maret A 3 (in Ca | a 


1S. Woinowsky-Krieger, a. a. O., S. 332. 
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fiir die Momentensummen, Ist c/a <1 und vernachlissigt man demzufolge das Glied c/2 a gegen 
den Logarithmus, so fiihrt die Gleichsetzung der beiden Momentensummen zu der Beziehung 


mee 
s=t=cyV2e = W158 Cs ee) 


die die Bedingung fiir die Aquivalenz der quadratischen und der kreisférmigen Lastflache aus- 
driickt. Die Frage nach der Gréfe der Biegemomente im letzteren Falle wird somit durch das 
Einsetzen des Wertes (7) fiir s und t, sowie von p = 7/2 und y = 0 in die Gleichungen (3) erledigt. 

Liegt eine Platte von polygonaler, aber nicht rechteckiger Begrenzung vor, so wird man. zur 
Grundlage der Differenzenrechnung auch entsprechend begrenzte Lastflachen wahlen miissen. 
Beispielsweise kénnte es sich dabei um Dreiecke, Sechsecke oder Rhomben handeln. Die Ableitung 
von Ubergangsformeln von derartigen Druckflachen zu der kreisférmigen reduziert sich stets auf eine 
Quadratur’. Hiermit wird aber auch der Ubergang zur rechteckigen Druckflache gesichert, die in 
der technischen Praxis am haufigsten vorkommt. 


3. Beispiele zur Beurteilung der Genauigkeit der Differenzenrechnung. Stellt man die aus den 
Ausdriicken (1) gewonnenen Momentenwerte Werten gegeniiber, die sich durch direkte Anwendung 
strenger Lésungen, meist in Form von einfachen Reihen, 
ergeben, so driicken sich die dabei begangenen Fehler in 
Bruchteilen eines Prozentes aus, sofern etwa s/a < 0,5, 
t/b < 0,5 bleibt und die Plattengestalt nicht ausgesprochen 
langlich ist. Da aber bei der Anwendung der Methode die 
Fehlerquelle ebensogut in der Differenzenrechnung liegen kann, 
ist es nicht uninteressant, sich von der Genauigkeit dieser 
letzteren an Hand einiger Beispiele ein Bild zu machen. 

Fir die Biegemomente in der Mitte einer durch eine Ein- 
zelkraft P belasteten freiaufliegenden quadratischen Platte 
nach Abb. 2 findet H. Marcus die Werte? 


m, = my = 0,263451 P, (8) 


wobei s = t = 0,125 a und vy = 0,3 vorausgesetzt wird. Die 


a a 
a 2 


NN 


BJea 


erste der Gleichungen (1) liefert im selben Falle, da? 1 = 2,669 ¥ 
(fiir a=b=—2 &) und P= m/2 (fiir s= t) anzunehmen sind, Abb. 2. aE oats einer zentrischen 
it, == fit, == 0,201 4 (9) 


Da die Formeln (1) gerade bei der quadratischen Platte sehr genau sind, kommt der Unterschied 
von 0,9%, zwischen den Werten (8) und (9) gréBtenteils auf das Konto der Differenzenrechnung. 


Fur die allseitig eingespannte Platte, sonst aber unter den gleichen Voraussetzungen wie zuvor, 
liefert die Differenzenrechnung * 


m, = my = 0,223 P., (10) 


wahrend die notwendige Korrektur des Wertes (9) sich auf Grund vorhandener analytischer L6- 
sungen fiir die eingespannte Platte leicht ausrechnen laBt°. Stellt man dabei nur 7 Glieder, der das 
Einspannungsmoment darstellenden unendlichen Reihe in Rechnung, so ergibt sich die Einwirkung 
der Randeinspannung auf die Biegemomente in Plattenmitte zu — 0,042 P, und man findet fiir die 
letzteren Momente selbst die Werte 


m, =m, = 0,219 P. (11) 


Die Ergebnisse (10) und (11) weichen somit um 1,8% voneinander ab. Im Sonderfall von Stahl- 
betonplatten ist diese Differenz immer noch erheblich kleiner als der von der Unsicherheit des 
y-Wertes des Materials herriihrende Fehler in der Bestimmung der Momente. 


1 S$. Woinowsky-Krieger, a. a. O., S. 332. 

2 H. Marcus, Die Theorie elastischer Gewebe und ihre Anwendung auf die Berechnung biegsamer Platten, 
S. 70, Berlin 1932. 

3S. Woinowsky-Krieger, a. a. O., S. 336, 334. he 

4 Der von H. Marcus, a. a. O., S. 161, angegebene Momentenwert 0,163094 P ist indessen fehlerhaft und 
das ihm zugrunde liegende Moment.0,20352 P der freigestiitzten Platte unverstandlich. Der EinfluB der Ein- 
spannung ist nach Marcus, — 0,040426 P (S. 161); zieht man ihn vom Wert(8) ab, so erhalt man das Re- 
sultat (10). 

5 Dana Young, J. Appl. Mech. 6 (1939), S. A—114. 
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4. Die nur in den Eckpunkten abgestiitzte quadratische Platte mit zentrischer Einzellast. Die 
Differenzenrechnung ist zur Behandlung dieses Problems sehr wohl geeignet; es mége daher auch 
als Grundlage fiir die Anwendung des im Abschnitt 2 beschriebenen Verfahrens gewahlt werden. 

Die Biegemomente im Mittelpunkt « = y = 0 der quadratischen Platte mit freien Randern 
nach Abb. 2 berechnen sich fiir sy = t) = 0,125 und y = 0,3 zu! 


m, = m, = 0,34461 P. (12) 


Nach Einsetzen dieser Werte in die erste der Gleichungen (3) erhalt man wegen yy = 77/2 den Wert 
x = 4,767, wahrend die zweite Gleichung mit y) = 0 den Wert 4 = 0 liefert. 

Geht man mit den so erhaltenen Werten x und y und auBerdem noch mit y = 0,3 in die Glei- 
chungen (3) ein, so findet man fiir die Lastflache in der Form eines beliebigen Rechteckes die Be- 
ziehungen 


2 
m, + m, = 0.1034 (in a=) =e 0,492,|P | 
(13) 
m,— my =0,0557 yp P. 


Die hier enthaltenen GréBen ~ und y sind nach wie vor durch die Ausdriicke (2) bestimmt. Die 
Losung ist hiermit hergestellt, und es sind nur noch einige Worte iiber die Rolle zu sagen, die darin 
die Querdehnungszahl spielt. 

Da die Randbedingungen der Platte nur fiir y = 0,3 erfiillt sind, so wurde eben dieser Wert der 
Querdehnungszahl auch fiir die Faktoren 1 + » und 1 —y in den Gleichungen (3) beibehalten; er 
liegt somit auch den Beziehungen (13) zugrunde. 

Mit abnehmendem » vergréBern sich die Durchbiegungen der Platte nicht unerheblich, somit 
auch die Kriimmungen; die Biegemomente in Plattenmitte werden aber durch die Zahl vy nur wenig 
beeinfluBt. So sind beispielsweise die Biegemomente in der Mitte einer gleichférmig durch eine 
Last Q belasteten, nur an den Ecken aufruhenden quadratischen Platte gleich 


m, = my, = 0,1090Q, (14) 


wenn man vy = 0,3 annimmt? und sie haben den Wert 0,1109 Q falls mit y = 0,25 gerechnet wird?®. 
Weicht also die tatsachliche Querdehnungszahl von dem Wert y = 0,3 ab, so erscheint es richtiger, 
die Biegemomente trotzdem nach den Gleichungen (13) zu berechnen, als etwa die Momentensumme 
in (13) mit dem Faktor (1 + »)/1,3 zu vervielfachen. 


5. Quadratische Fundamentplatte bei gleichformiger Bodenpressung. Es ist naheliegend, die 
soeben gewonnenen Ergebnisse auf den Fall einer zentrisch belasteten quadratischen Fundament- 
platte von konstanter Biegesteifigkeit anzuwenden. Trotzdem man iiber verschiedene Theorien 
verfiigt, die, je nach der Beschaffenheit des Baugrundes, die Verteilung der Bodenreaktion iiber die 
Fundamentsohle mehr oder weniger zutreffend angeben, wird in der Praxis noch sehr oft mit einer 
gleichférmigen Bodenpressung bei zentrischer Auflast gerechnet. Fiir die Bemessung der Platten- 
bewehrung ist dies fast stets eine sichere Annahme. Vom Standpunkt der Bodenziffertheorie 
wiirde sie einer absolut steifen Platte entsprechen. 

Die Auflast P mége sich nach wie vor gleichférmig tiber eine Rechteckflache nach Abb. 2 ver- 
teilen. Die zugehdérigen Biegemomente in Plattenmitte ergeben sich dann durch eine einfache 
Superposition der Werte (13) mit dem Wert (14), wo offenbar Q——P zu setzen ist. Die Momenten- 
summe in den Ausdricken (13) verringert sich nun um den Wert 0,218 P, und das Endergebnis 
lautet 


Mm, + my = (0.1034 (in sos —9) 2 0,274 Pe 
m,—my = 0,0557y P. 


(15) 


Die Ausdriicke gelten streng genommen nur fir y = 0,3, praktisch aber fiir ein beliebiges vy. Wie 
oben ausgefiihrt, ist der EinfluB der Querdehnungszahl auf die Biegemomente in Plattenmitte in 
der Tat nur gering und eine iibertriebene Genauigkeit ware gerade bei der Bemessung von Funda- 
mentplatten fehl am Platze. 


1H. Marcus, a. a. O., S. 188. 
2H Marcus aan Os lols 
3 B. Galerkin, Gesammelte Werke Bd. 2, S. 22, Moskau 1953 (Russisch). 


XXIII. Band 1955 i -Kri : U i i 5 
Woinowsky-Krieger: Uber ein Verfahren zur Bestimmung der Biegemomente usw. 300 


Verteilt sich die Auflast anstatt auf eine Rechteckfliche gleichférmig auf cine Kreisflache mit 
dem Halbmesser c, so hat man in den vorstehenden Gleichungen s und t durch 1,753 ¢ zu ersetzen 
und schlieBlich noch py = z/2 und y=0 anzunehmen. Die Biegemomente in der Plattenmitte 
ergeben sich dann gleich 


mM, = my = (0,1034 In — 0,038) P. (16) 
Der fiir die Anordnung der Stahlbewehrung wichtige Verlauf der Biegemomente langs der Symme- 


trieachsen ee 0 und y = 0 der Platte ist auf Abb. 3 fiir c = 0,05 a undc = 0,1 a als Beispiel dar- 
gestellt. Die Auftragung der Momentenkurven wird durch die folgenden Uberlegungen erleichtert ; 


0056P 
QO56P 


Abb, 3. Biegemomente m, in den Schnitten x = 0 bzw. y = 0 einer quadratischen Fundamentplatte unter zentrischer, 
auf eine Kreisfliche vom Halbmesser ¢ verteilter Einzellast P und gleichformiger Bodenpressung. Links Momente fir 


c/a = 0,10, rechts fir c/a = 0,05; » = 0,3. 

1) Die Biegemomente in der Nahe der Plattenrander werden durch die Art der Lastkonzen- 
tration nur wenig beeinfluBt solange die Druckflache klein bleibt im Verhaltnis zu der Flache der 
ganzen Platte. Die betreffenden Momentenwerte lassen sich also den Ergebnissen der Differenzen- 
rechnung entnehmen. 

2) Am Rande des Druckkreises ist das tangentiale Biegemoment um etwa v P/4 7 und das 
radiale Moment um etwa P/4 7 kleiner als der Scheitelwert beider Momente m, = m, in der Plat- 
tenmitte}, 

3) Der Durchschnittswert der tangentialen Spannungsmomente, also etwa der Momente m, im 
Schnitte x = 0 1aBt sich aus dem statischen Moment der an einer Plattenhalfte angreifenden Krafte 
in einfacher Weise berechnen. 

6. Schlu8bemerkung. Allen bisherigen Ausfiihrungen lag die Annahme einer gleichférmigen 
Verteilung der Last auf eine Druckfliche von gegebener Form zugrunde. Vom technischen Stand- 
punkt aus ist noch der Fall einer Einspannung der Platte, etwa einer Fundamentplatte langs der 
Randlinie einer Saule, also der Lastflache, von Bedeutung. Anstatt der Biegemomente in der Mitte 
der Lastflache miissen alsdann die Einspannmomente um die Sadule herum als die maBgebenden 
Plattenmomente angesehen werden. Im Sonderfalle einer kreisrunden Saule ist es nicht schwer von 
den ersteren auf die letzteren zu schlieBen, doch soll hier auf diese Frage nicht eingegangen werden. 
Wesentlich verwickelter liegen die Verhaltnisse bei der rechteckigen Saule.. Bei einem Material 
ohne FlieBgrenze wiirden dabei die Plattenbeanspruchungen mit der Annaherung an die Saulen- 
ecken ins Unendliche gehen. Bei einem nicht ideellen Material dagegen hatte man hier mit einem 
elastisch-plastischen Deformationszustand der Platte rechnen miissen. 


(Eingegangen am 26. November 1954.) 


Anschrift des Verfassers: Prof. Dr.-Ing. S. Woinowsky-Krieger, Quebec, Que. Canada, 
200 Ave. Marguerite Bourgeoys Apt. 8. 


1 4, Nadai, Die elastischen Platten, S. 63, Berlin 1925. 
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Nichtlineare Schwingungen und kinetische Instabilitat 
bei Saiten und Stiben 


Von E. Mettler 


1. Einleitung. Das Problem der Stabilitat erzwungener Schwingungen elastischer K6rper ist in 
letzter Zeit viel bearbeitet worden; man zahlt heute etwa 30 Veréffentlichungen iiber diesen Gegen- 
stand!, Trotzdem kann man nicht sagen, daB das Problem erschépfend gelést ist. Vielmehr muf 
auffallen, daB die theoretisch so viel erérterte Instabilitatserscheinung nie praktisch beobachtet 
worden ist (mit einer Ausnahme, die noch besprochen wird), und im Zusammenhang damit drangt 
sich die Frage auf, wie sich das merkwiirdige Phanomen eigentlich 4uBern muf oder was man bei 
einer experimentellen Verwirklichung als sichtbare Erscheinung zu erwarten hat. 

Die letztere Frage 1a8t sich natiirlich mit der bisher fast allein geférderten linearisierten Stabili- 
tatstheorie nicht vollstandig beantworten. Denn diese kann zwar bestimmte Schwingungen als 
stabil oder instabil nachweisen, sie sagt aber gerade in dem zur Diskussion gestellten Instabilitats- 
fall lediglich aus, was man nicht beobachtet (namlich die instabile Schwingung). Dagegen kann die 
lineare Theorie ihrem Wesen nach nicht angeben, welche Bewegung an Stelle der instabilen Schwin- 
gung wirklich eintritt. Das leistet nur eine nichtlineare Theorie, die allerdings bisher nicht tber 
einige wenige Ansiatze hinaus gediehen ist, wobei sie sich auf die Biegeschwingungen von Saiten 
und Staben beschrankte. Die folgenden Ausfiihrungen beschaftigen sich zunachst in Abschnitt 2 
mit einer Zusammenstellung und Deutung dieser letztgenannten Ergebnisse, die somit das Verhalten 
einer elastischen Schwingung in der Umgebung eines Instabilitétsbereichs erster Art wiedergeben. 
In Abschnitt 3 werden einige dariiber hinausgehende Untersuchungen tiber die Kippschwingungen 
des Stabes angestellt mit dem Ziel, entsprechende Ergebnisse fiir einen Bereich zweiter Art zu ge- 
winnen. 


2. Biegeschwingungen von Stab und Saite®. a) Freie und quererregte Schwingungen. 
Zur Vorbereitung des folgenden betrachten wir einen geraden Stab mit konstanter Biegesteifigkeit 
E J, Massenbelegung 4 und Lange I. Er sei beiderseits ge- 
Sabachse lenkig gelagert und trage die Querbelastung p(x, t) sowie die 
Wp hx eingepragte Langskraft P = const. Die Langenkoordinate sei 
ae wie tiblich x, die Zeit t, tiefgestellte Indizes x und t bedeuten Ab- 
leitung nach diesen GréBen. Die Querschwingung mit der Aus- 
lenkung w gehorcht dann der bekannten Differentialgleichung 
ee Te Wa ees a S We x + Wey == 0D (1) 
ax worin die Langskraft S (als Druckkraft positiv gerechnet) in 

Abb. 1. Zur Berechnung der Dehnung zwei Anteile zerfallt; die Vorlast P und die Zusatzkraft A P, 

die durch die Dehnung der Stabachse bei der Ausbiegung ge- 
weckt wird. Diese Dehnung ist nach Abb. 1 gleich (ds — dx)/dx = yl + wi—1 2 w;/2, die mitt- 


lere Dehnung somit 


1 
1 Ww; 
Cmittel = ie [ 5° dx . 


j 0 
Der Ausdruck 1 


AP © EF enn = "yy | wh dx (2) 
0 


ist also eine plausible Naherung fiir A P, und (1) geht mit S = P — AP in eine Integrodifferential- 
gleichung tber. 


1S. die Literaturzusammenstellung bei E. Mettler, Publ. sci. techn. Min. Air Paris 281 (1953), S. 77, wo das 
Schrifttum bis 1951 genannt ist. 

* Die in diesem Abschnitt vorgebrachten Dinge stammen im wesentlichen aus einem Vortrag ,,Bericht 
iiber die Theorie der nichtlinearen Schwingungen elastischer Kérper“, den der Verf. auf der GAMM-Tagung 
1952 in Braunschweig hielt und der damals nicht mehr veréffentlicht wurde. 
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Das Integralglied (2), das hier eine Deutung als Mittelwert gefunden hat, wurde schon von 
mehreren Autoren mit verschiedener Begriindung eingefiihrt. Zuerst gab es Kirchhoff bei seiner 
Berechnung der freien Schwingungen einer ,,schwach gespannten” Saite an. Die Gleichung dieses 
Problems erhalt man aus (1) und (2) fir E J = 0 und p =O. Ferner benutzen mehrere neuere 
Arbeiten? dasselbe Glied (2), von denen hier nur die Untersuchungen von Schlesinger und Burgreen 
hervorgehoben seien. Diese behandeln die freien Schwingungen der Saite bzw. des Stabes experi- 
mentell und bestatigen die von dem nichtlinearen Ausdruck (2) herrithrende Amplitudenabhangig- 
keit der Eigenfrequenz durch Messung, was umgekehrt als Bestiétigung der Brauchbarkeit der 
Naherung (2) aufgefaBt werden darf. 

Von Bedeutung fiir die weitere Diskussion ist ferner der Fall der erzwungenen Schwingungen 
p =+ 0. Dieser ist bisher fiir die spezielle Querbelastung 


- Ux 
P eg (Ore 


behandelt worden? und fiihrt mit dem Ansatz 


w = v(t) sin ae (3) 
von (1), (2) auf die Duffingsche Gleichung 
v+a?v(1+yv") =f cost, (4) 


worin @, die erste Biegeeigenfrequenz des Stabes, § und y Konstante bedeuten. Von der Differen- 
tialgleichung (4) ist als wichtigstes Integral eine nicht harmonische periodische Lésung bekannt, 
deren Amplitude durch das charakteristische Am- 
plituden-Frequenz-Diagramm nach Abb. 2 wieder- 


gegeben wird. Auf den ausgezogenen Asten ist nee ee 
die Lésung stabil, auf dem gestrichelten Ast in- at ye 
stabil. Zwischen den stabilen Asten finden die be- va 


kannten Kipperscheinungen statt. 


Amplitude 


b) Langserregte Schwingungen aq) Die 
Differentialgleichung. In der vorliegenden 
Arbeit steht allerdings, wie in der Einleitung D) 
auseinandergesetzt, nicht der im vorhergehenden , : ; Pe eae 

2 Abb. 2. Amplituden-Frequenz-Diagramm des einfachen nichtlinearen 
Abschnitt behandelte Fall zur Debatte, sondern Schwingers und des querangeregten Stabes. 
der Stab mit pulsierender Langskraft, denn er 
zeigt die charakteristischen Instabilitaten. Um auf seine Gleichung zu kommen, liegt es nahe 
auch hier die Gleichgewichtsaussage (1) heranzuziehen, darin aber p = 0 zu setzen und dafiir die 
Bedingung P = const aufzuheben, was natiirlich nur méglich ist, wenn man ein Lager langsver- 
schieblich macht. Setzen wir also wieder S = P— AP und jetzt P = P, + P , coswt, so kommt 
nach (1) und (2) 


I 
EI oenen | Pat P, cos wt — ws) an + pan = 9. (5) 
\ 0 

Die Giiltigkeit dieser Gleichung ist durch die bekannte — itibrigens auch fiir (4) bestehende — Be- 
dingung eingeschrankt, daB fiir die Langsschwingungen im Stab keine Resonanz entstehen darf, 
was erfiillt ist, wenn die Erregerfrequenz w geniigend tief unter der kleinsten Langseigenfrequenz 
liegt. 

Fiihrt man in (5) den Ansatz (3) ein,* so entsteht nach kurzer Rechnung mit den Abkiirzungen 


TU 2 Pp— Py eS = a 
aa(G) ao: pep 2 P;r= i Ihe 
yeh i en eee 
Viera AS], (PPE fur den Stab, y =—( | 4 P, fiir die Saite , 


1G. Kirchhoff, Vorles. tiber Math. Physik (Mechanik), Leipzig 1876, S. 446. 

2 K. Schlesinger, Z. techn. Physik 12 (1931), S. 33; G. F. Carrier, Quart. appl. Math. 3 (1945), S. 157 und 
7 (1949), S. 97; S. Woinowsky-Krieger, J. appl. Mech. 17 (1950), S. 35; D. Burgreen, J. appl. Mech. 18 (1951), 
8. 135. 

3 E. Mettler, Z. £. angew. Math. u. Mech. 31 (1951), S. 263. Die Rechnung ist dort fir P = 0 durchgefiihrt, 
was aber vollig unwesentlich ist. : ; 

4 Es geniigt, wenn wir uns hier auf die Grundschwingung beschranken. 
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die beherrschende Differentialgleichung 

v+owv(l—ecoswt+yv) =0, (6) 
deren Verwandtschaft mit (4) in die Augen fallt. Gleichung (6) wurde fiir den Stab auf einem 
anderen Wege von Weidenhammer' sorgfaltig abgeleitet; seine Gleichung (33) stimmt mit (6) iiber- 
ein, wenn man die Dampfung weglaBt. Dies zeigt erneut die Brauchbarkeit des Integralausdrucks 
(2) fiir die Langskraft, den wir daher spater (Abschnitt 3) ohne weitere Priifung auch fiir ein anderes 
Problem verwenden wollen. Ferner hat Raman? die Beziehung (6) bei seinen Untersuchungen iiber 
die schwingende Saite mehrfach verwendet. ; 

f) Ergebnisse der linearen Theorie. Wir diskutieren nun die Lésungen von (6) zunichst 
unter der Voraussetzung der gewohnlichen Stabilitatstheorie, da sehr kleine Ausschlage vorliegen 
und deshalb das nichtlineare Glied y v2 weggelassen werden darf. Dann geht (6) in die Mathieu- 
Gleichung iiber, deren Lésung folgende Aussagen ergibt: ee. 

Fiir Wertepaare w, ¢, deren Bildpunkte in der (w, é)-Ebene auBerhalb von gewissen Instabilitats- 
bereichen liegen, ist die gestreckte Stabform, oder, wie man auch sagen kann, die reine Langs- 
schwingung stabil. Gibt man dem Stab eine anfangliche Biegeauslenkung, so schlieBen sich be- 
schrinkte Biegeschwingungen an, die von der in (6) nicht beriicksichtigten, aber in Wirklichkeit 
stets vorhandenen Dampfung zum Abklingen gebracht werden. Die Amplitude wird also Null. 


stabil 


— — /nstab/| 


(a) (a7) Y, 
a 2 Db Y 0 


Abb. 3. Darstellung des Zusammenhangs zwischen der Schwingungsamplitude A, der Erregerfrequenz w und der dimensionslosen Amplitude ¢ 
der erregenden Kraft bei langsangeregten Biegeschwingungen und Kippschwingungen. 


M =20, fiir Instabilitatsbereich erster Art, 
@y = @, + wy fiir Instabilitatsbereich zweiter Art. 


Fir Punkte (w, ¢), die gerade auf der Grenzlinie eines Instabilitatsbereichs liegen, gibt es perio- 
dische Biegeschwingungen mit unbestimmt bleibender, endlicher Amplitude, die von der Dampfung 
nicht beseitigt werden. Die Dampfung bewirkt, sofern sie klein ist, lediglich eine unbedeutende 
Verlagerung der Grenzlinien. Die Lésungen von (6) sind jetzt die Mathieu- Funktionen. 

Befindet man sich mit den Parameternw, ¢ schlieBlich innerhalb eines Instabilitatsbereichs, 
so ist die gestreckte Stabform instabil und die entstehende Biegeschwingung hat keine endliche 
Amplitude, sondern zeigt anwachsende Tendenz. 

Es handelt sich bei der Mathieu-Gleichung um sogenannte Instabilitatsbereiche erster Art. Zur 
Erérterung unsres grundsatzlichen Problems beschranken wir uns auf den ersten und gréBten Be- 
reich, Auf seinen Grenzen haben die periodischen Lésungen die Frequenz w/2. Wir geben den Be- 
reich wieder in Abb. 3, wo der eben erlauterte Zusammenhang zwischen w, ¢ und der Amplitude A 
der Biegeschwingung bildlich dargestellt werden soll. Abb. 3a zeigt den Bereich in der tiblichen 
Darstellung in der (w, ¢)-Ebene. Er hat eine Spitzeim Punkt ¢e = 0,~ = Wy) = 2 @, und die Grenzen 

ORs ie ey 

O; oes) een (7) 
Abb. 3b gibt eine raumliche Darstellung im (w, ¢, A)-System, u. zw. gehoren die punktierten Linien 
zum linearen Fall. Da die Amplitude auBerhalb des Bereichs Null ist, auf den Grenzen unbestimmt 
wird und innerhalb des Bereichs nicht existiert, kann man sie nur durch zwei lotrecht stehende 
Wandflachen darstellen, die sich iber den Grenzkurven der (@, €)-Ebene ins Unendliche erstrecken 
und deshalb abgebrochen gezeichnet sind. 


1 F, Weidenhammer, Ing.-Arch. 20 (1952), S. 315. 


> C. V. Raman, Phys. Rev. 35 (1912), S. 449 und Proc. Indian. Assoc. for the Cultivation of Sci. 6 (1912), 
S. 1 sowie 11 (1914), S. 43. 
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Abb. 3c zeigt einen Schnitt e = const durch das System Abb. 3b. In dem so entstehenden 
Amplituden-Frequenz-Diagramm stellen die beiden punktierten lotrechten Geraden und die aus- 
gezogenen Teile der w-Achse das vor, was man in der Schwingungslehre die ,,Resonanzkurve“ oder 
auch den ,,Frequenzgang der Amplitude“ oder englisch » response curve“ nennt. 


y) Nichtlineare Biegeschwingungen. Nimmt man nun in (6) das nichtlineare Glied y v? 
mit, so bleibt natiirlich die gestreckte Stabform im unveranderten Bereich instabil. Aber die Aus- 
sage tiber die Amplitude der Biegeschwingung ndert sich, und zwar derart, daB jeder Frequenz w 
endliche, beobachtbare Amplitudenwerte zugeordnet werden, womit wir also der Antwort auf die 
in der Einleitung gestellten. Frage niher kommen, 

Es 148t sich streng zeigen!, dafs (6) auch fiir y + 0 periodische Lésungen besitzt, und diese 
lassen sich durch Stérungsrechnung bequem berechnen, wobei diese Rechnung dadurch eindeutig 
zu machen ist, da8 man vorschreibt, die Lésung solle fiir y =0 in die entsprechende Mathieu-Funk- 
tion tibergehen. Fiir den bisher besprochenen Instabilitatsbereich haben die periodischen Lésungen 
wiederum die Frequenz w/2, und zwischen w, ¢ und der Amplitude A der Biegeschwingung ergibt 
sich die Beziehung 

@ 


s € 3y See 2 
alee ae 9 = A Lh + A ee (8) 


Die GréfBe A ist allerdings nur eine Naherung fiir die wirkliche Amplitude?, die ja bei einer nicht 
harmonischen Schwingung nicht unmittelbar anzugeben ist. Doch geniigt es fiir unsre Zwecke 
vollig, A als Amplitude anzusehen. 

Wenn wir nun die Beziehung (8) zwischen w, ¢ und A geometrisch deuten, so finden wir zunachst, 
da} (8) fur y = 0 in (7) ttbergeht. Fir y + 0 wird also Gleichung (8) im Koordinatensystem der 
Abbb. 3b durch zwei (ausgezogen gezeichnete) Flachen dargestellt, die auf den alten Grenzlinien in 
der (w, €)-Ebene fuBen, aber nicht mehr lotrecht aufsteigen, sondern nach rechts abbiegen. Die 
Wande des iiber dem Instabilitatsbereich stehenden Keiles neigen sich nach rechts derart, daB die 
linke Wand den Bereich wie ein Dach iiberdeckt. Besteht Instabilitat der gestreckten Stabform 
[Punkt (w, ¢) im Instabilitatsbereich], so tritt eine Biegeschwingung der Frequenz w/2 ein, deren 
Amplitude durch die Héhe A der oberen Flache iiber der (w, ¢)-Ebene gegeben ist. Steigert man 
nun w soweit, dai der Punkt (w, e) die rechte Grenzlinie des Instabilitatsbereichs uberschreitet, so 
bleibt die Biegeschwingung mit steigender Amplitude erhalten. 

Sehr wichtig ist ferner, daB die Schwingung, die zu der unteren Flache in Abb. 3b gehért, zwar 
theoretisch existiert, aber, wie die genauere Untersuchung zeigt, nicht stabil ist, also praktisch 
nicht auftritt. Die zur oberen Flache gehérende Schwingung dagegen ist stabil. 

Abb. 3c zeigt auch fiir den nichtlinearen Fall einen Schnitt ¢ = const aus Abb. 3b und weist 
iiberzeugend auf die enge Verwandtschaft der Amplituden-Frequenz-Diagramme (Resonanzkurven) 
des quer- und des langsangeregten Stabes hin. Abb. 3c erscheint als ein eckiges Gegenstiick von 
Abb. 2. 

Die schwingende Saite ist das alteste und einzige® Objekt, an dem man bisher Versuche und 
Beobachtungen iiber langsangeregte Biegeschwingungen angestellt hat. Nach den klassischen 
Versuchen von Melde* haben sich hauptsachlich Raman und Krishnaiyar® mit diesen Saitenschwin- 
gungen beschaftigt. Der Verf. hatte ferner Gelegenheit, Versuche von Uhde® zu sehen, welche die 
Charakteristik der Abb. 3c deutlich zeigten, wobei sich lediglich gegenitiber Abb. 3c eine erheblich 
gréBere Mannigfaltigkeit von Kurven ergab, da nicht nur die Grundschwingung der Saite in der 
Nahe des ersten Instabilitatsbereichs zu der besprochenen Erscheinung fiihren kann, sondern ebenso 
eine Oberschwingung in der Nahe eines héheren Instabilitatsbereichs. Die Versuche zeigen ferner, 
da die Verwandtschaft zwischen der langsangeregten Saite und dem Duffingschen Schwinger noch 
weiter geht. So wie bei letzterem die Schwingung bei einer bestimmten maximalen Amplitude, die 
von der Dampfung abhangt, abreifst und “om oberen auf den unteren Ast der Kurve Abb. 2 »kippt™, 
liegen die Dinge auch bei der Saite. Krishnaiyar hat gezeigt, daB auch bei der Saite der Punkt 


1 Die folgenden mathematischen Ergebnisse stammen von Weidenhammer a. a. O. 

2 A ist der Koeffizient des ersten Fourier-Glieds der hier nicht wiedergegebenen Losung. 

3 I. Utida und K. Sezawa, [Rep. Aeron. Res. Inst. Tokyo Imp. Univ. 15 (1940), S. 193] haben zwar Ver- 
suche mit langsangeregten Staben angestellt, doch lassen sich Versuchsbedingungen und Ergebnisse nicht in 
unseren Zusammenhang einordnen. 

4 F’, Melde, Pogg. Ann. 109 (1860), S. 193 und 111 (1860), S. 513. 

5 C. V. Raman, a. a. O. und N. C. Krishnaiyar, Phys. Rev. 14 (1919), S.490 sowie Phil. Mag. 43 (1922), S.503. 

6 A. Uhde, Erregung von Seilschwingungen durch axiale Krafte, Diplomarbeit Techn. Hochschule Braun- 
schweig (Wohlerinstitut) 1948. 
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des AbreiBens von der Dampfung bestimmt wird. Es ist nicht daran zu zweifeln, daB sich Stabe 
mit pulsierender Axiallast ahnlich verhalten. 


3. Kippschwingungen des I-Tragers. a) Die Differentialgleichungen. Nachdem im vor- 
hergehenden Abschnitt ausschlieBlich Instabilitaten erster Art behandelt worden sind, wird jetzt 
die Frage aufgeworfen, wie sich eine elastische Schwingung in der Nahe einer Instabilitat zweiter 
Art verhalt. Als konkretes mechanisches Problem, das am kiirzesten zu einer solchen Instabilitat 
fiihrt, wahlen wir den I -Trager in Gabellagerung mit Axiallast S und pulsierenden Endmomenten 
M = M, coswt nach Abb. 4. Bezeichnen wir mit 


w(x, t) die waagerechte Verschiebungskomponente der Stabachse, 

(x, t) den Verdrehwinkel des Querschnitts, 2.3 
EJ die Biegesteifigkeit des ganzen Querschnitts bei Biegung um die y-Achse (Lotlinie), 
E Jy die Biegesteifigkeit eines Flansches bei Biegung um die y-Achse, 

GJ, die Torsionssteifigkeit, 


a die Querschnittshéhe, 
Up den polaren Tragheitsradius des Querschnitts, 
ll die Massenbelegung des Stabes, 
7 
3 
pa 
Abb. 4, Der kippende Stab. ! : 


so gehorcht die Bewegung des Stabes den Differentialgleichungen! 


9 
Bn Oia 6 Jib Meee 0) 


mit den Randbedingungen 


WO) =) =O 5 =a = fir x«=0 und x=l. 


Dazu tritt noch eine Gleichung fiir die Bewegung in der lotrechten Lastebene, die von (9) unabhangig 
ist und hier nicht interessiert. Voraussetzung fiir die Giiltigkeit von (9) ist, da® die Erregerfre- 
quenz w wesentlich unter der kleinsten Eigenfrequenz der Biegeschwingung in y-Richtung liegt?. 

Wir wollen ferner voraussetzen, daB die Langskraft S keinen eingepragten Lastanteil hat, 
sondern erst durch die Verlangerung der Stabachse infolge der Ausbiegung w bei festgehaltenen 
Enden entsteht. Dabei wird man annehmen diirfen, daf die Langskraft lediglich mit w zusammen- 
hangt, nicht aber mit der lotrechten erzwungenen Schwingung, die durch die groBe Biegesteifigkeit 
des Profils Mein gehalten wird. Somit wird man die Kraft wie im vorhergehenden Abschnitt ge- 


* Eine Ableitung dieser Gleichungen aus den Gleichgewichtsbedingungen am Stabelement gibt S. Woi- 
nowsky-Krieger, Ing.-Arch. 13 (1942), S. 197, ferner K. A. Reckling, Ing.-Arch. 20 (1952), S. 137. 

* Diese Bedingung ist identisch mit der Forderung, da das Biegemoment, welches von den eingepragten 
Endmomenten M im Stab hervorgerufen wird, unabhangig von x und gleich M ist. Die Bedingung kann, da 
wir sogleich [Gleichung (13)] die Erregerfrequenz in einen ganz bestimmten Bereich verlegen, in vielen Fallen 
nur unvollkommen erfiillt werden. Wie man auch dann, und sogar bei allgemeinerer Belastung und Lagerung 
als der hier angenommenen, die Form der Gleichung (9) retten kann, erértert Reckling a. a. O. 


fe Ts . = . . . oe . . ‘a ‘ 
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nahert aus dem bewahrten Ausdruck fiir die mittlere Dehnung ableiten und 


I 
BRC, 
Sepa If Cees (10) 


0 


S= 


in (9) einsetzen. Machen wir ferner den Ansatz 


w = v,(t) sin > ; 


9 2 (11) 


bn ee 
sin > » 


v 
P 
der die Randbedingungen erfiillt, so kommt aus (9) und (10) nach kurzer Rechnung das System ge- 
wohnlicher Differentialgleichungen fiir die Zeitfunktionen v, und v, 


Vy + WY + EW, M, v2 cosMt +yorv}—0, | 


op (12) 
; Vy + W2 Vy + € Wy Wy Vy COS Mt =0, f 
worin 
2 IF 
Tt lib) eer : : 
O, = ie oe die tiefste Biegeeigenfrequenz, 
: a 
_ GJa+ ( i ) E Jr 
O2 =F; r die tiefste Torsions- Eigenfrequenz bedeutet und 
P 
M. : IENe : : : 
s=y = mit Mxipp= i) Wy We Ll Up [statisches Kippmoment, aus (12) 
Amp durch 0, = ¥, =0 sofort zu gewinnen], 
F 
Carag 


gesetzt ist. 


b) Ergebnisse der linearen Theorie. Lat man in (12) das Glied mit y weg, so entsteht 
ein lineares Gleichungsystem, fiir das die Eigenschaften der Lésungen eingehend erforscht sind!. 
Im Gegensatz zur einfachen Mathieu-Gleichung kommen jetzt Instabilitaten zweiter Art vor, auf 
die wir unser besonderes Augenmerk richten wollen. Wir beschranken uns dabei wie schon oben 
auf die Diskussion der bedeutsamsten Instabilitat, die jetzt dann eintritt, wenn die Erregerfrequenz 
in die Nahe des Wertes 

Op = Oy Os (13) 


kommt. Um spater den EinfluB des nichtlinearen Gliedes beriicksichtigen zu konnen, miissen wir 
hier kurz wiederholen, wie sich das Stabilitatsverhalten des Systems mathematisch ausdriickt. 

Nach der bekannten Floquetschen Theorie hat (12) ein Fundamentalsystem von Lésungen der 
Form 


vp = eee PR) (fr = 1,2) (14) 
Vi, = eer! py (—s), 
worin wie wblich 
s=wt (15) 


gesetzt ist und 9, die charakteristischen Exponenten, pf, periodische Funktionen der Periode 2 7 
bedeuten. Zunachst soll (14) etwas umgeformt werden. Dazu fiihren wir ys, als komplexe Fourier- 


Reihe 


+0 ; 
Pfr = at Crry et 7 § 
i) —"—co 
ein, setzen 
Oy == u h, 


und nehmen eine Vereinigung der in (14) untereinanderstehenden Partikularlésungen mittels zweier 
Integrationskonstanten ~, und y, vor, indem wir 


1 i aoe 
Vgp == (Vr + Vir) cosy, + + (Vp — Vj) sin y, 


1 FE. Mettler, Ing.-Arch. 17 (1949), S. 418. 
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bilden. So entsteht ein System von vier Partikularlésungen 
-+ (oe) 
vp = orn cos[(h, + n)s+y] (f,r=1,2), (16) 
=— 


die im allgemeinen linear unabhangig sind. 


Die Stabilitatsfrage beantwortet man nun wieder am besten an Hand des zugehérigen Instabili- | 


tatsbereiches (zweiter Art) in der (w, ¢)-Ebene. Abb. 3a paBbt auf ihn ohne Anderung, wenn man 


dort den Wert (13) fiir w) nimmt. Die Grenzlinien haben, wie man mit den bekannten Methoden | 


der linearen Stabilitatstheorie ausrechnet, die Gleichungen 
ee 2 1 
w= Wy Yoo Oy + $5004 (8+ 4G) + (17) 


Liegt der Punkt (w, ¢) auBerhalb des Bereichs, so sind h, und hy, reell, und die Funktionen (16) 
werden durch reelle cos-Reihen gebildet, bleiben also beschrankt. Dasselbe gilt dann fiir die Kipp- 
schwingungen des Stabes, woraus folgt, daB die erzwungene lotrechte Biegeschwingung stabil ist. 


Liegt der Punkt (w, ¢) dagegen innerhalb des Bereichs, so werden h, und h, komplex, so daB die | 


cos-Funktionen in (16) Faktoren der Form e’* mit 0, > 0 abspalten, was offensichtlich Instabili- 
tat bedeutet. Wenn der Punkt (w, ¢) schlieSlich auf eine Grenzlinie des Bereichs fallt, bleibt die 
allgemeine Form der Lésungen (16) erhalten, aber es wird h, = h,, so dafs die Funktionen (16) nicht 
mehr linear unabhangig sind. Zur Vervollstandigung des Fundamentalsystems treten jetzt andere 
Lésungen hinzu, die aber hier nicht interessieren. Die Funktionen (16) wollen wir in diesem Sonder- 
fall ,,Grenzfunktionen“ nennen. Sie entsprechen den Mathieu-Funktionen, welche die lineare Glei- 
chung (6) ,,auf den Grenzen“ lésen. Die besondere Bedeutung der Grenzfunktionen ist nicht so 
augenfallig wie die der Mathieu-Funktionen, die sich durch ihre Periodizitat vor allen andern Lésun- 
gen auszeichnen und deshalb den zwanglosen Ausgangspunkt fiir eine Integration der nichtlinearen 
Gleichung (6) bilden. Wir miissen aber genau so die Integration der nichtlinearen Gleichungen (12) 
an die Grenzfunktionen anschlieBen, deren Bedeutung durch folgende Uberlegung klar wird. 

Wenn man die Wirkung einer beliebig kleinen geschwindigkeitsproportionalen Dampfung be- 
ricksichtigt!, so stellt man zunachst fest, da sich die Grenzlinien des Instabilitatsbereichs ent- 
sprechend wenig verschieben. Im Fall der Stabilitaét, d.h. auBerhalb des Bereichs, bringt die 
Dampfung die Kippschwingungen zum Abklingen, so da8 nur die stabile erzwungene Biegeschwin- 
gung ubrigbleibt. Innerhalb des Bereichs besteht nach wie vor Aufschaukelung der Kippschwin- 
gungen, also Instabilitat. Auf den Grenzen aber 1aBt auch die Dampfung eine beschrankte, weder 
auf- noch abklingende Kippschwingung bestehen. Diese zwischen Abklingen und Aufschaukelung 
stehende, in einem erweiterten Sinne stationdre Bewegung wird durch die — ohne Beriicksichtigung 
der, Dampfung definierten — Grenzfunktionen gut wiedergegeben, und zwar gut aus Stetigkeits- 
griinden, weil man ja in praktischen Fallen stets mit kleiner Dampfung zu rechnen hat. 

Wenn wir nun an die Lésung des nichtlinearen Gleichungssystems (12) gehen, wobei die um- 
standlich zu beriicksichtigende Dampfung nach wie vor weggelassen sei, so kann es sich dabei nur 
um eine bestimmte Aufgabe handeln, namlich die Fortsetzung der Grenzfunktionen v;,(s) im nicht- 


A linearen Gebiet y + 0. In Formeln: die Lésungen vs = v,(s, y) 
von (12), die wir suchen, miissen die Bedingung 
vj(s, 0) = vps) (18) 


erfiillen. Etwas anschaulicher lat sich die Situation und die 
anschliobende Kurven B¥fgabe an dem Amplituden-Frequenz-Diagramm der Be- 


aus der nichtinearen Wegung erklaren. Wir bezeichnen wieder mit A die ,,Ampli- 

Theorie tude“ der Kippschwingung, wobei unter dieser Gréfe, die ja 

bei einer nicht periodischen, raumlichen Bewegung erst defi- 
c niert werden mifte, vorerst irgend ein Maf fiir den gréBten 


Ausschlag aus der Ruhelage verstanden sei. Dann 1a48t sich 


a das genannte Diagramm im linearen Fall, wo die Grenzfunk- 
TE pee tionen nur bis auf einen konstanten Faktor bestimmt sind, durch 
elemente aus der die lotrechten punktierten Linien der Abb. 3c darstellen. 
inearen Theorie 


Aber die lineare Theorie gilt genau genommen nur fiir un- 
Abb. 5, Zum RGN Belienon zur nichtlinearen endlich kleine Amplituden, d.h. sie kann von der Amplituden- 

: kurve A(w) nur Anfangspunkte und die Ausgangsrichtung 
auf der w-Achse, zusammen also Linienelemente, liefern. Die nichtlineare Theorie mu diejenigen 
Kurven A(q) errechnen, die an diese Linienelemente anschlieBen (Abb. 5) 


1G. Grammel, Ing.-Arch. 20 (1952), S. 170. 


Pay 
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c) Berechnung der nichtlinearen Kippschwingungen. Zur Durchfiihrung bedienen 


wir uns der Stérungsrechnung!, und zwar in méglichst engem AnschluB an die Bezeichnungen der 
allgemeinen linearen Stabilitatstheorie?. Dort wurde 


a) Oy Wy 


See a Vs cade a Goh 4 bg (19) 
gesetzt, woraus wegen (13) 
1+a=1 
folgt. Wir filhren ferner neben der unabhangigen Variabeln s eine zweite 
$= {wt+y 3 (20) 


ein. Dann stehen in den Reihen (16) die Funktionen cos (s + 7s). Wir suchen daher die Lésungen 
vy des nichtlinearen Systems (12) in der etwas erweiterten Form 

ve = > Am, cos (ms + ns) 
au berechnen, was den geforderten Ubergang (18) erméglichen wird. Setzen wir (15), (19), (20), in 
(12) ein, so es nach Paes Rechnung 


= ev 
q Sea ess ok ie ve + Oy +ENG%,coss+ygu=0, 
21) 
ev = Ov aay CHS ( 
Coe 20d ae Gees + G3 Vz + € Gi Jo Vy COS S =0, 
Die lineare Stabilitatstheorie (y = 0) setzt nun vy, q und q als Potenzreihen von ¢ an. Wir 


suchen im folgenden entsprechende Potenzreihen von ¢ und y, fiihren aber zur formalen Abkiirzung 
der Rechnung den Kunstgriff ein, von y zunachst den Faktor ¢ abzuspalten und alle Gréfen nach 
Potenzen von € allein zu entwickeln. Die verlangten Potenzen von y kommen dabei, wie sich so- 
gleich zeigen wird, von selbst heraus. Ferner spalten wir von den Lésungen v; einen Faktor A ab, 
unter dem wir den Koeffizienten des Gliedes cos s in der Funktion v, verstehen wollen. Alles in 
allem setzen wir somit 


vp =A Beane ea Deere -) ea) | 


ay 4 i eh. | 
und auBberdem 
; Ae)". (23) 
Fiihrt man (22) und (23) in (21) ein und benutzt die Abktirzung 
av Cv 
Dio) = 23 +2 “sat aa 
so kommt rasch durch Vergleich der Potenzen von ¢ die Folge von Differentialgleichungen fiir die of 
D(vp?) + ae y=0 (f=1,2), (24.0) 
2 (1 oe () tn) Oe cy ep 
D(vf?) + af of) 4.2 (dE +24 +4! ee a ee | 
24.1 
— qi 92 v0) cos s — @? Tv fare yl | eet 
— qx dav") cos s ian y 26 
92,,(0) 2p) 
a ‘ ov (2 2 (1) @ Ove 
D(of?) + geo? + (24° + 4%) ae +2 : Mt gg gg) a 
| (0) pep) pap) 
av (1 i ‘ 
+0 q4? 49%) 2 4 2 2 £24 aS + ISL | ra 


Os Os 
— 919,09) coss—3 GT vo) fir f=l 
— qi Jz V\) cos s fit pay ta, 


1A. Lindstedt, C . R. 47 (1883), S. 1276 und 1353. C. H. Miiller u. G. Prange, Allgemeine Mechanik, Han- 
nover 1923, S. 311 ff. 

2S. FuBnote 1 von S, 359. 

3 Der Index r an p, hat keine Bedeutung mehr und wird weggelassen. 
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usw. Die Lésung von (24.0) ist 
i= COs Sr, 


(25.0) 
D0) c) cos (s — s), 

worin c®) eine Integrationskonstante bedeutet, wahrend cos s verabredungsgemaB den Faktor Eins 

hat. Die Funktionen (25.0) miissen in die Differentialgleichungen (24.1) eingesetzt werden. Eine 

Lésung von (24.1) ohne Sakularglieder ist dann nur mdglich, wenn die kritischen Glieder mit cos s 

in der ersten bzw. cos (s —s) in der zweiten Gleichung (24.1) zum Verschwinden gebracht werden. 

Dieser Gedankengang fiibrt nach einer elementaren Rechnung auf die ,, Hauptgleichungen“ 


= Hee 
2g et 


(26.1) 
und die Lésungen von (24. 1) 
vf) = 4 a cos (2 s—s) + = con35, | 
(25.1) 


v(") = % cos(s + s) + c}) cos (s — s) 


mit einer neuen Integrationskonstanten c/), Eine solche Konstante kann nur in der Funktion vf) 
auftreten, nicht aber in v{!), da hier gemaf unserer Verabredung iiber den Faktor A kein Glied cos s 
mehr vyorkommen darf. Weiterhin liefert Einsetzen von (25.1) in (24.2) und Herausheben der 
kritischen Glieder die Hauptgleichungen 


BONE eo 2 ae 
2 qy qQ + q@) —e tt Onl ==), | 
(26.2) 


— _ — 2 
[2 go (q — 9) + (GQ) — @)] oO) + 2 q, (q — GW) el —E ol) — 0, | 


Die Lésung von (24.2) wollen wir schon nicht mehr anschreiben. Der Fortgang der Rechnung ist 
ja ohne weiteres zu iibersehen. 

Die Hauptgleichungen (26.7) sind Bestimmungsgleichungen fiir die Koeffizienten q) und gq) 
(j = 1, 2,...); ihre sukzessive Aufloésung ergibt 


= 3 
@=2 O44 1, 


| 27.1) 
1 3 ( 
(3) pacpsa yaa 0 ls 
q z(t +40) + pal, | 
5@) = 22) 2 3 @,@ 15 ay 
Ta ei gay 2) ame ee 
(27.2) 


24 


usw. Offensichtlich liegt nun in den Ergebnissen (25.0), (25.1), (27.1), (27.2) usw. insofern eine 
Unbestimmtheit, als die Integrationskonstanten c®), c() usw. zunachst frei geblieben sind. Diese 
Liicke wird aber zwangslaufig dadurch geschlossen, da wir fiir die berechneten Reihen (22) den 
Anschluf (18) an die Grenzfunktionen verlangen. Setzt man in (27.7) I’ = 0, so stehen dort genau 
a ape ais der linearen Stabilitatstheorie. Diese Theorie legt die Konstanten cj) durch die Vor- 
schrift 


(1) il 2 2 3 
2 c q qi 93 (oy r 1 15 r 
g y= A (1 oy) 16 ( (0)? cl 4 392 q2 c(°) ae 39 (qi ae q;) 956 nh i | 


aq‘) : 
70) ga (jen 2) (28) 


fest, und dieselbe Vorschrift miissen wir daher auf die Ausdriicke (27. j ) anwenden. Die Gleichungen 


(28) sind Bestimmungsgleichungen fiir c, c™,...; ihren Wurzeln geben wir wie in der linearen 
Stabilitatstheorie den Index G. Das Ergebnis der Auflésung von (28) ist dann 


¢() Sa /* 9 
© Pes 


il oe) a 
@—la—rssyer. | 


(29) 
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SchlieBlich hat man (29) in (25.0), (25.1), (27.1) und (27.2) einzusetzen und erhalt, wenn man zu- 


gleich (22) und (23) beriicksichtigt, als endgiiltiges Resultat die Anfangsglieder der gesuchten 
Lésungen 


3 


v= Acoss ted h |/h Ceo ag I on 35 ee 
qe 32 : 


v,=+ Aye cos (s — s) Jee [q. cos (s + s) + (1 — 2 q,) cos (s — 5) (30) 
+e By A cos (s—3) +--+, | 


ea 2 Ey eae 3y 2 3 IES) og Bae oe 
ieee ea hte ya 4 +3528 + 49) ge? At FS Vg At: (31) 


om ej 8y. a, 15m eee 
G=HtG nde + 4 + 35% (—2 4) —seey ne 4 ne et. . (32) 


d) Deutung der Ergebnisse. Die Gleichungen (30) geben die Kippschwingungen des Stabes in 
ihrem zeitlichen Verlauf wieder. Dabei ist A als Koeffizient des ersten Gliedes von v, eine Naherung 
fiir die Amplitude der Biegung in w-Richtung, ebenso wie A V41/92 A 
die Amplitude des Torsionswinkels } na&herungsweise angibt. Ge- 
nau genommen sind die Maximalausschlage von w und # ge- 
wisse Funktionen von A; wir nennen jedoch A weiterhin verein- é= const 
fachend die Amplitude. Am wichtigsten fiir die Diskussion der Be- 
wegung ist nun die Beziehung (31) zwischenw,¢ und A. Sie entspricht 
in ihrem Aufbau ganz der schon ausfiihrlich gedeuteten Gleichung (8), 
so da wir auch die an (8) angeschlossene geometrische Darstellung 
in den Abb. 3a bis c tibernehmen und auf ihre schon ausgefiihrte Dis- 
kussion verweisen kénnen. Das véllige Entsprechen der Frequenz- 
Amplituden-Beziehung beim Instabilitatsbereich erster Art und zweiter 


Art laBt es ferner als nahezu sicher erscheinen, da wie im ersten ; ha 
so auch im zweiten Fall die links gelegene Flache in Abb. 3b zu stabi- ae hee 
len, die rechts gelegene zu instabilen Bewegungen gehért. Allerdings __ he Coe x ae 


bleibt dies noch nachzuweisen. Doch muf der Beweis auf eine spatere 1) 6 schematische Darras tae 


Untersuchung verschoben werden. Amplituden-Frequenz-Beziehung bei 
‘ ° * oa ° Variation der Integrationskonstanten 
Stattdessen sei eine erginzende Bemerkung angeftigt. Die zu- 0 1 
(0) und c(1)- 


nachst etwas unanschauliche Bedingung (28) gewinnt man in der 
linearen Stabilitatstheorie aus der Abgrenzung zwischen den stabilen und den instabilen Gebieten 
in der (w, ¢)-Ebene. Im nichtlinearen Fall 1a%t (28) in der Darstellung des (w, A)-Diagramms eine 
andere wichtige Deutung zu, wobei wir der Kiirze halber lediglich den quadratischen Ausdruck 

q —] + eq 4 ¢q®) 
zugrunde legen wollen. Nach (27.1) und (27.2) ist q eine Funktion von c und c), Fiir cl) = cl, 
ce) = c®) wird nach (28) dq/éc) = 0. AuBerdem ist, wie man aus (27.2) sofort sieht, dg/dc = 0. 
Das bedeutet, daB q fiir c) = a), cl) — cl) einen stationdren Wert annimmt. Untersucht man 
dessen Charakter naher, so stellt man fest, daB auf der linken (rechten) Kurve von Abb. 3c, zu der 
die unteren (oberen) Vorzeichen in (29) gehéren, wegen 0°q/dc)” < 0 (>0), @q/dc” = 0 ein Maxi- 
mum (Minimum) von q vorliegt. Dies ist in dem Sinne zu verstehen, daB q bei konstantem A und ¢ 
entweder konstant bleibt, nimlich dann, wenn c©) = c) bleibt und cl) sich allein andert, oder 
kleiner (gréBer) wird, namlich bei Anderung von c). Dieses Ergebnis ist in Abb. 6, einer Wieder- 
holung von Abb. 3c, dargestellt und 14Bt sich folgendermaBen aussprechen: Die Schar von Ampli 
tuden-Frequenz-Kuryen, die man aus (27.1), (27.2) bei Variation der Integrationskonstanten c() 
und c@) bekommt, hat zwei durch die Bedingung (18) bzw. (28) ausgesonderte Grenzkurven (31), 
in deren (in Abb. 6 schraffierten) Zwischenraum keine Kurven fallen. Man kann auch sagen, dab 
auf dem linken (stabilen) Ast die Amplitude bei fester Erregerfrequenz @ ein Minimum hat. 

Nach der Einleitung besteht eine Aufgabe dieser Arbeit in der Schilderung der sichtbaren Er- 
scheinung der berechneten nichtlinearen Schwingungen. Die nun ausfiihrlich diskutierte Amplitude 
ist ein Teil dieser Erscheinung, Das Bild wird vervollstandigt durch eine Feststellung, die wohl 
das Merkwiirdigste an dem Verlauf der Kippschwingungen darstellt. Wir beschranken uns dabei 
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bewuBt auf das grob sichtbare Phanomen ohne numerische Feinheiten. Auf dieses stoBt man bei der 
Uberlegung, da ¢ praktisch stets eine sehr kleine Zahl sein wird und daB auch der zweite Parameter, 
dessen Potenzen in den Reihen (30) bis (32) stehen, nimlich y A’, aus Konvergenzgriinden klein sein 
muB. Ubrigens ist y A? = (A/2 i,,)?, wenn man mit i, den Tragheitsradius des Querschnitts in 
der Waagerechten bezeichnet, woraus folgt, daB die Amplitude héchstens die GréBenordnung des 
Tragheitsradius haben darf. Wenn wir nun aus jeder der Reihen (30) bis (32) das numerisch tiber- 
wiegende, namlich das erste Glied nehmen, kommen wir zu folgender Aussage. Die Erregerfrequenz 
w liegt bei oder besser gesagt etwas tiber der Schwebungsfrequenz w, + W,. Die Variable s geht 
nach (20) und (32) inw,t + y, die Variable s — s nach (15) und (13) inw, t —y tiber. Somit kann 
man in der Gesamtbewegung drei Komponenten unterscheiden. Die Biegeschwingung in der Lot- 
rechten, die hier die Rolle der erzwungenen Grundbewegung spielt, verlauft mit der Frequenz 
(1 + Wy des erregenden Moments. Als Nachbarbewegung kommt hinzu einmal die Biegeschwingung 
in der Waagerechten (wiedergegeben durch v,) im Takt der zugehérigen Biege-Eigenfrequenz w, 
und daneben die Torsionsschwingung (bestimmt durch v,), die entsprechend mit der Torsions- 
Eigenfrequenz jw, verlauft. Da die drei Frequenzen w, + W.,, W, und ~, im allgemeinen nicht kom- 
mensurabel sind, ruft somit das harmonische Erregermoment eine im ganzen nicht periodische Be- 
wegung des elastischen Kérpers hervor. 

Das ist ein Ergebnis, das der klassischen Elastokinetik véllig fremd ist. Es ist nicht nur theore- 
tisch von Bedeutung, sondern auch fiir den Praktiker von Interesse. Man denke etwa daran, da} 
man in der Technik haufig Schwingungsmessungen an Tragwerken ausfiihrt und aus den festge- 
stellten Befunden auf die erregenden Krafte Rickschliisse zieht. Dabei arbeitet man wohl stets 
mit dem Grundsatz der klassischen Schwingungslehre, daB erregende Kraft und erzwungene 
Schwingung dieselbe Frequenz haben. Wiirde man auf dieser Grundlage eine Balkenschwingung 
der oben diskutierten Art untersuchen, so mite man auf Erregerkrafte der Frequenzen w, und w, 
schlieBen, die in Wirklichkeit gar nicht vorhanden sind. 


Aus dem Institut fiir Mechanische Schwingungstechnik der Techn. Hochschule Karlsruhe. 
(Eingegangen am 1. Dezember 1954.) 
Anschrift des Verfassers: Professor Dr. E. Mettler, Karlsruhe-Durlach, Geigersbergstr. 12. 
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Uber eine mogliche Verbesserung der Ritzschen Methode 
Von Th. Péschl 


x Einleitung. Bei der Ritzschen Methode zur Lésung von Randwertaufgaben von Differential- 
gleichungen wird eine Folge von Naherungsfunktionen w;(x) (i =1,2,...n) ausgewihlt, die 
einzeln die Randbedingungen erfiillen, aber nicht der gegebenen Differentialgleichung geniigen. 
Diese werden mit zunachst unbekannten Koeffizienten a; multipliziert, und die Summe 


n 
w = 3) a, w;(x) 


i=l 
wird fiir die unbekannte Funktion w = w(x) in das Integral 


A= f F(w"’, w’, w, x) dx = Min. 
i 


eingesetzt, von dem die gegebene Differentialgleichung nach den Regeln der Variationsrechnung 
abgeleitet werden kann. Diese Differentialgleichung hat, wenn F die zweiten Ableitungen der 
unbekannten Funktion enthalt, die Form 
Gh) a Gl d (oF oF 
Luss dae (ea) dx (srr) ) 


T Ow ‘ 


wobei w’ = dw/dx usw. ist. 


Nach Ausfiihrung der Integrationen wird A eine Funktion der a;: 
A= A(a,; G35" «= Gn) 


und die unbekannten Koeffizienten werden nach den Regeln tiber gewéhnliche Maxima und Minima 
durch die Gleichungen erhalten: 


—_- = 0 A ad Nberameert ul 8 


Wenn die Funktion A eine quadratische Funktion von w und den Ableitungen ist, so ist dies ein 
System von n linearen Gleichungen, deren Auflésung keine grundsatzlichen Schwierigkeiten be- 
reitet. Bei statischen Aufgaben sind dies nicht-homogene Gleichungen, bei Schwingungs- oder 
Stabilitatsaufgaben sind die Gleichungen linear und homogen, und ihre Wurzeln stellen Naherungen 
fur die Eigenwerte dar. 

Bei den Anwendungen dieser Methode, wie sie gewéhnlich ausgefiihrt werden, ist folgender Um- 
stand bemerkenswert: Die Naherungsfunktionen w; werden nur gemaB den Randbedingungen 
(gegf. auch den Symmetriebedingungen u. dgl.) ausgewahlt, aber ohne jede Riicksichtnahme auf 
die besondere Form der gegebenen Differentialgleichung, deren Lésung gesucht wird. Wenn man 
also, wie dies in vielen praktischen Fallen zutrifft, die Lésung der Differentialgleichung fiir eine 
Folge von Zahlenwerten eines (oder mehrerer) Parameter zu ermitteln hat, die in der Differential- 
gleichung auftreten, so gehen diese Parameter nur in die Werte der Koeffizienten a; ein, mit denen 
die Naherungsfunktionen w; multipliziert werden, um durch die Summe dieser Produkte die Lésung 
aufzubauen. Diese Naherungsfunktionen selbst sind aber von den in der gegebenen Differential- 
gleichung auftretenden Parametern véllig unabhingig. Die Eigenart der gegebenen Differential- 
gleichung ist auf die Wahl des Funktionensystems der w; nur insoweit von EinfluB, als die Randbe- 
dingungen durch den AbleitungsprozéB der Euler-Lagrangeschen Gleichungen durch die bei der 
Produktintegration entstehenden Glieder als ,,natiirliche Randbedingungen“ mit dem zugehérigen 
Variationsproblem verbunden sind. 

Durch diesen Aufbau der gesuchten Lésung wird es verstandlich, das fiir einen groBeren Bereich 
der Parameterwerte im Einzelnen gréRere Abweichungen von den richtigen Lésungen auftreten 
kénnen und auch mehrfach festgestellt wurden. 

Fiir diese Methode sind von verschiedenen Forschern Verbesserungen vorgeschlagen worden, die 
eine Vereinfachung der Rechenarbeit mit sich bringen, so von B. G. Galerkin, C. B. Biezeno, 
R.Grammel, H. Hencky u.a. Eine etwas anders geartete Abanderung der Methode soll im folgenden 
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dargelegt werden, u. zw. beruht diese auf einer Verbindung der Ritzschen mit der Methode der 
Kollokation, die sich gerade in dieser Verbindung als besonders vorteilhaft erweist. Durch Ver- 
gleich mit den von L. Collatz in seinem bekannten Buche Numerische Behandlung von Differential- 
gleichungen, 2. Auflage 1955 (im folgenden mit L. C. zitiert) angefiihrten Methoden konnte fest- 
gestellt werden, daB bei einigen der durchgerechneten Beispiele mit Hilfe dieser Abanderung 
durch ein einziges Polynom fast dieselbe Genauigkeit erreicht werden kann wie sonst fiir zwei, 
manchmal auch fiir drei Naherungspolynome. Auch fiir die Galerkinsche Methode (u. a.) 1aBt sich 
diese Abanderung verwerten. 


2. Die Methode. Den Ausgangspunkt der folgenden Betrachtungen bildete die Frage, wie man 
es erreichen kann, da auch die Naherungsfunktionen w; selbst — und nicht nur die Koeffizienten 
a; — von den Parameterwerten abhangig werden, mit dem Ziel, auf diese Weise die Methode selbst 
zu vereinfachen und woméglich zu verbessern. Man wird also trachten, die gegebene Differential- 
gleichung selbst fiir die Wahl der Naherungsfunktionen zu verwerten. Dies geht nun freilich nicht 
fiir den ganzen Verlauf der gesuchten Lésung; es ist aber méglich, die Differentialgleichung fir 
einen Punkt (gegf. auch fiir mehrere Punkte) des in Betracht kommenden Intervalls durch die 
Naherungsfunktionen zu _ befriedigen. 

Es wird also gefordert, daB die gegebene Differentialgleichung fiir eine geeignet gewahlte Stelle x, 
des betrachteten Bereichs (oder deren mehrere) erfiillt sein soll, und zwar unter Verwendung der 
Naherungsfunktionen w;. Wir schreiben also die Bedingungen vor 

VA sal (Ga 1,2.65 

und betrachten diese als zusatzliche, zu den Randbedingungen hinzutretende Bedingungen, die den 
Naherungsfunktionen auferlegt werden. Wenn als Naherungsfunktionen Polynome gewahlt werden, 
so kénnen diese zusatzlichen Bedingungen dadurch erfiillt werden, da der Grad der Polynome durch 
weitere Glieder erhéht wird. Die Konstanten dieser Polynome werden dann nicht nur durch die 
Randbedingungen, sondern in gleicher Art auch durch diese zusatzlichen ,,Kollokationsbedingun- 
gen bestimmt. Auf diese Weise ist es méglich, die besondere Form der gegebenen Differential- 
gleichung mit den Randbedingungen in Verbindung zu bringen und zugleich zu erreichen, daB die 
in der Differentialgleichung enthaltenen charakteristischen Parameter auch in die Naherungs- 
funktionen w;(x) selbst — und nicht nur in die Beiwerte a; — eintreten. 

Selbstverstandlich ist dieses Verfahren auch in der Weise anwendbar, da} nicht nur die gegebene 
Differentialgleichung selbst, sondern auch deren Ableitungen nach x fir solche ,,Kollokations- 
bedingungen“, die also von der Form sind 


d d? 
Es (Phed}| =, E (Eee}}] = ‘aes 

herangezogen werden, wobei jedes Mal der Grad der Naherungspolynome um eine weitere Einheit 
erhéht wird. — Auch die Erweiterung auf mehrere unbekannte Funktionen und auf partielle 
Differentialgleichungen ist grundsatzlich auf dieselbe Weise durchfihrbar. 

Es mége nochmals bemerkt werden, da die hier angegebene Erweiterung des Ritzschen Ver- 
fahrens auch von B.G. Galerkin u. a. herriihrende Art der Konstantenbestimmung in sonst unver- 
anderter Form anwendbar ist. 


3. Beispiele. Es diirfte zweckmafBig sein, den geschilderten Vorgang an Hand einiger einfacher 
Beispiele zu erlautern, die so gewahlt werden mégen, daB ein Vergleich mit den nach anderen Metho- 
den erhaltenen Lésungen durchfiihrbar ist. Wir beschrinken uns hier auf gewohnliche Differential- 
gleichungen und wiihlen drei Beispiele aus L. C. und ein solches iiber die Verformung der Wand 
eines zylindrischen Behalters von konstanter Wandstarke, die durch Wasserdruck belastet ist. Wir 
k6énnen uns dabei jedes Mal auf ein einziges Naherungspolynom beschranken und wollen prifen, 
wie gro die Abweichungen gegeniiber den bekannten Lésungen sind, die durch zwei bzw. drei 
Naherungsfunktionen erhalten worden sind. 


a) Gegeben sei die Differentialgleichung (Siehe L. C. Seite 195/6): 


Mie) yO (ees reat 
mit dem zugehérigen Variationsproblem 


ie 
Jiyl=_f y?—(+#) y°—2y] dx = Min, 
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Die erste Naherungsfunktion, die den Randbedingungen entspricht, ware A (1 — x"). Die Hinzu- 
nahme einer Kollokationsbedingung verlangt die Erhéhung des Grades des Polynoms; wir setzen 


pee Oe Ses 


Als Kollokationspunkt liegt es nahe die Stelle « = 0 zu wahlen. Die Differentialgleichung gibt 
dann fiir diesen Punkt 2 C + 4A + 1 = 0, somit ist 


1 


Die vorgeschriebenen Randbedingungen liefern dann (fiir « = + 1) 


1 1 
Q==A sul + A)+D, also D =~ (1— A). 
Das erste Naherungspolynom lautet mithin 


y =A—5 (14 4) 2 et | 


2 2 2 
Setzt man dieses in den Ausdruck fiir J[y] ein, so liefert die Bedingung fiir das Extremum 
oJ 
0A = 0 
nach Ausfiihrung der elementaren Integrationen fiir A den Wert 
A= 0.927. 6 


mit einer Abweichung von etwa 1/2°% gegeniiber dem mit zwei bzw. drei Naherungsfunktionen 
von L. C. erhaltenen Wert 0,933... Das erste Naherungspolynom erhalt demnach die Form: 


y = 0,927 (1 — 0,987 x? + 0,0543 x*). 
b) Gegeben sei die Eigenwertaufgabe (L. C. Seite 196/7): 
SpE syle A x) ye yO) yl) = 0 
mit dem zugehérigen Variationsproblem 
1 
JUy] =f +) ¥2—A (1 + 2) 9°] dx = Min. 
Wir machen den Ansatz, durch den die Randbedingung y’(0) = 0 sogleich erfillt ist, 
y= A+ Os + Dx. 
Die Kollokationsbedingung fiir x = 0 lautet 


By Cd Palio ist C= A. 
Ferner gibt die zweite Randbedingung y(1) = 0 
Mee CRD 0. woraus D=—(1—) 4, 


so daf das erste Naherungspolynom die Form erhalt 
¥ =4)l—> | x : 
Man beachte, daB der Eigenwert A auch im Naherungspolynom vorkommt. (Ahnliches gilt auch fiir 
das in a) gegebene Beispiel.) ; 
Der Ausdruck J[y] liefert sodann nach Ausfiihrung der elementaren Integrationen fiir 2 die 
kubische Gleichung 
13 3 AN eee 


3 2 = ORY See 
336 7 Sy ee ¢ 


Der erste Eigenwert ist nach L. C. /, ~ 10/3. Wir setzen A = 10/3 + € und erhalten 


Ax“ +0,026... 


Die Abweichung gegen den von L. C. erhaltenen Wert 10/3 betragt daher weniger als 1%. 
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Die zugehérige Eigenfunktion ist wie bei L. C. 
y =4(1—3# ie — a (0,6 — x? + 0,4 x4). 


Die verbesserte erste Naherungsfunktion leistet daher ungefahr so viel wie die zwei oder drei sonst 
verwendeten Naherungsfunktionen. 

c) Formanderung einer zylindrischen Behalterwand nach Abb. 1. von der Héhe h und dem 
mittleren Halbmesserr. Ferner sei w = w(x) die 
Ausweichung der Wand an der Stelle x zufolge 
der Belastung durch den Wasserdruck, 6 = (x) 
die mit x im allgemeinen veradnderliche Wand- 
starke, dann ist das Biegemoment an der 
Stelle x 


v 


d*w 
ee eae 


und die Querkraft 


Die Differentialgleichung des Problems ergibt sich 
durch die Euler-Lagrangesche Gleichung des 
Variationsproblems 


Abb. 1. Verformung eines zylindrischen Wasserbehilters. 


d?w 


de 


1 
etd Wa We [ (3 DEES kdw—1Ew) Mie 
0 


Darin ist zur Abkiirzung gesetzt 


x 12h — 2%) 1285 (1 — 9) y 
eae a or r? s he Bec’ 


Die zugehérige Euler-Lagrangsche Gleichung ist dann 


a d*w 
pik Sea ey f= 
L[w] = de (0 7) kdbw—lé=0 
mit den Grenzbedingungen: 
f= 1: w=0, w’ =0 (eingespanntes Ende) , 


E=0: M,=0, Q,=0 (freies Ende) 
oder 
OF == 05 (0% ann ae Oe 
Fir eine Behalterwand mit konstanter Wandstarke §6 = 6, = konst. setzen wir noch 


=a 


k l w 
j2 


=H, 53 A> ay, 


Das Variationsproblem nimmt dann die einfachere Form an 
1 
1 
A= | iE W'? + xW2—é W) dé = Min. 
0 


und gibt als Euler-Lagrangesche Gleichung 


Lilw] = WY + 4xW—E=0 
mit den Randbedingungen 
é=1: Warr WwW’ =0, 
E—0: Wo = 08 Wii) 
Ein Polynom als Naherungsfunktion mu also nach dem iiblichen Verfahren entsprechend 
diesen Randbedingungen vom vierten Grade sein und die Form haben 


Wm all—se tae). 
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Geht man mit diesem Ansatz in den Ausdruck fiir A, so ergibt sich nach Ausfiihrung der elemen- 
taren Integrationen: . 


ees 
~ 1296 + 104% ° 


Bei diesem Verfahren geht der in der gegebenen Differentialgleichung enthaltene Parameter x 
nur in die Koeffizienten a, nicht aber in die Naherungsfunktion W selbst ein. Wenn man jedoch die 
Beziehung zwischen W1V und W, wie sie in der gegebenen Differentialgleichung vorliegt, auch noch 
fiir irgend einen geeignet gewahlten Punkt € heranzieht, als welcher sich im vorliegenden Fall der 
Punkt € = 0 empfiehlt, so erreicht man, daB die Form der gegebenen Differentialgleichung auch bei 
der Festlegung der Naherungsfunktionen beriicksichtigt sind, und da® x auch in diesen selbst 
auftritt. 

Bei dem vorliegenden Beispiel hat man demnach zu setzen 


w =a+bét cH+ d&, 


oe b +4c8+4+ 5d#,z 
eo = 2c + 20d&, 
ee 24cé + 60d & 


und als zusatzliche Kollokationsbedingung 


[WV 1% W—é] 7 te na — 0 


o— 


einzuftihren. Dadurch erhalt man 


eee 
ewe 
und durch Auflésung der verbleibenden Gleichungen fiir das erste Naherungspolynom 


I fee x I (2 
; eral) yes ag e +1)e]. 

Durch die angegebene Erweiterung wird, wie ersichtlich, erreicht, da der Parameter x auch in 
die Naherungsfunktion w; selbst eingeht. 

Fuhrt man jetzt diesen Ausdruck fur W; in das Integral A ein, so kann man die Integrationen 
mit elementaren Mitteln ausfiihren und erhalt fiir A einen quadratischen Ausdruck; aus diesem 
folgt durch Ableitung nach a eine lineare Gleichung, aus der sich a berechnen laBt. Es ergibt sich 

x + 120 
= 0, (2)89 2 - 3,432 22 + 685 x 4+ 3000° 

Aus der Gleichung fiir W(é) folgt u. a., da sich das Vorzeichen von W’(é) am freien Ende fiir 
x = 120 umkehrt; wahrend es fiir x < 120 negativ ist, wird es fir x > 120 positiv, was auch die 
exakte Lésung ergibt. Der arspriinglich angegebene Ansatz gibt fiir alle Werte von x das gleiche 
negative Vorzeichen der Ableitung, was offenbar irrig ist. 

Dieses Verfahren der Verbesserung kann man beliebig fortsetzen, indem man als weitere Bedin- 
gung die Ableitung der gegebenen Differentialgleichung nach & fiir den Punkt § = 0 hinzufiigt, 
also die Gleichung 


W(t) =a 


a 


d 

—, {wv Ww — =O) 

ent ae é} De 
und gleichzeitig den Grad des Naherungspolygons um eine weitere Einheit erhéht; man hatte also 
dann zu setzen 

Wear bf 4 ¢t* od £1 f° 

und erhielte durch Einsetzen in die vier Randbedingungen (wovon zwei bereits beriicksichtigt sind) 
und der zwei zusatzlichen Bedingungen vier lineare Gleichungen zur Bestimmung der vier Beiwerte 
b,..,ein Abhangigkeit von a, wobei auBerdem auch x auftritt. Durch Kinsetzen des so erhaltenen 
Ausdrucks in das Minimalproblem A = Min und nach Ausfiihrung der elementaren Integrationen 
ergibt sich dann durch die Bedingung 


oA 
Fn ae 


wieder eine lineare Gleichung fiir a und damit W = W (&, x) als erstes Naherungspolynom. 
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d) Als letztes Beispiel behandeln wir die nicht-lineare Randwertaufgabe (L. C. Seite 197/8) 
vr 3 
Ye Ps, (O) Ay (2) alee 


Bei der a.a.O. angegebenen Lésung wird deutlich, da auch bei Erfillung der Randbedingungen 
die erhaltenen Lésungskurven der gegebenen Diffentialgleichung — auch bei einem zweigliedrigen 
Lésungsansatz — nur sehr unvollkommen geniigen. Es zeigen sich vielmehr vor allem in den End- 
punkten 0 und 1 erhebliche Abweichungen zwischen den Werten von y” und 3 y?/2. Auch hier er- 
weist sich die Verbindung der Ritzschen Methode mit der der Kollokation von Vorteil. 

Wir beschrinken uns darauf, den Ansatz und die Durchfiihrung der einzelnen Schritte des Ver- 
fahrens anzudeuten, wobei auch hier bemerkt werden soll, daB sich dieses in mannigfacher Weise vari- 
ierenlaBt, so daf fiir die Anpassung an ein vorgelegtes Problem im allgemeinen verschiedene Méglich- 
keiten in Betracht kommen. In den meisten Fallen wird man als Naherungsfunktionen Polynome 
wihlen, wobei die Anzahl der Glieder von der Zahl der zu erfiillenden Bedingungen abhangt, 
wahrend der Grad der einzelnen Glieder in weitem Mafe willkirlich bleibt und grofe Mannig- 
faltigkeiten zula®t. Im Ganzen gehen wir darauf aus, die Lésung durch ein einziges Polynom 
darzustellen und bemerken nur, da die Einfiihrung mehrerer derartiger Polynome ohne Weiteres 
méglich ist, aber zu einer auBerordentlichen Vermehrung der Rechenarbeit fiihren wiirde. 

Fiir die vorliegende Aufgabe ist die Einfiihrung eines fiinfgliedrigen Polynoms erforderlich, das 
in der einfachsten Form gewahlt werden mége: 


y=AtBx4+C?+Dx+ Ex. 
Zur Festlegung der fiinf Konstanten A,. . . ., E dienen die zwei Randbedingungen und die zwei 


Kollokationsbedingungen an den Randpunkten, wahrend die fiinfte Konstante durch die Ritzsche 
Methode bestimmt wird. 


Durch die Randbedingungen erhalt man 
fur 40 Ae Al se dl. 
fiir Roa dase ==) B+C+D+E=-—3, 


und zufolge der Kollokationsbedingungen: 


fiir oe Oly a Ae aus yl aay GaN 
fir eee ledy eel 2C+6D+12E=-+. 

Daraus ergeben sich die Konstanten 
D =— 26,25 —2 B und 9b =.11,25 -— B 


und der Ansatz fiir die Naherungsfunktion 
y = (44 12 x? — 26,25 x8 4 11,25 x4) + B(x—2 x8 + Soe 


Man beachte, daB dieser Ansatz fiir alle Werte von B die Rand- und die Kollokationsbe- 
dingungen erfiillt. Die Festlegung von B kann nun entweder nach der Ritzschen Methode oder 
durch Einfiihrung einer weiteren Kollokation (etwa fiir x = 1/2) erfolgen. Hierzu fiihre man diesen 
Ansatz in das der gegebenen Differentialgleichung zugeordnete Variationsintegral 


1 
J = f (y2+9%) dx = Min. 
0 


ein und erhalt gemaB dem Ritzschen Verfahren 
1 


i 
op = 2 [fy Go dx +3 [yo dx <0. 
0 0 
Dadurch erhalt man fiir B eine quadratische Gleichung, in der jedoch der Koeffizient von B2 (bei 
Anwendung des Rechenschiebers) sehr nahe verschwindet. Man erhalt fiir B den Naherungswert 
B =— 17,65 
und damit den eingliedrigen Lésungsansatz der Aufgabe 
y =4—7,65 x + 12 x2 — 10,95 x3 + 3,60 x4, 


der schon recht gute numerische Ergebnisse liefert. 
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In dieser Gleichung kann man aber noch eine Korrektion dadurch vornehmen, da®B man durch 


den Ansatz 
B=— 7,65 + 8 
nachtraglich noch eine weitere Kollokationsbedingung etwa fiir den Punkt « — 1/2 einfiihrt. Ge- 
ma der angegebenen Gleichung erhalt man durch eine genauere Rechnung fiir B die Werte 
B, =— 8,622 und B, = — 43,094 
folgen. Damit ergeben sich fiir die beiden Lisungen der Differentialgleichung die Polynome 
yy = 4— 8,622 x + 12 x? — 9,006 x3 + 2,628 x4 


und 


Vo = 4— 43,094 x + 12 x2 + 59,938 x3 — 318,44 x4 
und als Funktionswerte fiir die erste Lésung 


Yea 4 e591 — 9,045 = 6.341 1 
ys —43,094 | —27,847 —12,063 22,309 +33,344 
fie 24 90,024 108,282 18,774 1,5 

3 

ek. 24 40,87 122,718 60,726 1,5 


Abb. 2. Verlauf der Lésungen der nicht-linearen Randwertaufgabe nach Beispiel d). 


Die so erhaltenen Lésungskurven sind in Abb. 2 dargestellt. Man beachte, daf fiir y = 0 gleich- 
zeitig auch y”’ = 0 und y’” = 0 ist, so daB die zweite Lésung in y = 0 keinen gewohnlichen Wende- 
punkt sondern einen Flachpunkt besitzt und durch diese Punkte in nach oben konkaver Kriim- 
mung durchlauft. 


25 
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4, Schlu&bemerkung. Durch die iibliche Form der Ritzschen Methode ergeben sich Naherungs- 
funktionen, die wohl die vorgegebenen Randbedingungen erfiillen, aber in vielen Fallen der vor- 
gelegten Differentialgleichung nur in unvollkommener Weise geniigen. Wenn man jedoch die 
Naherungsfunktionen so umgestaltet, daB sie auver den Randbedingungen auch noch Kollokations- 
bedingungen fiir einen oder mehrere Punkte des vorgegebenen Bereiches erfillen, so gelingt es oft 
schon mit einer so veranderten Naherungsfunktion eine Verbesserung der Werte der gesuchten 
Funktionen und ihrer Ableitungen im ganzen Bereiche zu erhalten. Die Methode lat sich in 
mannigfacher Weise ausgestalten und ist fiir viele Randwertaufgaben fiir Differentialgleichungen, 
die aus einem Variationsproblem entspringen und auch fiir mehrere unabhangig und abhangig 
Verinderliche anwendbar. Es ist zu vermuten, da® durch die zusatzlicheén Kollokationsbedin- 
gungen der Wert der Minima der Integrale herabgesetzt und dadurch die Annaherungen ver- 
bessert werden. 

Zur vollstandigen Lésung des Problems wiirden sodann noch die Erledigung der Konvergenz- 
frage und auch eine Fehlerabschatzung gehéren. 


(Eingegangen am 20. Januar 1955). 
Anschrift des Verfassers: Prof. Dr. Theodor Péschl, Karlsruhe i. B. Wendtstr. 5 Fernruf 52 480. 
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Reine ebene Biegung bei einem nichtlinearen Elastizititsgesetz 


Von F. Jindra 


1. Einleitung. In einer friiheren Mitteilung! wurde der ebene Spannungszustand bei einem 
nichtlinearen Elastizitatsgesetz fiir kleine Verzerrungen durch Einfithrung der Airyschen Span- 
nungsfunktion behandelt. In kartesischen Koordinaten x, y lassen sich die Gleichgewichtshe- 
dingungen des ebenen Spannungszustandes identisch erfiillen, wenn man die Spannungen von 
einer Spannungsfunktion F(x, y) in bekannter Weise ableitet: 

(od Or er 
oe Go? 5 aa Sy ony: (1) 
Mit dem Kompressionsmodul K und dem Schubmodul G geniigt die Spannungsfunktion der nicht- 
linearen Differentialgleichung 


ety, Oe __ 1 [@F (th) | OF M(t) @F ae(e8)] _ 
_ | TK PS) + 3G Blo) 4r| Re || ayers" oat crys Oney Os ay | Z (?) 
mit den Abkiirzungen 
1 1 /@F | @F 
ae Pater joer 
a 9G? M4") he \ ox oy Ox? dy? | * 


Die sogenannte Kompressionsfunktion k(s,) und Schubfunktion g(t) beschreiben das nichtlineare 
Verhalten des Materials und werden als Potenzreihen der dimensionslosen invarianten Spannungs- 
groBen (3) eingefiihrt. Dabei beschrankt man sich im wesentlichen wie in der ersten Mitteilung auf 
den Sonderfall, daf diese Funktionen in der Form 
Ks) = 1, g(t) =1+ 82% (g. = konst.) (4) 
dargestellt werden kénnen. 
Die vorliegende Mitteilung beschaftigt sich mit den Anwendungen der Methode der Airyschen 
Spannungsfunktion auf das ebene Problem der reinen Biegung. Wir untersuchen die Biegung eines 
Rechteckstreifens, die leicht in geschlossener Form gelést werden kann, und die eines Kreisring- 
sektors, die sich aber nur durch ein Naherungsverfahren behandeln ]48t. Die numerische Aus- 
wertung der Ergebnisse ergibt wesentliche Anderungen der Spannungsverteilung gegeniiber der 
linearen Theorie und einen erheblichen Abbau der 4 
Spannungsspitzen durch kleine Abweichungen des Elasti- 


zitatsgesetzes von der Linearitat. 

2. Der Rechteckstreifen. Im folgenden betrachten 2 | S 
wir einen Rechteckstreifen von der Héheh und be- = 
liebiger Lange, der an beiden Enden durch Biege- M 
momente M beansprucht wird (Abb. 1). Die als klein 


anzusehende Dicke des Streifens senkrecht zur Zeichen- 
ebene ist gleich Eins gesetzt. 

Wir gehen hier ahnlich wie in der linearen Theorie von der Voraussetzung aus, dafi der Span- 
nungszustand der reinen Biegung durch die einzige Spannungskomponente o,(y) beschrieben 
wird, also 


Abb. 1. Durch Momente belasteter Rechteckstreifen. 


er , Fr er 
On = Fa = oa(y) » Ee Bs ? YP ee oy 8 : (5) 
Nach (3) wird 
1 aot Nee 
OK aes 0 0G? 
Dann vereinfacht sich (2) unter Beachtung von (4) zu 
@ ee 
a (o. A o:) = 0 (6) 


1 F, Jindra, Ing.-Arch. 22 (1954), S. 121. 
20% 


o74 Jindra: Reine ebene Biegung bei einem nichtlinearen Elastizitatsgesetz. Ingenieur-Archiv 


mit der Stoffzahl 


= gk a 
4=GKLOG” (7) 

die hier allein die Abweichung von der Linearitat kennzeichnet. Die Integration ergibt sofort 
on +fAoi =Cy¥+D, (8) 


wobei C und D Integrationskonstanten bedeuten. Wir setzen weiter 
D=0, (8a) 


Der Fall D +0 entspricht einer zusatzlichen gleichmaSigen Beanspruchung in der x-Richtung 
(welcher Fall uns hier aber nicht interessiert). 

Mit Hilfe der Cardanischen Formel lésen wir die kubische Gleichung (8) auf. Als einzige reelle 
Wurzel finden wir die Biegespannung 


om 3|/ (ley rVay tie + Cy—1Or 428 | ) 


mit der Abkiirzung yu = a ; 


Fiir die Randspannung an der Stelle y = h/2 folgt 


(o.)y=n2 = 2 Yn (Cn + VOR Ean + Ch—yOm +p). (10) 


Die zugehérige Spannungsfunktion F wird dann 


ae Se ace a). iff Ape 
P= [ford =x0 7/4 (19 cry? 16 x —9 Cy CRF Bp) Cy + Ore + 2p | 


+ (19 C2 y?—16 pp + 9 Cy (Oy F 2p) Cy — (ey 2 u|.| 


Die Konstante C errechnet sich aus der Forderung, daB das Moment der in einem beliebigen Quer- 
schnitt wirkenden Spannungen o, das vorgeschriebene Biegemoment M ergeben muB: 


M =f oydy =a ie io (Cy? +p) — Cy (@x 4 2p] [Cy +/+ 2p 
+[6(C2y? +p) + Cy (COP 2 yp] /Cy—[Cy +2 u oe 
Nach Einsetzen der Integrationsgrenzen erhalten wir 
oes relia (cm +4) —ChY@h + 8p] /Ch + Ch + 8p 
he 140 C?R (12) 
+[6(CR +4) + ChYCR + Bu] Jon— yore + Buf 
und nach einiger Umformung die Bestimmungsgleichung fiir C 
10125 A? h? C> — 4900 A? M (2802 M2 + 27 h4) C4 | (13) 


+ 9576 A h> C> — 117600 A h? M C2 + 2240 h3 C— 26880 M =0. |{ 


Fir 7 =0 gehen unsere Gleichungen in die bekannte Lésung der linearen Theorie iiber. 
Aus (8) und (13) finden wir 


1 
Onn ys M = 7, C* h’ (14) 
und somit die Randspannung 
M 
(0%)y =n/2 = 675 . (15) 


Hierin bezeichnet C* die Integrationskonstante der linearen Theorie. 
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Abb. 2 zeigt Se Verteilung der Spannung o, lings der y-Achse fiir einen beliebigen Querschnitt, 
und zwar fiir M/h ie 200 kg/cm und A = 5.10-§ cm4/kg?. Gestrichelt eingetragen ist die Lésung (14) 
der linearen Theorie. Die Biegespannung o* ist iiber dem Querschnitt geradlinig verteilt, wie es die 


16, 


t (x) Yi 5 
1200 


Kg/ené « /| 
1000 L ce Kyles 


7000 


800 


600+ 


400 


200 


M/n? 
0 50 100 750 ‘kg/cm 200 


Abb. 2. Verlauf der Biegespannung. Abb. 3. Randspannung im Rechteckstreifen. 


elementare Biegelehre voraussetzt. Man erkennt einen starken Abbau der Spannungsspitze an den 
Randern y = + h/2 infolge der Nichtlinearitaét. Die Werte der Randspannung (ox)y—4/2 fiir 
verschiedene Annahmen von der Stoffzahl J in Abhangigkeit von der GréBe M/h? zeigt Abb. 3. Um 


AOz den Einflu® der Nichtlinearitat besser zu uber- 
- 085 Or sehen, ist schlieBlich in Abb. 4 die Anderung 
der Spannung Jo, =o,—o; an der Stelle 
y =+h/2 im Verhaltnis zur Spannung o% 
der linearen Theorie nach (15) abhangig von 
-~9.20 M/h? dargestellt. 
c 3. Der Kreisringsektor. Als zweites Bei- 
ne Ry spiel behandeln wir den in Abb. 5 dargestell- 
Zw ten krummen diinnen Stab von der Gestalt 
~ G15 iP eines Kreisringsektors mit den radialen Be- 
4 
0,10 nO 
ne 
ME 
-0,05 i + 
jet 
oa Mn? 
0 50 100 150 kg cm 200 


Abb. 5. Durch Momente belasteter Kreisringsektor. 


Abb, 4. Anderung der Biegespannung am Rand. 
grenzungen a und b, der an beiden Enden durch Biegemomente M belastet ist. Die tberall 


konstante, gleich Eins gesetzte Dicke senkrecht zur Zeichenebene sei klein gegentiber den iibrigen 
Abmessungen, so daB wir den Spannungszustand als einen ebenen auffassen dirfen. 
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Wir fiihren ebene Polarkoordinaten r, y ein, indem wir den Mittelpunkt 0 zum Koordinaten- 
ursprung nehmen. Da in jedem Querschnitt dasselbe Biegemoment ubertragen wird, ist dieser 
Spannungszustand von Winkel y unabhangig. Das Problem hingt also nur vom Radius r ab. Die 
Spannungskomponenten sind dann bekanntlich durch die Beziehungen gegeben: 

1 dF _ @F Ae 
LE Slee res Oo ee? Tro =O. ( ) 


Damit kann die Differentialgleichung (2) in der folgenden Form angeschrieben werden: 


oi ee Se OB! 1 @F 1 dF 1 1 d /dF dg(t)\ _ 
(a3 ry ir) ox oo) Peca) (get eal) are ole dr )=0 (ig 


mit den Abkirzungen 


se : 
0 dr? ' r dr 


9 2 [/d?F\2 1 dF dF 1 /dF\2 
0 = 9G (a=) “or dr? dr > \\ GP) || ~ 
Beachtet man noch (4), so kommt nach kurzer Umformung 


(2A? P’2—2Ar FP! FA F247) 2 F442 7 F’—F) 8 F’2424(27 FF) FF’ 
+2 F"4 (APF? — 3171 PF" F421 F2—1) F’4+rF’ =0, (19) 


wobei wieder die Stoffzahl 


1 (@F , 1dF 
ral 


ad & K 
A=GKLOG (20) 
eingefiihrt wird. Striche bedeuten Ableitungen nach r. 
Dazu treten die folgenden Randbedingungen: 
a= ie Pp Sw mel fale (21) 
b b 
fo,dr=0 und fo,rdr=—M. (22) 


Mit Hilfe von (16) driicken wir die Bedingungen (22) in der Spannungsfunktion aus. Wir bilden 


zunachst 


b b 
ff “er dF \b 
| oar = | Gade = (Fr) =0. (23) 
Durch partielle Integration der zweiten Gleichung (22) erhalten wir 
b b b 
‘ wdzE dF |b “dF 
[corde = [Garde lary [ae 
und wegen (21) folgt 
FP =M. (24) 


Die nichtlineare Differentialgleichung (19) vierter Ordnung 1a8t sich einmal integrieren. Man 
findet 


(2 Ar? B28 — 2 Aer BE BY A 8) Be Fk eh ee ee (25) 


mit einer willkirlichen Integrationskonstanten C, wie man leicht durch Einsetzen nachprifen kann. 
Als Naherungsverfahren zur weiteren Integration der Differentialgleichung (25) wenden wir 


das Galerkinsche Verfahren an. Fiir die erste Ableitung F’ der Spannungsfunktion machen wir den 
sehr einfachen, zweigliedrigen Ansatz 


¥ =(r—b)(r—a) (q+ a0), (26) 
wo c, und cy noch zu bestimmende Freiwerte sind. Jede der zwei sogenannten Koordinatenfunk- 


tionen befriedigt nur die homogenen Randbedingungen (21) und (23) 
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Setzt man die Naherungsfunktion (26) in (25) ein, so erhalt man 
D(r; ¢, 2) =A [12 7% — 11 (6 +a)r4+ (3027+ 2ba4- 34°) r3 + 2ba(b +a) r? 
—ba(b + a)*r-+ b? a? (b + a)] c3 
+ A[71 r§ — 82 (b + a) r® + (30 b? + 53 b a + 30 a?) r4— (3 b§ + 5 b?a+5 ba?+ 34%) r3 
—ba(2b?+3ba+2a*)r? + 2b? a2 (b+ a)r—ba®] ec, 
+ A[135 r§ — 186 (b + a) r> + (82 b? + 189 b a + 82 a?) r4 
— (11 b? + 53 b 2a + 53 b a? + 11 a) 
+ba(2b®+5ba + 2a?) 24 2b a2 (b + a)r—Bb a] rec 
+ A [84 r> — 135 (b + a) r4 + (71 D2 + 1866 a 
+ 71 a®) r3 — 2(6 63 + 41 ba + 41 ba? + 63) r? 
+ ba (116? + 30ba 4 lla?) r—3 Ba (b+ a)] P3 
+ (3 r?— ba)re, + (8r—3b—3 a) rc, — Cr’. 
Nach der Galerkinschen Vorschrift verlangt man dann die Orthogonalitat dieses Ausdrucks zu den 
beiden Koordinatenfunktionen 


b 
f D(r3 ¢,, ey) (r —b) (r —a) dr = 0, 
‘ (27) 


Oe 


D(r3 ¢,, €3) (7 — b) (r — a) rdr = 0. 


Die Auswertung ergibt 


A(Og9 CE + Moy Cf Cg + yy Cy C3 + og 2) +O Cy +O, ¢,—aC =D, (28) 
A(Bao ct + Bar Ct Cy + Bie C3 + Bosc) + Bier + Be e2—fByp C = 9 | 
mit den Abkiirzungen 
Xo = 132 (o + 1) (e— 1)? [11 (e? + 1)— 8] a’, 
Boy = 11 (o— 1)? [343 (9% +- 1) 116 g (e? + 1) + 90-67] a*; 
Oy. = 132 (9 + 1) (oe —1)? [25 (e* + 1) — 149 (0? +1) + 200?) a, 
os = (oe — 1)? [983 (9% + 1) + 493 @ (e* +1) + 269 9? (0? + 1) + 206 9?) a®, 
a — 924 (9 +1) [3 (e* + 1)—e]@, 
OX = 66 [38 (ot + 1) + 230 (0? + 1) + 180? ] a’, 
A = 462 [3 (0? +1) + 4¢] 


und 
Bao = 33 (0 — 1)? [35 (ot +1) + 200 (0? +1) + 2g%] a8, 
Bor = 396 (9 + 1) (0 —1)? [8 (ot + 1) —4.0 (o? + 1) + 6%] a, 
Bie = (@ — 1)? [2857 (o° + 1) + 1547 0 (04 + 1) + 835 0? (0? + 1) + 610°] a5, 
Bos = 4 (o + 1) (0 — 1)? [217 (08 + 1) — 920 (0% + 1) + 164 9? (0? +1) — 116 g° Ja, 
B, = 66 [30 (o8 +1) +270 (¢? +1) + 26¢%]a, 
Bz = 132 (9 + 1) [15 (gt + 1) — 49 (0? + 1) + 139] a’, 
Py = 462 (0 +1) [2 (2 +1) +01, 
wobei das Halbmesserverhaltnis b 
9 == (29) 


gesetzt ist. 
Eliminiert man aus (28) die Integrationskonstante C, so entsteht nach kurzer Rechnung 


A (530 Cf + O21 ct Cy + O12 cy’c¢ + dog 3) + 0, + bo Cg = 0 (30) 
mit den weiteren Abkirzungen 

dg9 = 33 (9 — 1)? [17 (of +1) + 78¢@ (o? + 1) + 138 9") a? , 

dy, = 11 (2 + 1) (e — 1)? [178 (o* + 1) + 501 9 (@? + 1) + 778 9] a*, 

b19 = (o — 1)? [1971 (e8 + 1) + 7207 @ (o4 + 1) + 13867 0? (o? + 1) 4 15422 Oars 

dos = (0 + 1) (0 — 1)? [638 (08 + 1) + 1675 0 (9* + 1) + 2815 979 + 1) + 2416 9? ] a®, 

6, = 66 [6 (ot + 1) + 31 @ (oe? +1) + 660°] a, 

dy = 66 (9 + 1) [14 (ot + 1) + 400 (0? +1) + 6707] @. 
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Nun haben wir noch die Bedingung (24) zu erfiillen. Sie verlangt 


b 
M =f ¥dr=—5 (0-19 2q+et)aal@ 
oder 
P 12 M 
iN a aE ae ac pe ek (31) 


Hat man die Beiwerte c, und c, aus (30) und (31) bestimmt, so kénnen die Naherungswerte fiir 
die Spannungen gemaf (16) berechnet werden. Fir die Spannungskomponenten ergibt sich dann 


m= (r— 0) (ra) (Test ey), (32) 


Oy, = (2r—b—a)c, + [3 r?—2(b+a)r+ba]c, 
und somit fiir die Randwerte der Tangentialspannung 
(09)r=a = — (0 — 1) (¢, a + €, a°), | 
(Gp)r=0 = (0 — 1) (¢, a + @ cy a”). 


Zum Vergleich sei hier noch die bekannte Lésung! der linearen Theorie (A = 0) hinzugefiigt. 
Op, Or Sie lautet 


(33) 


px EM) re ee aaa ae 

He of = — Ay (Ing + atin ® +S ice 

g/em (34) 
2 

7000 Hee ee ote nz + ot — at) 


mit 
b \2 
N = (b? — a*)? — 4B? «(In 2) : 


Die Verteilung der beiden Spannungen ist lings eines Halbmessers 
in Abb.6 dargestellt (und zwar fir 9 =3, M/a? = 200 kg/em und 
A= 5-10 cm*/kg?). Die Lésungen der linearen Theorie nach (34) sind ge- 
strichelt eingezeichnet. Man erkennt eine starke Abnahme der Tangential- 
spannung og, am Innenrand, wahrend die radiale Spannung o, gegeniiber 
der linearen Theorie kaum geandert wird. 


SchlieBlich soll noch bemerkt werden, dah 
wir eine gréBere Genauigkeit fiir die Nahe- 
rungswerte der Spannungsverteilung durch 
/ einen Ansatz erhalten kénnen, in dem die 


500 |- 


gen fiir sich erfiillen (wie es ja beim Galerkin- 
schen Verfahren tatsachlich sein muf). Die hier 
gezeigte Rechnung hat mehr qualitative Bedeu- 


des Vorganges bei einem nichtlinearen Elastizi- 
tatsgesetz wieder; doch ist nach Erfahrungen 
bei anderen Problemen dieser Art anzunehmen, 
daB auch zahlenmafSig den Ergebnissen kein 
Abb. 6. Spannungen im Kreisringsektor. allzugroBer Fehler anhaftet. 
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(Eingegangen am 31. Januar 1955.) 


Anschrift des Verfassers: Dr.-Ing. F. Jindra, Stuttgart-O, Gerokstr. 70. 


' Vel. z. B. S. Timoshenko, Theory of Elasticity, 1. Aufl., S. 58, New York 1934. 


Koordinatenfunktionen alle vier Randbedingun- — 


tung; also gibt sie nur noch das Grundsatzliche _ 
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| Stahileichtbau von Maschinen 

; Von Dipl.-Ing. K. Bobek, Konstruktionsdirektor der Allg. Elektrizitits-Gesellschaft Berlin, Dr.-Ing. 
; A. Hei$, Alfred Treves, Maschinen- und Armaturenfabrik KG., Frankfurt/Main, und Dr.-Ing. 
; _ Fr. Schmidt, Oberingenieur der MAN Augsburg. (Konstruktionsbiicher. Rise. Prof. Dr.-Ing. K. Koll- 
: mann, Karlsruhe, 1. Band.) Zweite neubearbeitete Auflage. Mit 243 Abbildungen. VII, 
q 183 Seiten Gr.-8°. 1955. Ganzleinen DM 24,— 
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Inhaltsibersicht: I. Grundlagen des Stahlieichtbaues. Von Dipl.-Ing. K. Bobek, Berlin: Leichtbau allgemein 
Vorteile des Stahlleichtbaues. Allgemeine Konstruktionsgrundsatze. Axioitavertikcon: Grundlagen var Taaughatdecinenes 
a Ti. Stahlleichtbau von Elektromaschinen. Von Dipl.-Ing. K. Bobek, Berlin: Die besonderen Bedingungen und Grund- 
saize des Elektromaschinenbaues. Kleinmaschinen. Grofmaschinen. — Il. StahIschweiSbau von Werkzeugmaschinen, Von 
Dr.-Ing. A. Hei8, Frankfurt/Main: Die Entwicklung des StahlschweiSbaues von Werkzeugmaschinen. .Grundformen von 
Werkzeugmaschinen-Gestellen. Entwicklungsrichtungen fiir die Gestaltung yon Werkzeugmaschinengestellen in Stahlschweib- 
bau. Bemessungforderungen im Werkzeugmaschinenbau. Das SchwingungsmeGprinzip nach Kettner. Elementegestaltung 
im. StahlschweiBbau. Erklarung der VergréS¢rung der Dampfung durch Scheuerwirkung. *Konstruktionsbeispiel aus dem 
Zellenbau. Konstruktionsbeispiele aus der Plattenbauweise. Kurze Zusammenfassung iiber den Stand der Anwendung des | 
StahischweiSbaues und Ausblick auf die zukiinftige Entwicklung. — IV. Stahlleichthau fiir Verbrennungemaschinen. Von 
Dr.-Ing. Fr. Schmidt, Augsburg: Allgemeines iiber Stahlleichtbau fiir Verbrennungsmaschinen. Richtlinien fir die 
Kenstruktion. Entwurf yon Gestellen und Grundplatten. Weitere Bauteile yon Verbrennungsmotoren. Allgemeine Richt- 
. linien-fir den Leichtbau. — Schrifttum zu I—IV — Sachvyerzeichnis. 
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Gestaltung von Fahrzeuggetrieben 
Von Dr.-Ing. Hans Reichenbicher, Osnabriick. (Konstruktionsbiicher. Hrsg. Prof. Dr.-Ing. K. Koil- 
mann, Karlsruhe, 15, Band.) Mit 164 Abbildungen im Text und auf 4 Tafeln. VII, 155 Seiten 
Gz.-8°.. 1955. Ganzleinen DM 26,40 


Inhaltstbersicht: Umfang der Darstellung. — Grundforment der Getriebe. — Wahl der Getriebestufen. — Ubersetzung 

. des gréBten Fahrganges. Der Schnellgang. Der kleinste Gang: Geschmeidigkeit. Steigfahigkeit. Beschleunigungsvermégen. — 
Der Riickwartsgang. — Der Zwischengang: Geometrische Stufung. Auslegung nach den Geschmeidigkeitsschaubildern. Aus- 
legung nach dem Beschleunigungsvermégen. Zusammenfassung zur Zwischengang-Ubersetzung. — Reine Zahnradgetriebe. — 
Standgetriebe: Entwicklung. Bedienung der Schaltung. Unmittelbar von Hand geschaltete Getriebe. Schalterleichterungen. 
FuBschaltung. — Umlaufgetriebe: Momente, Drehzahlen, Wirkungsgrad. Formschliissige Kupplung und Freilauf. Kraft- 
-schliissige Bremsen und Kupplungen. Uberblick iiber die Umlaufgetriebe, Bandbremsen. — Einige Einzeltcile und ihre Be- 
rechnung: Zahnrader. Lager. Wellen. Naben. — Féttinger-Kupplung. — Grundgesetze der. Fottinger-Kupplung. — Aus- 
-legung, Bauarten: Normale Ausfithrung. Stérglieder-Teilentleerung. Lage der Fottinger-Kupplung in Getriebe. Leistungs- 
verzweigung, Hydra-Matic-Getriebe. — Einfluf der Féttinger-Kupplung auf die Fahrleistungen. — Vergleich Strémungs- mit 
Reibungskupplung. — Strémungswandler. — Grundgesetze des Strémungswandlers: Bauarten «der Strémungswandler. 
Hinfacher Wandler. Einkreislaufer, dreiteilig. Einkreislaufer, vierteilig. Einkreislaufer, finfteilig. Einkreislaufer mit 
Leistungsteilung. Radialwandler. — Einflu8 des Strémungswandlers auf die Fahrleistungen: Zwei Grundcharakteristiken 
des Strémungswandlers. Vergleich zwischen Stsimungswandler, Féttingerkupplung und reinem Zahnradgetriebe in bezug 
auf Steigfahigkeit und Beschleunigungsvermégen. Vergleich in bezug auf Brenustoffverbrauch bei Vollast. Brennstoffver- 
brauch bei Teillast. Zusammenfassung. — Schrifttum. Sachverzeichnis. 
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Technische Kompositionslehre 


Anleitung zu technisch-wirtschaftlichem und verantwortungsbewubtem Schaffen. Von Dr. sc. techn. 
Fritz Kesselring, Zollikon (Schweiz). Mit 162 Abbildungen. IX, 394 Seiten Gr.-8°. 1954. 
: ¢ Ganzleinen DM 39,— 


Inhaltstibersicht: Einfithrende Gedanken. Erster Teil: Die allgemeinen Voraussetzung én 
des technischen Schaffens. 1. Die ethischen Voraussetzungen des technischen Schaffens: Die Philosophie und 
ihre Bedeutung fiir eine Sinngebung der Technik. Von der Kunst als dem grofen Ausgleich. Uber die Lebensgestaltung des 
Ingenieurs. Das Verantwortungsproblem der Technik. II. Die wissenschaftlichen Voraussetzungen des technischen Schaf- 
fens: Ziele und Methoden der Wissenschaft. Von der Schénheit und dem Nutzen der Mathematik. Die Bedeutung von Wahr- 
scheinlichkeitsrechnung und Statistik, Einige grundsatzliche physikelische Zusammenhdnge. — III. Die wirtschaftlichen 
Voraussetzungen des technischen Schaffens: Vom Wesen und der Bedeutung wirtschaftlicher Betrachtungen. Die Herstell- 
kosten bestehender technischer Erzeugnisse. Die Herstellkosten in Entwicklung befindlicher Erzeugnisse. Das Problem der 
Rationalisierung. — Zweiter Teil: Die Grundlehren des technischen Schatfens, — IV. Er- 
findungslehre: Die Geschichte einiger Erfindungen. Der Weg zur Erfindung. Der Weg zum Patent, — V. Cestaltungslehre : 
Hinige tibergeordnete Gestaltungsprinzipien. Die Gestaltung mit Hilfe eines konvergierenden Naherungsverfahrens. Die Ge- 
staltung mit Hilfe mathematischer Methoden. — VI. Formungslehre: Der Stoff und seine Veranderung durch Felder, Krafte 
und Impulse. Die physikalischen Grundlagen der Formung, Die technisch-wirtschaftliche Formung von Einzelteilen und 
ihrer Verbindung. Anhang: Erlauterungen zu den Tafeln I bis VI. Sachverzeichnis. 
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Die Grundgesetze der Warmeiibertragung 
Dritte, véllig neubearbeitete Auflage von Dr.-Ing. Ulrich Grigull, Farbenfabriken Bayer A. e 
Leverkusen. Mit 190 Abbildungen. VIII, 428 Seiten Gr.-8°. 1955. Ganzleinen DM 37,50 


Inhaltsiitbersicht: Erster Teil: Wiirmeleitung in festen Kérpern. A. Die mathematischen Grundlagen. ee 

B. Uber die Lésung von Randwertaufgaben. — C, Die zeitlich veranderlichen Tem peraturfelder ohne Warmogaplien, rE D. Die - 
zeitlich konstanten Temperaturfelder ohne Warmequellen. — E, Die zeitlich konstanten Temperaturfelder mit Warme- 
quellen. — F. Verschiedene Sonderfille. — G. Warmeleitung in verdiinnten Gasen,— Zweiter Teil: Konyektive Warme- 
ubertragung. A. Fliissigkeits- und Encrgiebewegung. — B. Das Ahnlichkeitsgesetz der Warmeiibertragung. — Cc, Warme- 
iibergang bei erzwungener laminarer Strémung. — D. Warmeiibergang bei erzwungener turbulenter Strémung. — E. Warme- 
iibergang bei querangestrémten Kérpern.. — ¥. Warmeiibergang hei hohen Geschwindigkeiten. — G. Warmetibergang bei 
freier Konyektion; — H. Warmetibergang bei Kondensation. — I. Warmetibergang bei Verdampfung. — K. Stoffiibertragung. 

—Dritter Teil: Wirmestrahlung. A. Einfihrung. — B. Die Strahlung des schwarzen Kérpers. — C. Der Strahlungs- 
austausch zwischen festen Kérpern. — D. Strahlung von Gasen und Dampfen. — Anshan. Einheiten, Umrechnungstafeln, 
Stoffwerte. — Schrifttum. — Namen- und Sachverzeichnis. i 


Der rasche Fortschritt auf dem Gebiet der Warmeiibertragung erforderte einé durchgreifende Neubearbeitung des in zweiter. 
Auflage 1933 von Gréber und Erk herausgegebenen Buches. Der Teil ,,Konvektive Warmeiibertragung‘‘ wurde véllig neu ge- 
schrieben. Aber auch die beiden anderen Teile des Buches, die Warmeleitung und die Warmestrahlung, wurden an vielen Stellen 


~~ om 


ergiinzt und auf den neuesten Stand gebracht. : z 


Die Lehre von der konyektiven Warmeiibertragung verdankt ihren Fortschritt im wesentlichen zwei Methoden: der An 
wendung der Ahnlichkeitstheorie undder Betrachtung der wandnahen Schichten als einen besonderen Teil der Strémung, also. 
der Prandtlschen Grenzschichtlehre. Von beiden Methoden wurde in der neuen Auflage ausgiebig Gebrauch gemacht, die 
Ahnlichkeitstheorie wurde besonders ausfiihrlich behandelt. Daneben wurde die~Analogie zwischen Impuls-, Warme- und 
Stoffaustausch zu einer gemeinsamen Darstellung der drei Vorginge benutzt, uni dadurch besonders die Austauschvorgange 
in turbulenter Strémung berechnen zu kénnen. Das Buch gibt damit feinen Uberblick tiber die heute tiblichen Lésungsmethoden 
und zeigt die erzielten Ergebnisse der Theorie. Die Verbindung zur praktischen Anwendung liefern Abschnitte tiber Versuchs- 
ergebnisse und Gebrauchsformeln. Gerade bei der Warmetibertragung erfordert die richtige Anwendung einer Gleichung ein 
gewisses MaB von physikalischem Verstandnis; daher wurde bei der Darstellung hierauf besonderer Wert gelegt. Cemaf ihrem 
Titel beschrankt sich auch die 3. Auflage auf die ,,Grundgesetze“‘, da nur auf diese Weise das Gesamtgebiet in einem ene be- 
handelt werden konnte. 


Das Buch wendet sich in erster Linie an den Ingenieur, und zwar ebenso an den Praktiker wie an den ie tie cain ar- 
beitenden. Auch der Student der héheren Semester wird das cate als Studienhilfe heranziehen kénnen. ~ 2 BS 
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Das Cross-Verfahren. 


Die Berechnung biegefester Tragwerke nach der Methode des Motcenénaucgletthe: Von Dr.-Ing, 
Johannes Johannson. Zweite , vollig umgearbeitete Auflage von Dipl.-Ing. Gin’ ter Raczat. ‘Mit 
161 Abbildungen und 38 Becbepbeipicles: XII, 168 Seiten Gr.-8°. 1955. DM 17,40 — 


Inhaltsiibersicht: Grundlagen: Phystkallsche Deutung ie Cross-Verfahrens, Aufgabenbereich. Aligeugeine Verein- 
barungen fiir Berechnungsform und Darstellung. Die wichtigsten Zeichen. Rechnerische Verfolgung der mechanischen Vor- 
ginge. Steifigkeiten und Ubertragungszahlen. Verteilungs- oder Ausgleichzahlen. Zusammengesetzte Drehsteifigkeiten. 
Teilverformungsbilder. — Die allgemcine Handhabung des Cross-Verfahrens bei Durchlauftragern und einfachen Rahmen- 
werken mit und ohne Stabdrehungen unter auSeren Lasten: Allgemeine Regeln. Beispiele fur einmaligen Momentenausgleich. 
Beispiel fiir die Anwendung der Cross’schen Iteration in der augfiihrlichen Form, Kiirzere Schreibweise der Cross’schen Iteration. 
Beispiele fiir Rahmenwerke mit Stabdrehwinkeln (verschiebbaren Knoten). — Einige Verfahren zur Abkiirzang des Berech- es | 3 
nungsganges bei Durchlauftrigern und Rahmenwerken ohne Stabdrehwinkel: Allgemeines. Das zweistufige Verfahren nach ; ee 
Dasek. Umgehung der Cross’schen Iteration. — Behandlung von Rahmenwerken mit ein- oder mebrfacher Symmetrie: Be- oh aH 
lastungsumordnung. Ansatz der Steifigkeiten, wenn Stabdrehwinkel nicht vorhanden sind. Ansatz der Steifigkeiten, wenn _ 
Stabdrehwinkel vorhanden sind. Regeln fiir die Zusammensetzung der Ergebnisse bei unsymmetrischer Belastung. Beispiele. - ; 
— Ein besonderes Verfahren fiir eingeschossige Rahmenwerke: Allgemeines. Beispiele. — Die Handhabung des Cross-Ver- 
fahrens bei groferen Stockwerkrahmen, Rahmentrigern und bei Durchlauftragern auf élastisch senkbaren Stiitzen: Allgemeines. 
Beispiele. Durchlauftriiger auf elastisch senkbaren Stiitzen. Durchiauftrager mit gegebenen Auflagesenkungen, Neben-— 
spannungen in Fachwerken. — Einflu8 von Temperaturinderungen, Schwinden und Normalkraften: GleichmaBige Temperatur- 
inderung. UngleichmaSige’ Temperaturainderung. EinfluB von Langeniénderungen durch; Normalkrafte. Beispiel. Einflu8- 
linien: Allgemcines, Beispiel. — Verdnderliche Triigheitsmomente: Ausgangsmomente, Steifigkeiten, Ubertragungs- fund 
Ausgleichszahlen. Beispicle. — Anhang: Tabellen I-4X. — Namen--und Sachyerzeichnis.: 
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